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A reasoning mind—man has always found 

- this his most useful weapon when venturing 

_ into unknown fields. Today’s frontier—one 
of the greatest unknowns ever faced by man 
—lies in outer space. 

RMI, America’s first rocket family, is 
uniquely qualified to help our country con- 
quer the new frontier. For this project will 

require more powerful and more efficient 
rocket propulsion systems than any yet de- 
signed. And RMI’s experience in the rocket 
_ power field extends back over 15 years... 
- years in which RMI scientists and engineers 
poo‘ed their reasoning power and knowledge 


Power for Progress 
® 


to develop entirely new methods of propul- 
sion and to produce rocket powerplants for 
dozens of military vehicles. 

With this backlog of experience—and with 
an increasing staff of imaginative, reasoning 
engineers—RMI continues to blaze trails in 
many fields of rocket power: supersonic 
manned aircraft ... versatile helicopters... 
catapults and test sleds . . . and missiles for 


defense and space exploration. 

Engineers, Scientists—Perhaps you, too, 
can work with America’s first rocket family. 
You’ll find the problems challenging, the 
rewards great. 
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FASTEST—Lino-Writ 4 is the fastes )- 
graphic paper you can use—and the paper 
with the widest range. As shown by this 
photograph, you can record frequencies from 
60 to 2000 cps on the same roll. 


TOUGHEST—Snap 


SPLICE-FREE TOO!— Whatever speed you require, what- 
ever thickness is best for you, whatever length roll you 
can use... you can get Lino-Writ splice-free, at no extra 
cost. Lino-Writ | for low frequency recording. Lino- 
Writ 2 for mid-range recording. Lino-Writ 3 for high 


Du Pont Lino-Writ 4 springs back like new 
... with its valuable record intact! Its all-rag 
base not only makes it rugged, but also makes 
it whitest, easiest to read. 


THINNEST— You get more paper on a standard 
roll. This roll holds 475 feet of Lino-Writ 4. 
Roll with same outside diameter holds only 
250 feet of standard-weight paper. Translu- 
gence permits quick duplication in standard 
equipment. 


frequency recording. New Lino-Writ 4, the fastest, 
toughest, thinnest photorecording paper you can use. 

E. I. du Pont de Nemours & Co. (Inc.), Photo 
Products Dept., Wilmington 98, Delaware. In Canada: 
Du Pont Company of Canada (1956) Limited, Toronto. 


REG. U.S, pat. OFF. 


Better Things for Better Living 
... through Chemistry 


DU PONT OSCILLOGRAPHIC PRODUCTS 
for Functional Photography 


Photography with a purpose...not an end in itself, but a means to an end. 


DECEMBER 1957 


1219 


| 
| 
4 600 2000 
2 
j 
| 
| 
| 
! 
| 
| 
+ 
M 


publication of the 


AMERICAN ROCKET 


Rescarch and Development aid 


IRWIN HERSEY—DIRECTOR OF PUBLICATIONS 


EDITOR 
ASSISTANT EDITOR 


ART EDITOR 


MARTIN SUMMERFIELD 
STANLEY BEITLER 
JOHN CULIN 


ASSOCIATE EDITORS © 


ALI BULENT CAMBEL, Northwestern University 
IRVIN GLASSMAN, Princeton University rat 
M. H. SMITH, Princeton University a 


CONTRIBUTORS 


MARSHALL FISHER, Princeton University 
GEORGE F. McLAUGHLIN 


ADVERTISING & PROMOTION MANAGER 
WILLIAM CHENOWETH 


ADVERTISING 


D.C. Emery & Associates James C. Galloway & Co. 
155 East 42 St., New York, N. Y. 6535 Wilshire Blvd., Los Angeles, Calif. 
Telephone: Yukon 6-6855 Telephone: Olive 3-3223 


_ Jim Summers & Asssociates R. F. and Larry Pickrell 
35 E. Wacker Dr., Chicago, Ill. 318 Stephenson Bidg., Detroit, Mich. 
rae —_ Telephone: Andover 3-1154 Telephone: Trinity 1-0790 


Harold Short ak Louis J. Bresnick 
Holt Rd., Andover, Mass. mi 304 Washington Ave., Chelsea 50, Mass. 
Telephone: Andover 2212 Telephone: Chelsea 3-3335 


AMERICAN ROCKET SOCIETY 


> 


OFFICERS 

Robert C. Truax 


President 
Vice-President George P. Sutton | 
Executive Secretary aie James J. Harford 
Secretary A. C. Slade 
Treasurer Robert M. Lawrence 
General Counsel Andrew G. Haley 

Terms expiring on dates indicated _ 7 
Krafft Ehricke, 1959 ; Milton W. Rosen, 1957 
Andrew G. Haley, 1957 a H. S. Seifert, 1958 
S. K. Hoffman, 1958 : John P. Stapp, 1959 
H. W. Ritchey, 1959 K. R. Stehling, 1958 


Wernher von Braun, 1957 


TECHNICAL DIVISION CHAIRMEN 


David G. Simons, Human Factors John F. Tormey, Propellants and 


John J. Burke, Instrumentation Combustion 
and Guidance Brooks T. Morris, Ramjet 
Edward N. Hall, Liquid Rocket William L. Rogers, Solid Rocket 


Krafft A. Ehricke, Space Flight 


Scope of JET PROPULSION 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance- 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. As 
used herein, the term “jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct, and thus it includes air-consuming 
engines and underwater systems as well as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other aspects of jet-propelled flight, 
such as flight nani, guidance, telemeter- 
ing, and research instrumentation. In- 
creasing emphasis will be given to the scientific 
problems of extraterrestrial flight. 


Information for Authors 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, 
typical maximum length is 22 pages of text 
(including equations), 1 page of references 
1 page of abstract, and 12 illustrations. 
Fewer illustrations permit more text. and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes. They may be devoted 
either to new developments requiring prompt 
disclosure or to comments on previously pub 
lished papers. Such manuscripts are usual! 
published within two months of the date of 
receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal! 
or as a special supplement if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Do not type equations; write 
them in ink. Identify unusual symbols or 
Greek letters for the printer. References are 
to be grouped at the end of the manuscript 
and are to be given as follows: for journal 
articles: authors first, then title, journal, 
volume, year, page numbers; for books: 
authors first, then title, publisher, city, edi- 
tion, and page or chapter numbers. Line 
drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 


JET PROPULSION is published monthly by 
the American Rocket Society, Inc., and the 
American Interplanetary Society at 20th & 
Northampton Sts., Easton, Pa., U. 8S. A. 
Editorial offices: 500 Fifth Ave., New York 
36, N. Y. Price: $12.50 per year, $2.00 
per single copy. Entered as second-class 
matter at the Post Office at Easton, Pa. 
Notice of change of — should be sent to 
the Secretary, ARS, at least 30 days prior to 
publication. Opinions expressed herein are 
the authors’ and do not necessarily reflect the 
views of the Editors or of the Society. 
© Copyright 1957 by the American Rocket 
Society, Inc. 
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This is 


NATIONAL 


4 
NO RT H E; RN Test a rocket motor’s performance at 100,000 feet altitude... deter- 
mine how extremes of temperature affect explosive and propellant 
charges...develop igniters of maximum sturdiness and reliability. 
ase These are typical activities at National Northern Division, West Han- 
newest APé& Ps division over, Mass., latest addition to the AMERICAN POTASH family. Here, a 
eee qualified staff of scientists and technicians conducts research, devel- 


opment, field testing and production of propellants for rockets and 
missiles, explosives, pyrotechnics, detonators, igniters, squibs, gas 


Dedicated to research, generators, fuzes and other items. 
° ° Extreme environmental conditions are simulated in the laboratory 
development and production mM or encountered in actual field tests on the firing range at West Han- 
over or the 1900-acre test range at Halifax, Mass. Here, too, AP&CC 
ro cket ower will extend its facilities for the investigation and development of 
p TRONA* boron, lithium and perchlorate chemicals in high energy fuels. 
For more than ten years National Northern has demonstrated its 
ability to work in close coordination with others. As an integral part 


of AMERICAN POTASH & CHEMICAL CORPORATION it is ready to go to 
work for you! 


We invite inquiries leading to research and production contracts. 


NATIONAL NORTHERN’S CURRENT FIELDS OF ACTIVITY: 


ROCKET-MOTOR-IGNITER DEVELOPMENT e GAS GENERATORS e SURVEILLANCE TESTING OF 
EXPLOSIVES AND AMMUNITION . EFFECT OF ALTITUDE ON PERFORMANCE OF EXPLOSIVES, 
PROPELLANTS, TRACERS AND INCENDIARY ELEMENTS e¢ DESIGN OF ELECTRIC DETONATORS AND 
th cot e@ AIRCRAFT VULNERABILITY TESTS AND EVALUATION e EXPLOSIVE METAL FORMING 


& Potassium Chlorates 


American Potash & Chemical Corporation 


Nitrate 
Perchlorate 


3030 WEST SIXTH STREET, LOS ANGELES 54, CALIFORNIA 
99 PARK AVENUE, NEW YORK 16, NEW YORK 
NATIONAL NORTHERN DIVISION 
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Where Miles are Measured 
in Thousandths of an Inch 


As ballistic missiles go farther and faster, we becomes an 
increasingly critical problem. Walls of rocket cases, for example, must 
be still stronger to withstand higher pressures, yet still thinner to cut 
‘hardware’ poundage to a minimum. 

Obtaining optimum wall thinness throughout im rockets now being 
developed is one of M. W. Kellogg's current defense assignments. It is 
a complex engineering task and demands a highly specialized 
knowledge of metallurgy, metal fatigue, heat transfer, corrosion, 
weight-strength relationships, and welding techniques. 

The M. W. Kellogg Company has been closely associated with the 
development and fabrication of rocket cases and Components since 
1947. Inquiries are invited from organizations wishing to put Kellogeg’s 


Specialized background to work on specific PrOBIEMS. 


DEFENSE PRODUCTS DIVISION 
THE M. W. KELLOGG COMPANY 


711 THIRD AVENUE, NEW YORK 17, N. Y. 
Lan A SUBSIDIARY OF PULLMAN INCORPORATED 


The Canadian Kellogg Company Limited, Toronto « Kellogg International Corporation, London 
Companhia Kellogg Brasileira, Rio de Janeiro « Compania Kellogg de Venezuela, Caracas 
Kellogg Pan American Corporation, New York « Societe Kellogg, Paris 
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em Petroleum 


by Phillips 

_ Phillips solid propellants for rockets, made from low 
cost petrochemicals, come in a wide range of im- 
pulses for a variety of applications where required. 
_ These new solid propellants can be made relatively 
smokeless. They operate successfully at tempera- 


by impact or explosion, they are also capable of be- 
ing stored for long periods without deterioration. 
Exhaust gases are noncorrosive and relatively low in 
temperature. The facilities of the Phillips operated 
_ Air Force Plant 66, near McGregor, Texas, are com- 
_ pletely equipped for development, testing and man- 
_ ufacture of rocket motors. Consult Phillips about 
your propulsion needs. 


We invite you to write and discuss your 
technical and production problems with us. 
Or call one of our regional representatives: 
Washington, D. C.—E. L. Klein, EXecu- 
tive 3-3050; Los Angeles—R. O. Gose, 
GRanite 2-0218; Dayton, Ohio—C. W. 
Strayer, YOrktown 3263. 


PHILLIPS 
PETROLEUM COMPANY 


Bartlesville, Oklahoma | 


Address all inquiries to: 


Rocket Fuels Division, Bartlesville, Oklahoma 


| 
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Jet-Age Strateg’ 


PUMPS 


provide dependable 
high-speed fuel flow 


Today, in-flight refueling is the key to America’s 
global air defense — saves time that can speed a 
plane to any possible target with the Sunday punch 
of retaliation, in the event of sneak attack. Where 
seconds can make the difference, America’s 
fighting aircraft are being air-refueled in 
increasing numbers. At the Propulsion 
Research Division of Curtiss-Wright, highly 
trained and experienced specialists are 
developing a family of compact, high 
performance pump assemblies which can 
transfer fuel from tanker to striker aircraft 
at extremely high rates of speed. Units 

are available to meet a variety of fuel flow 
and pressure requirements. Refueling pumps 
are a good example of the contribution 

of the Propulsion Research Division to a 
PROPULSION RESEARCH DIVISION wide variety of accessory power requirements 


C for our air forces, and for industry. 


CORPORATION * SANTA MONICA, CALIFORNIA 
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Now 


As prime producers of unsym-Dimethylhydrazine . . . 
under the tradename DIMAZINE .. . Westvaco has compiled 
and kept up-to-date all the data you need on the storage 
and handling of this advanced rocket fuel. 


The third edition of this useful book is just off the press and © 
available now. For your copy, return the coupon below. 


DIMAZINE is Westvaco’s trademark for unsym-Dimethylhydrazine 


Gentlemen: Please forward a copy of your book, 
‘ Storage and Handling of Dimazine. 


NAME 


TITLE 


COMPANY. 


ADDRESS 


CITY ZONE__ STATE. 


Chior-Alkali Division 
FOOD MACHINERY AND CHEMICAL CORPORATION 


161 East 42nd Street, New York 17, New York 


e 
BECCO® peroxygen chemicais * FAIRFIELD® pesticide compounds « FMC® organic chemicals NIAGARA® fungici ana 
industrial suiphur « OHIO- APEX® piasticizers and resins + WESTVACO® alkalis, Pp Pp barium and 9 
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PRECISION FORGINGS 
WYMAN-GORDON 


LENGTH 29” 


LENGTH 49” 


.-. Close tolerance .. . low-draft ... 
wit ; no-draft — these are relative terms and con- 

siderable confusion has been created because 
they have been loosely used. Close tolerances on large forg- 
ings have been pioneered by Wyman-Gordon over the years. 
The degree of precision which is practical now becomes a 
question of economics. 


The forgings illustrated are typical of Wyman-Gordon close 
tolerance forgings that have been produced in quantity. In 
solving problems of tolerances, design and metallurgical 
quality, there is no substitute for Wyman-Gordon experience. 


WYMAN-GORDON COMPANY 


Established 1883 


FORGINGS OF ALUMINUM @ MAGNESIUM @ STEEL @ TITANIUM 
WORCESTER 1, MASSACHUSETTS 
ILLINOIS DETROIT, MICHIGAN 
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TECHNIQUES ond DEVELOPMENTS 


in oscillographic recording 


CIRCUIT DESIGN AND TYPICAL USES 
OF THE “150” CARRIER PREAMPLIFIER 


One of the most frequently used plug-in front ends for 
Sanborn 150 Series oscillographic recording systems is the 
Model 150-1100 Carrier Preamplifier, since with it a “150” 
system can record such variables as force, temperature, 
strain, pressure, displacement, velocity, flow, acceleration — 
or any variable which can be 
expressed as a suitable input 
signal by a transducer. The “1100 
Carrier” will operate with a 
variety of different transducers 
and bridge circuits, which will be 
mentioned later on. 

In the block diagram (Fig. 1), 


MIXER DEMOD. DRIVER AMP. 


apush-pulloscillator 
provides an excitation 
voltage of about 5 v. 
to the transducer at a 
standard frequency of 
2,400 cycles or optional 
frequencies of 600 and 
1,200 cycles, using 
plug-in components. 
This excitation voltage also feeds the Balancing. Calibration 
and Zero Suppression circuits. (The Balancing controls allow 
correction of resistive and reactive signal leakage from the 


t 
bid 
~ 


em 


transducer, so that at zero load the net signal to the Pre- 
amplifier is zero. The Zero Suppression feature permits 
bucking out a large static load so that a small part of the 
load can be expanded over the full recording chart. The 
Gage Factor control allows the zero suppression range to be 
made equivalent to some convenient transducer load, or the 
full load rating of the transducer, and also causes the cali- 
bration signal to represent 2% of that load.) Transducer 
output is fed to the transformer through the Gage Factor 
potentiometer, across which the Balancing-Calibration-Zero 
Suppression circuits develop a voltage effectively in series 
with the transducer output. The mixer receives a suppressed 
carrier AM signal and re-inserts a carrier component, to 
make its output a conventional AM signal whose modulation 
represents the transducer load. The modulation signal 
(whose amplitude and polarity represent magnitude and 
direction of transducer output) is recovered by the demodu- 
lator and fed to the output amplifier, which in turn excites 
the Driver Amplifier and recording galvanometer of a 
“150” system. 

Transducers which may be used with the Carrier Pre- 
amplifier include strain gage half-bridges or full-bridges, 
commercial resistance or reactance bridges, differential 
transformers and resistance thermometer bridges. The trans- 
ducer chosen should provide at least 18.0 microvolts per volt 
of excitation at the minimum load to be recorded, for a one 
cm. deflection; impedance should be 100 to 1000 ohms. With 
strain gages, normal operation provides sensitivities of 50, 
20 or 10 micro-inches per inch for each cm. on the recording, 
depending on the number of active gages. With resistance 
thermometers, if 1°C. or 2°F. per cm. stylus deflection is 
sufficient sensitivity, the user can construct his own resist- 
ance thermometer by including a 3 ohm coil of copper wire 
in one arm of an equal arm 100 ohm bridge. 


Helpful information about the use of transducers with the 150-1100 Preamplifier is contained in the followi ing 


Le Sanborn RIGHT ANGLE arlicles (reprints on request): Coupling Differential Transformers ug. and 
ie Nor. 1956; Filter Networks for use w uh Force Dvnomometers, Nov 1956: Calibration with 1-, 2 or 4-arm 
: Strain Gage Bridges, Aug. 1955; Theoretical and Actual Applications of Bridge Circuits, May and Aug. 1954 


traffic studies . . 


with the versatile ‘‘1100 Carrier” 


Today, Carrier Preamp-equipped Sanborn “150” systems are being used for 
frequency response tests of process control system components; to record 
shaft deflections of fluid mixing equipment; in infrared research . . . vehicular 
. submarine hull vibration measurements. Applications are 
limited only by the transducers available. 
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These are applications of only one 150" front-end; eleven more inter- 
changeable, plug-in. Preamplifiers increase the scope of Sanborn 
oscillographic recording systems to meet an almost infinite variety of 
research, production and field testing requirements. All Sanborn “150” 
direct writing systems record inkless traces in true rectangular coordi- 
nates; all provide 1% linearity; Basic Assemblies — equipped with your 
choice of Preamps — are available from one- to eight-channels 
packaged in vertical cabinets, portable cases, specially mod 
fled housings. 


SANBORN COMPANY 
INDUSTRIAL DIVISION 
175 Wyman St., Waltham 54, Mass. 


Technical data and help with your 
oscillographic recording problem are 
always available from Sanborn. 
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SERGEANT 
NIKE-HERCULES 
JUPITER JR. 


JUPITER SR. 


(XM~17; XM-20, XM-33) 


POLARIS O 
POLARIS A 


MANY OTHER CLASSIFIED PROJECTS 


A SMALL EXPERIENCED ORGANIZATION GEARED TO 
HANDLE YOUR DEVELOPMENT AND PROTOTYPE REQUIREMENTS 
FOR STATIC AND FLIGHT TESTS IN THE SHORTEST POSS- Ee 
-| TIME. 


CALL OR WRITE 


SILVER ‘CREEK, NEW YORK 


JET PROPULSION 


BUILDERS OF MORE LARGE,THIN WALL, HIGH 
OTHER COMPANY IN AMERICA, 
RV-A-10 
VEHICLES 
I 
— 


ON 


many possible © 


BAFFLE PLATES 
BEARING SURFACES 
CRUCIBLES 


FEED ROLLS 
NOZZLES 
PUMP SHAFTS 
GAS TURBINE PARTS | 
DUCT LININGS 
EXHAUST PORTS 
THERMOCOUPLE TUBES 
THERMOCOUPLE WIRE 
GENERAL WIRE 


applications 


MECHANICAL SEALS 


HATIAVA WOK” 


The three types of ROKIDE spray 
coatings — “A” aluminum oxide, 
“ZS” zirconium silicate and ‘Z”’ zir- 
conium oxide — are hard, adherent, 
crystalline refractory materials. They 
can be applied to underlying parts, 
particularly metals, of a wide variety 
of shapes and sizes. 

Parts on which ROKIDE spray coat- 
ings are especially useful are those re- 
quiring thermal or electrical insula- 


tion ... refractoriness .. . resistance 
to wear or corrosion .. . high melting 
points . . . excellent mechanical 


strength ... dimensional stability . 
relative chemical inertness. 

The coatings are proving their 
great protective values in a rapidly 
growing number of applications. In- 
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cluded are not only the newer fields, 
such as the manufacture of missiles, 
ram jets, nuclear reactors, etc. In- 
dustry in general, including foundries, 
metal melting, heat and power plants 
— plus makers of automotive en- 
gines, electrical and electronic prod- 
ucts, heavy and light chemicals, 
ceramics, and pumps — are also 
safeguarding many different parts 
with these new Norton developments. 

Facilities for applying ROKIDE 
coatings are maintained at Norton 
Company, Worcester, Mass., and at 
its plant 2555 Lafayette Street, Santa 
Clara, California. 

For the latest ROKIDE Bulletin, 
write to NORTON COMPANY, 731 New 
Bond St., Worcester, Massachusetts. 


ROKIDE* Spray Coatings 
Widest Variety of Parts 


Pad ue 


Protect 


Applying ROKIDE Coating to pump shaft. 
- 


WNORTONY 


NEW PRODUCTS 


Galaking better products...» 
to make your products better 


NORTON PRODUCTS 
Abrasives « Grinding Wheels 
Grinding Machines « Refractories 


BEHR-MANNING DIVISION 
Coated Abrasives *« Sharpening Stones 
Behr-cat Tapes 


*Trade-Mark Reg. U. S. Pat. Off. and Foreign Countries 
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HYDROSPINNING 


Big Nozzles for Big Rockets Now avaiaste 


A new Hydrospinning Division has 
been formed at Diversey which uses 


Behind the tailstock you see a Diversey craftsman —_ the latest and largest equipment to 


ready to take a cut that will trace the bottom flange a 
and taper of the O.D. of the nozzle. As nozzles go, 


this one is big ... real big. It is the nozzle for one of oe die 
the large rocket motors. SEND rai 


At Diversey you have the LARGEST FACILITIES 
exclusively devoted to your missile metal machining 


; problems. You work with fast, precise and progres- FREE 
sive technical people who know what works and what BOOKLET 
won't. Bring your big missile problems to Diversey. 


LEADERS IN CONTOUR MACHINING 


wers ENGINEERING COMPANY 


10550 WEST ANDERSON PLACE 


FRANKLIN PARK, ILLINOIS « A Suburb of Chicago 
FROM NOSE TO NOZZLE, FROM ni TO FIN, amen TURNED PARTS—WITH PRECISION BUILT IN 
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The burning program for a rocket-powered missile mov- 
ing along a rectilinear path is investigated in the light of 
the thermal effects induced on the missile skin. Part 1 
considers the problem of determining the thrust-time 
relationship which minimizes the difference AG between 
the final and initial values of an arbitrarily specified fune- 
tion G(t, m, h, V, T) of time, mass, altitude, velocity and 
skin temperature, for the case where the so-called induced 
drag is negligible with respect to the zero-lift drag. It 
is shown that the totality of extremal ares is composed of 
zero-thrust sub-ares, sub-arcs to be flown with maximum 
engine output and variable-thrust sub-ares. For prob- 
lems not involving time, an explicit solution is obtained 
for the optimizing mass flow as a function of the local co- 
ordinates of the missile. In Part 2 closed form solutions 
are derived under the assumptions of isothermal atmos- 
phere, negligible drag, negligible irradiation and constant 
Stanton number. Particular problems such as mini- 
mum increase in skin temperature are treated within the 
general frame of the present theory. The boundary 
value problem is investigated, and useful criteria are sup- 
plied for constructing extremal paths under various types 
of boundary conditions. Several numerical examples 
illustrate the effect of important design parameters on 
the physical nature of the optimum burning program. 


Introduction 


i ew analysis of the optimum burning program for a rocket- 
powered vehicle has received considerable attention in 
recent years. The problem can be stated as follows: A 
rocket engine, capable of delivering a variable thrust, is ap- 
plied to a vehicle of given configuration moving along a path 
of prescribed geometry. It is required to determine the 
thrust programming technique which minimizes the difference 
AG = G, — G, between the final and initial values of a speci- 
fied function G of the coordinates of the vehicle. 

Concerning vertical flight, previous investigations have 
heen carried out by Hamel (1)? and by Tsien and Evans (2) 
for problems where G = —m. A broader treatment is due 
to Miele (8) for problems where G = G(h,V,m,t). With 
regard to horizontal flight, the papers by Hibbs (4) and by 
Cicala and Miele (5) must be mentioned in connection with 
problems where G = —X, X denoting horizontal distance. 
More general questions, of the form G = G(X,V,m,t), have 
heen considered by the author in (6). 

In all these investigations, an optimum burning program 
was determined regardless of thermal effects induced on the — 
missile skin. Mathematically speaking, the only auxiliary | 
conditions accounted for in the variational problem were the 
equations of motion. 
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It is to be noted, however, that in the upward flight phase 
of sounding rockets or of long-range ballistic rockets and, in 
general, for all vehicles traversing the earth’s atmosphere at 
high speed, important heat transfer phenomena occur from 
and to the surrounding medium. Specifically, the missile 
skin receives heat energy from the boundary layer by con- 
vection, and from solar, terrestrial and interstellar sources 
by radiation; in turn, the skin emits radiant energy according 
to the well-known Stefan-Boltzmann law. 

The increase in temperature of the skin is functionally re- 
lated to the maneuver performed by the vehicle, i.e., func- 
tionally related to the thrust program for the case where 
the geometry of the flight path (or an equivalent condition) 
is prescribed. For missiles designed for very high Mach 
numbers, the temperature rise may be such as to dictate 
structural design, because of the loss of strength of materia!s. 
Thus, it becomes important to determine a thrust program 
which is optimum in the light of heat transfer phenomena. 
Mathematically speaking, the auxiliary conditions to be con- 


sidered for the Mayer problem are the equations of motion 


and the heat balance equation (7). 


Fundamental Hypotheses and Equations of 
Motion 


The following hypotheses are used throughout the paper: 


(a) The rocket-powered vehicle is ideally regarded as a 
particle of mass m variable with the time ¢. 

(b) The thrust 7’, is tangent to the flight path. 

(c) The equivalent exit velocity V, of the rocket engine is 
regarded as a constant, independent of the operational con- 
dition of the engine. 

(d) The engine is capable of delivering all mass flows 
8 bounded between a lower value B = 0 and an upper value 
B= Bmax. 

(e) The aerodynamic lag is disregarded; i.e., the aero- 
dynamic forces are calculated as in unaccelerated flight. 

(f) The acceleration of gravity is a constant. 

(g) The trajectory is straight and inclined at an angle 6 
with respect to a horizontal plane. 

(h) The “induced” drag is negligible with respect to the 
zero-lift drag.* 

In view of the hypotheses (a, b, g and h) the dynamic be- 
havior of the rocket-powered vehicle is represented with the 
following set of differential equations 


(1] 
Js=V + gsin0 + (D — Ts)/m = 0........[3] 


where h is the altitude, V the velocity, g the acceleration of 
gravity and D the drag. The dot denotes derivative with 
respect to time. 

3 Notice that the hypotheses (g and h) limit the applicability 
of the present theory to vertical or near-vertical trajectories, as 
is the case with a sounding rocket and with some types of ballis- 
tic missiles. 
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_ The Thrust Function 
ae Because of hypothesis (c) the thrust 7's is expressed as a 


linear function of the mass flow B 
In turn, hypothesis (d) implies that the inequality 
[5] 


4 be satisfied by the class of all arcs investigated in the present 

analysis. 

__ The limitation [5] which the size of the engine imposes on 
the mass flow 8 has a critical importance in determining the 
behavior of the solution, as has been recognized in (3, 5, 6). 
Several devices have been proposed in the literature in order 
to handle the above inequality. The one considered most 

appropriate is based on the idea of a parametric representation 

of the engine characteristics (3 and 6). Such a technique is 
developed in this article. 

7 The mass flow 8 is represented as a function of a parameter 

q@ having the properties: 

(a) For — © < a< a, the mass flow is B = 0. 

(b) For as < a < + ©, the mass flow is B = Bmax. 

(c) For a, < @ < a, a one-to-one correspondence is 
assumed between 6 and a. 

With the above scheme, a is considered as the independent 
parameter of the rocket engine and is allowed to vary between 
_—oand+o. The mass flow 8 becomes a dependent quan- 
tity, varying between 0 and £,,,x, according to the scheme 
of Fig. 1. In turn, the thrust 7’, is limited between 0 and 
4 Tx max according to Equation [4] and Fig. 1. Notice that 
the condition dB/da = 0 represents either a coasting flight 
— ora flight with maximum engine output. On the other hand, 

_ dB/da # 0 represents any other operating condition of the 

engine intermediate between the two limiting ones. Notice 

also that @ is only a parameter and that there is no necessity 
of attributing to it any special physical meaning. 


Atmospheric Properties 


In Part 1 it is assumed that the absolute temperature 
_ 7 of the atmosphere in which the missile is flying is an arbi- 
_ trarily specified function of the altitude h. The distribution 
of static pressure p and speed of sound a vs. altitude are 
~ consequently determined. The particular sub-case of an iso- 
_ thermal medium is considered in Part 2. 


The Drag Function © 


Because of the hypothesis (h) the total drag of the missile 
_ D is identical with the zero-lift drag. The latter can be 
_ decomposed into pressure drag and friction drag. The pres- 
_ sure drag depends on velocity and altitude, even accounting 
for compressibility and viscosity effects. The friction drag, 
in turn, depends on velocity, altitude, and on the distribution 
of skin temperature of the missile. In the present analysis 
the effects of variations of frictional drag caused by variations 
in the distribution of skin temperature are neglected. As a 
- consequence, the drag function is represented as 


[6] 


Heat Transfer Equation 


_ In accordance with the established practice, the calcula- 
_ tion of the transient temperature at a point of the missile 
skin is carried out by assuming that heat transfer is polarized, 
in the sense that the speed of propagation of thermal effects 
is infinite in the direction normal to the skin and zero in 
any direction tangential to the skin. 
In the light of the above hypotheses, a one-dimensional 


Fig. 1 Parametric representation of the engine characteristics 


analysis becomes possible for each point of the missile skin, 
according to the differential equation 


= Q+ Ge — [7] 


where the first member denotes heat gained by the unit area 
of missile skin per unit time. Terms appearing on the right 
have the following significance: Q is the heat transferred 
from boundary layer to skin per unit time and unit area; 
eG, is the heat received by the skin per unit time and unit 
area because of irradiation associated with solar, terrestrial 
and interstellar sources; ¢o%7 is the radiant heat emitted 
by the missile skin per unit time and unit area.‘ 

In Equation [7], 7 is the skin temperature; C, = C,(T), 
p, and 6, are, respectively, the specific heat (referred to the 
unit mass), density and thickness of the skin; «€ = ¢(7') is 
the emissivity, assumed equal to the absorptivity; o is 
the Stefan-Boltzmann radiation constant; G,» is the incident 
irradiation associated with solar, terrestrial and interstellar 
sources and depends, in general, on the altitude h and on the 
angle of incidence between surface and solar rays. For a 
given missile design and for a given point of the missile skin, 
the parameter Q is a complicated function of free-stream 
Reynolds number, free-stream Mach number and ratio of 
wall temperature to free-stream temperature. Since the 
Reynolds number and the Mach number depend on velocity 
V and altitude h only, one concludes that, for the present 
problem, Equation [7] has the form 


€(oxT 4 Gs) sand Q 
Ds 


[9] 


A given point of the missile skin is considered in this article, 
and Equation [8] is applied to determine the transient tem- 
perature at such a point. In the case in which the loca- 
tion in question is the one where the most unfavorable thermal 
effects are expected, it may be anticipated that the predicted 
skin temperature is higher than the actual one. In fact, the 
one-dimensional hypothesis embodied in Equation [8] rules 
out any alleviation effect which the most critical parts of the 
missile skin may receive from the surrounding elements of 
matter. 


‘ Equation [7] neglects heat transfer phenomena between the 
skin and other components of the missile. 
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PART 1. 


In this section of the paper, the optimum burning program 
is studied in connection with an arbitrary distribution of at- 
mospheric temperature r vs. altitude h, a drag function of the 
form D = D(V,A) and a thermal function of the form ¢ = 
¢| hyv,T). 


Variational Approach 


Generalized Solutions 


The set of differential equations [1—3, 8] is now considered, 
where Tx = BV,, B = B(a) and D = D(V,h). The time ¢ 
is assumed as the independent variable. The dependent 
variables are the mass m, the altitude h, the velocity V, the 
skin temperature 7' and the engine parameter a. The num- 
ber of differential equations is four while the number of un- 
known functions is five. One degree of freedom is left and 
an optimum requirement can be, therefore, imposed. 

After prescribing the initial coordinates t; = 0, m; = m(0), 
h. = h(0), V; = VO), T; = T(0), and some, but not all, of 
the final coordinates t;, m,;, h;, V,;, 7, the variational prob- 
lem, of Mayer type, is stated as follows: Among all sets of 
functions m(t), h(t), V(t), T(t), a(t) satisfying Equations 
||-3, 8] and the prescribed end-conditions, to determine the 
special set which minimizes the difference AG between the 
final and the initial value of an arbitrarily specified function 
G = G(t,m,h,V,T) of the coordinates of the missile. Par- 
ticular cases of this general problem are, for instance, the 
case G = T, minimizing the increase in skin temperature; 
and the case G = —m, minimizing the propellant consump- 
tion. 


Euler Equations 


A set of variable Lagrange multipliers \,(d), A2(t), As(é), 
\,(2) is introduced and the following expression formed 


where J,, Jo, J3, J, denote, respectively, the first members of 
Equations [1-3, 8]. Since the unknown functions are five 
in number, one must write five Euler equations. These are 
written 


d | oF oF 

where 2, = m, 2 = h, 2, = V, 2; = T,2; = a. The explicit 


form of Equation [11] is 


V.B D 
m 
° D 
m 
D, 
A; = —)e sin 6 + As + hg, {14] 
V 
da m 
oD oD 
ah’ 
ov 
Fr av’ Yr [18] 


First Integral 


Since the fundamental function F, given by Equation [10], 
is formally independent of the time, the following first inte- 
gral holds 
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7B —D 
+ AV sin 6+ A; — gsin a| 
[19] 


where C, is an integration constant. 


Discontinuity of the Eulerian Solution 
The Euler Equation [16] is particularly interesting because 
it shows that the extremal arc is discontinuous, being generally 
composed of sub-ares of equation 


da 
and sub-ares of equation 
m 


According to the parametric representation indicated in 
Fig. 1, Equation [20] represents either a coasting flight or a 
flight with maximum engine output. Equation [21], on 
the other hand, represents a flight condition with a continuous 
variable-thrust, as shown in the following sections. 


Erdmann-Weierstrass Corner Conditions 


In view of the discontinuous character of the solution, the 
Erdmann-Weierstrass corner conditions (8) must be applied. 
These are continuity conditions to be satisfied at every corner 
of the discontinuous extremal solution by the six quantities 


F OF OF OF OF S. OF 

Om Oh OV OT O21 
Since F is independent of @, the continuity of OF /0a@ is 
inherently verified at all corner points. The analogous re- 
quirement for the remaining five quantities [22] leads to 


where the subscript (—) denotes a condition immediately 
before the junction and the subscript (+). a condition im- 
mediately after the junction. The Erdmann-Weierstrass 
conditions require, therefore, (a) that the multipliers 4,(4 = 
1,2,3,4) be continuous at all junction points, and (b) that the 
integration constant C,; be the same for all sub-arcs of the 
discontinuous extremal solution; this circumstance implies 
that the quantity A, — (A3V./m) be zero immediately before 
and immediately after each junction point. 


Boundary Conditions® 


For the problem under consideration, the boundary condi- 
tions include a number of fixed end-point conditions plus a 
number of natural conditions. The latter must be obtained 
from the following general transversality condition, which is 
to be identically satisfied (8) for all systems of variations 
(6t, dz) consistent with the prescribed end-conditions 


[6G + Adm + + ASV + AGT — = 0........ [25] 
Closed Form Expressions for the Distribution of 
Lagrange Multipliers 


As the analysis of (9) shows, closed form expressions can be 
derived, yielding the Lagrange multipliers \,(k = 1,. . .4) at 
all points of the variable-thrust sub-are, defined by Equation 


5 The present article is a condensed form of the investigation 
described in (9), where more extensive analytical details are 
available. 
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[21]. In this connection, extensive manipulations yield the 


following results 


D isin 


(V¢e, — ¢)sin0 


( V + gtsin 0 


) | (7 + gf sin *) 
.) | exp 


Ve, - 
. [29] 


where C, is an integration constant. Notice that these equa- 
tions hold regardless of the boundary conditions and what- 
ever the shape of the function G whose difference AG is to 
be minimized. 


Problems Where No Time Condition Is Imposed 


For the particular case where the G-function has the form 
= G(mh,V,T) and the final time instant is free of choice 
(6t; = arbitrary), the transversality condition [25] yields 
C; = 0. 


Problems Not Involving Heat Transfer _— 
14 


If the G-function has the form G = G(t,m,h,V) and the final 
temperature is free of choice (67', = arbitrary), the transvers- 
ality condition [25] yields \y, = 0. Because of this con- 
dition at the final point, the Euler equation [15] leads to 
A; = 0. at all instants of time. As a consequence, Equation 
[29] reduces to 


D 


C V + gtsin 0 
exp | - + = =o. . [30] 
and there fore yields, for @ = 


by the writer in (3). 


= 
m/2, a result already obtained 


Problems Not Involving Heat Transfer Where No Time 
Condition Is Imposed 


For the sub-case where the G-function has the form G = 
G(m,h,V) and both the final skin temperature 7’, and the 
final time instant ¢; are free of choice, the simultaneous con- 
ditions A; = 0 and C,; = O must be applied. As a conse- 
quence, Equation [30] leads to 


D 
D+ mg sin @ — V (2 +.) [31] 
For 6 = 2/2, Equation [31] supplies a result which is im- 
plied in the equations already derived by Tsien and Evans 
(2), for one particular set of boundary conditions. 


1234 


- [15, 29] yield for C,; = 0 


[27] 


Optimum Mass Flow for Problems Where No Time 
Condition Is Imposed 


An explicit expression for the optimum acceleration (or 
for the optimum thrust) can be obtained if the Lagrange 
multipliers are eliminated. In this connection, Equations 


V 
32 
g Vs + 


+ VV.D,, — 


V(D + V.D,) — V.[D + mg sin 0] 


[33] 
VV. — Veen) sin + — PPr [34 
o- Ve. 
D + 2V.D, + V2 Dy. tee 
V(D + V.D,) — VD + mg sin @) 
[36 
y= Ve, 
Dre = = i 
02 
Cry = — Pun = = = [38 
ov? oVoh oVol 
= 


The drag D and its derivatives depend on velocity and alti- 
tude only. In turn, the thermal function ¢ and its deriva- 
tives depend on velocity, altitude and skin temperature. 
As a consequence, it is concluded that the y-function has the 
form yw Y(m,h,V.T). The instantaneous mass flow at 
points of the variable-thrust sub-are of the extremal solution 
is consequently given by 


VB = sin + 


In conclusion, the optimizing mass flow 8 can either be 
B = 0 or B = Bmax or can have the value which satisfies 
Equation [39]. The latter holds whatever be the distribu- 
tion of atmospheric properties vs. altitude, the drag function 
D and the thermal function ¢; it also holds for all G-functions 
of the form G = G(m,h,V,T) provided the final time instant 
is free of choice. 


D 
mg 


PART 2. Closed Form Solutions 

In Part 1, an explicit solution has been obtained for 
the optimum mass flow. For the variable-thrust sub-arc, 
this optimum mass flow is defined by Equations [32-39]. 
Its form, however, is such that analytical solutions of the 
system composed of heat transfer equation and equations of 
motion are out of reach. Numerical analyses are in order and 
digital computers must be used. Prior to any attack of the 
problem with digital computing equipment, however, it 
appears desirable to achieve some qualitative understanding 
of the physical nature of the optimum thrust program, 
even if drastic hypotheses are necessary to attain such an 
objective. 


In the following section, closed form solutions are derived 
for problems where G G(mh,V,T). The following hy- 
potheses are used in combination with those previously listed: 


(a) An ideally isothermal atmosphere is assumed, i.e., an 
atmosphere where the speed of sound is constant. 


Additional Hypotheses 
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(b) The drag of the missile is considered negligible with 
respect to thrust and weight (D = 0). ht ot K; m(L — log m) = Ce 

(c) Irradiation terms are neglected in the heat balance = 
equation. where m = m/m, is a nondimensional mass, M = V/a the 

(d) Use is made of a modified Stanton number St, based Mach number, u = V,/a the ratio of exit velocity to atmos- 


on free stream density p, free stream velocity V and on a pheric speed of sound, and 7? = fg sin 6/a a nondimensional 


conventional adiabatic wall temperature 7’, corresponding time. The terms appearing on the right of Equations 
to free-stream conditions; such a Stanton number‘ is regarded [48-50] are integration constants and 
as a constant; moreover, constant values are also assumed V. Bs 
for the thermal properties of air and skin and for the recovery yp eee. [51] 
gm, sin 6 
factor r. 
In view of hypothesis (a) the distribution of relative For the ideal case of a pulse-burning sub-arc (K; = ©) 
density o vs. altitude h is given by Equations [48-50] imply the constancy of the three quanti- 
o = p/p = exp(—Yh)............-. 140] ties, h, 7, and M + u log m. 
where / is a nondimensional altitude, p is the absolute density For B = 0 the equations of motion yield the obvious results 
o! air at altitude h, and po is the absolute density at sea level. 
Because of the neglect of irradiation terms, the thermal func- m = (%........ 
tion ¢ takes the form 
where Cz, Cs are integration constants. 
where Q is the heat input from boundary layer to skin per : 
unit time and unit wetted area. The latter can be rewritten Viaslsiiie ties Gelicten 
as 
For the variable thrust sub-arc, the optimum acceleration — 
Q reduces to 
43 
C,oV(T, — T) 143] 
II 
C, pp a? 
[56] 
where M is the free-stream Mach number, 7 the free-stream C, p, 96, 3y sin 8 


temperature and y the ratio of specific heat at constant 
pressure C,, to specific heat at constant volume C,. Simple 
manipulations yield the following useful form for the thermal 


At this point a change of independent variable is operated in— 
the sense that the relative density o is now assumed as the 
new independent variable. Equation [55] is rewritten as 


function 
<i dM 1 
T = | Kk, — 57] 
Kel | kar]. [45] M ( x) do 
; admitting, therefore, the general integral 
[46 
Ce Psd, [58] 
W/o 
[47] where Cy is a constant. The differential relationship between 
7 7, time and relative density is stated as 
Integration of the Equations of Motion . . do 
exp (Kyo) 
Sub-Are Flown With Maximum Engine Output 
1 
Because of the hypotheses of the previous paragraph the = [60] 
sub-are B = Bmax is identified by : / 
The integration process for Equation [59] leads to’ 
[48] 
A(o) = 3 Wo (-1)" (62] 


n\(3n + 1) 


0 


® For a given missile design and for a given point of the missile 
skin, the Stanton number St is a complicated function of free- 


where C\, is a constant. The distribution of mass along the _ 


stream Reynolds number Re, free-stream Mach number M and variable-thrust sub-are is calculated from : 

ratio of wall temperature 7’ to free-stream temperature r. In 

spite of this, the following sections of this paper refer to the ideal = : M+i 

case where the Stanton number is a constant, independent of M, m = Cz exp | — es 

Re T/r. It cannot be denied the above is 

quite drastic in its nature. Nevertheless, the writer feels that “— Si . 

the essential facts of the mechanism of heat transfer are retained where C12 is a constant. 

in the present idealized scheme. The latter, on the other hand, ee 

has the merit of leading to closed form solutions and therefore 7 A discussion of the manipulations leading to Equations [61, — 
to clear-cut information on the qualitative nature of the opti- 62] can be found in (9). In many practical cases, the term 
mum burning program. Moreover, as the subsequent analysis Kyo is at most of order one. As a consequence, no appreciable 
shows, the optimum burning program is affected only in a minor error is introduced into the results by neglecting the terms of © 
way by the hypothesis concerning the Stanton number. exponent larger than 3 in the expression for A(c). —_ 
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2 Mach number-altitude relationship along the variable- 
thrust sub-arc 


Discussion of Results 


As the previous analysis shows, the thermal parameter 
Kz, defined by Equation [56], is systematically present in the 
solutions. Such a parameter depends on the thermal prop- 
erties of air and skin and on the Stanton number, which char- 
acterizes the convective process from boundary layer to skin. 
The lowest value for the above thermal parameter is K,; = 0 
It corresponds to the ideal case of a skin of infinite thermal 
capacity. Increasing values of K, are associated with de- 
creasing thickness of skin. 


Mach Number-Altitude Relationship 


The ratio M/ Mp, of Mach number at altitude to Mach num- 


ber at sea-level is plotted in Figs. 2 and 3 as a function of the 
nondimensional altitude 4, and of the thermal parameter 
Ky. For Ky < 3 (skin of relatively high thermal capacity) 
the Mach number increases monotonically as the altitude 
increases. On the other hand, for K, > 3 (skin of relatively 
low thermal capacity) the Mach sianaet initially decreases, 
has a minimum at the altitude where o = 1/(3K;,) and after- 
wards increases as the altitude increases; in other words 
after entering the variable-thrust sub-arec, the velocity of the 
missile must be decreased to limit the amount of heat trans- 
ferred from boundary layer to skin in the low altitude region. 


Acceleration-Altitude Relationship 


The ratio V/V» of optimum acceleration at altitude to 
optimum acceleration at sea level is indicated in Fig. 4 as a 
function of the altitude ), for two values of the thermal 
parameter K, The diagram shows that the acceleration in- 
creases with such rapidity, as ) increases, that there is a 
very definite altitude above which the optimum burning 
program may become impractical from a structural stand- 
point. Consider, for instance, a missile whose thermal con- 
stant is K,; = 0.2 and assume that the initial acceleration 
is Vo/g = 0.5. Assume also that the maximum permissible 
acceleration is (V/g)max = 15. As Fig. 4 shows, the optimum 
burning program can only be used at altitudes below’ 
approximately 100,000 ft. Analogously, if it is assumed that 
the maximum permissible acceleration is (V/g)max = 40, 
one concludes that the optimum burning program can only 
be used at altitudes® below 130,000 ft. 


Integration of the Heat Transfer Equation 


After assuming the relative density o as the independent 


® The figure in question is based on an ideal isotherm: al at- 
mosphere where a = 1000 ft sec~ a 
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Fig. 3 Mach number-altitude relationship along the variable- 
thrust sub-arc 


S|<. 


Fig. 4 Acceleration-altitude relationship along the variable- 
thrust sub-arc 


variable the heat transfer equation [8] can be rewritten as 


» 
— 3K,T + + K.M*) =0........ [64] 


where 7’ = 7'/r is the ratio of skin temperature to free-stream 
temperature and K, is defined by Equation [56]. For a pre- 
scribed Mach number-density relationship M(o), the above 


JET PROPULSION 


iy . 
22 50 
VA 
18 
/) 
M / Z| 
3 
> 
14 
10 
4 5 Wa 
? 
Z| 
| 
; 3 4 5 7 8 p 
E 
> 
10,000 
J 
1000 A 
100 
w 
su 
al 
| 
co 
is 
a 
ta 
ta 
D 


equation can be regarded as linear in the unknown function ” | 
7(c), yielding the following general integral | 
T = 1+ [const — f M? exp (—3K,a)da] exp [65] | | 
heh, | heb, for the 
Sub-Are Flown With Maximum Engine Output ‘ i | 
For this type of sub-are the temperature equation must | “7 
be dealt with, in general, by approximate methods. An J pong eon gli 
important exception is represented by the pulse-burning case i « | | 
in which Equation [65] reduces to | | 
+ T 
Fig. 5 Altitude where the skin temperature is stationary 
Coasting Sub-Are (variable-thrust sub-arc) 


For the coasting sub-arc, Equation [65] yields - 


2 
T =1+4+ + Cy exp (3Kio) + Kz C(a).... {67] 


C(c) = log + [log (3K,) — Ei( —3Kyo)] exp [68] 
where Cy, is a constant and Ei(—3Ky,c) is the so-called ex- 


po:ential-integral function, defined as 


Ei( -3K,o) = f 


@ 


@ K n 
0.5772 + log + [69] 
1 


Variable-Thrust Sub-Are 


For the variable-thrust sub-are the skin temperature- 
density relationship is supplied by 


T = 1+ [Cis — 3K2K,C.2A(o)] exp (3Kyc)..... [70] 


where C}; is an integration constant. 


4 
7 | if 
An Upper Limit for the Increase in Skin Temperature _ 
4 14 
As the subsequent analysis shows, an upper limit can be 
detected for the increase in skin temperature of a missile ¥ 
which flies according to the optimum burning program. As- Fig. 6 Skin temperature at ; = © and at stationary points 
sume, for instance, that at the moment of entering the vari- (variable-thrust sub-arc) 


able-thrust sub-are, the flight condition is represented by 
h =0,0 =1,M = = 1. The integration constants 
Cy and C\; are subsequently calculated and Equation [70] 
rewritten as 


32 T 32 
T 3K,[A(1) — A(o)] exp [K(30 2)].... (71 
— = — A(a)] exp — =I 
To investigate whether Equation [71] has any stationary 24 ia 424 
points, the condition d7/do = 0 is imposed, yielding 
3K,0"/* exp (Kio) [A(1) — A(o)] = 1........ [72] 
416 
Equation [72] is plotted in Fig. 5 from which the following 
conclusions are drawn: KMS 
(a) Fora thermal constant A, < 0.99 the skin temperature 
is a monotonically increasing function of the altitude. As | tos 
a consequence, an upper limit to the skin temperature is | | 
supplied by the ideal condition which the missile would at- a 
tain should the optimum burning program be continued Kemi 
up toh = ©, ie, ¢ = 0 08 06 
ea Ai) Fig. 7 Adiabatic wall temperature and skin temperature at 
K2M¢? points of the variable-thrust sub-arc 
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Fig. 8 Adiabatic wall temperature and skin temperature at 
points of the variable-thrust sub-arc 

(b) For a thermal constant K, > 0.99 there are two alti- 

tudes at which the skin temperature is stationary, the lower 
altitude }, corresponding to a relative maximum for the tem- 
perature 7, the higher altitude he to a relative minimum 


T:. The temperature rise between sea level and the altitude 
defined by Equation [72] is supplied by 


= o~** exp (2Ki(o — 1)].......... [74] 
where the density ¢ is to be consistent with Equation [72]. 
In Fig. 6 the temperatures defined by Equations [73, 74] 
are plotted as a function of the thermal constant Ky. As the 
graph indicates, for 0.99 < K, < 1.4, the stationary condition 
reached at h = hy is not critical and the absolute maximum 
for the temperature is reached ath = ©. To emphasize the 
significance of the above results, Figs. 7 and 8 have been 
prepared. They yield the nondimensional skin temperature 
T = T/r and boundary layer temperature 7, = T',/7 as a 
function of the relative density o at points of the variable- 
thrust sub-are. For a thermal constant K, < 3 (for instance 
K, = 0.2) the optimum Mach number is a monotonically 
increasing function of the altitude (Figs. 2 and 3). As a 
consequence, the boundary layer temperature 7, increases 
with the altitude: Heat energy is continuously transferred 
from boundary layer to missile skin, as shown in Fig. 7 
On the other hand, for a thermal constant K, > 0.99, i.e., 
for a skin of relatively low thermal capacity, the Mach num- 
ber-altitude diagram presents a minimum point. The dia- 
gram of the boundary layer temperature, therefore, has also 
a minimum point. The skin temperature rise becomes such 
that an inversion occurs in the heat transfer process. This 
interesting phenomenon is shown in Fig. 8 which refers to a 
thermal constant K,; = 1.4: Heat energy is released from 
skin to boundary layer in the altitude interval oes 
0.112 < o < 0.657. ; 


Solution of the Boundary Value Problem 


In the previous sections, general solutions have been de- 
rived for the problem of minimizing the difference AG be- 
tween the final and the initial value of an arbitrarily specified 
function G = G(m,h,V,T). It has been shown that the ex- 
tremal arc is composed of sub-arcs B = 0, sub-arecs B = Bmax 
and sub-ares to be flown with regulated thrust. 

Now a particular form of the G-function is considered, 
namely, G = T (minimum increase in skin temperature). In 
this connection, the boundary value problem is investigated: 
It consists of determining the special combination of sub-arcs 
which satisfy a set of prescribed end-conditions, more specifi- 


1238 


maximum thrust coasting sub-arc 


sub-arc 
h h 
h 
B 
variable thrust 
sub-arc 


variable thrust 
sub-arc 


I 
@ 


Fig. 9 Combination of sub-arcs for minimum increase in skin 
temperature 


cally: 2; = 0,m; = 1,h; = 0, M; = 0, 7; = 1; m; = given, 
M, = given, h, = given. 

Two main types of extremal trajectories may exist (9). 
They are indicated in Fig. 9 with the roman numerals I and 
II. The trajectories of type I include: An initial sub-are 
IA flow with maximum thrust, a central sub-are AB flown with 
variable thrust, and a final sub-are BF flown with maximum 
thrust. The trajectories of type II differ from the previous 
ones in so far as the final sub-are BF involves a coasting flight. 

In the following sections the ideal case of an engine ca- 
pable of delivering all mass flows between 0 and © is con- 
sidered, for simplicity. The reader, however, may find an 
extensive treatment in (9), for the case where the maximum 


» The» Function 


The following function of the local coordinates of the missile 
is now defined 


exp (Kyo) 
[A(o 


Ve — A(1)] — yM[M + w log m].. [75] 


Let the symbol w, denote the value of w calculated at the 
final point. As the analysis of (9) points out, the final point 
F must be reached: (a) With a coasting flight if w, < 0; 
(b) with engine operating at maximum output if w, > 0; 
(ec) with regulated thrust ifw, = 0. 


Determination of Corner Points 


A fundamental aspect of the boundary value problem is to 
determine the corner points A and B of the extremal solution. 
In general, an indirect procedure must be used (9). For the 

_ particular case Bmax = ©, however, a direct approach is 


possible. 


Case (l)—The Final Sub-Are BF Is Flown With Maximum 
Engine Output 


As the analysis of (9) points out, the Mach numbers at 
points A and B are supplied by 


ou, = Alo,) — AQ) 
A My + wlog my 
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Fig. 10 Extremal trajectories of type II for different values of 
the thermal parameter K,; 


Case (I)D—The Final Sub-Are BF Is Flown With 
Engine Shut-Off 


As it is shown in (9) the following transcendental equation: 
determines the relative density at point B ; 


-¥ M,? + - log | u log m; + \ M,? + + 
oF f 

(Kyon) 


V 


Once op is known, the Mach numbers at points A and B are 
determined from 


Ma = \ M ,? log =< oy exp [K,(1 — op)]... [79] 


With the object of illustrating the previous approximate 
theory, several numerical examples have been carried out. 
These examples are now described. 


Numerical Examples 


Burning Program for Minimum Temperature 
Increase in a Sounding Rocket for Several Values of K, 


The problem of minimizing the temperature increase in a 
sounding rocket is now considered in connection with the 
following set of conditions 


m, = 1, hi = 0, T,=1, M,=0 


m; = 0.4, hy = 12.88, M, 


= 0.19, u = 7, 


Several values are considered for the thermal parameter 
Ky, namely, Ky = 0.2, Ky = 0.4 and Ky, = 0.6. The first 
step is to determine the sign of the function w, defined by 
Equation [75], at the final point. Due to the fact that 
w, <0, one concludes that the final point is to be reached with 
a coasting flight. The second step is to determine the altitude 
hy, where the transition from regulated thrust to coasting 
occurs. In this connection, the transcendental Equation 
78] must be solved. The Mach numbers M, and My are 
calculated from Equations [79 and 80]. 
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Fig. 11 Extremal trajectories of type II for different values of 
the thermal parameter K, 


The subsequent step is to determine the distribution of 
Mach number, mass, time and temperature with the equa- 
tions of paragraphs under “Integration of the Equations of 
Motion” and “Integration of the Heat Transfer Equation.” 
The Eulerian paths include a pulse-burning sub-are IA a 
central sub-are AB flown with regulated thrust and a final 
coasting sub-are BF, and are described in Fig. 10, from which 
the following conclusion is drawn: The optimum burning 
program is quite insensitive to changes in the thermal properties 
of air and skin. The temperature rise 7 — 1, on the con- 
trary, strongly depends on the value of the thermal parameter 
Ky, (Fig. 11). 

In the light of this statement, the writer feels that the hy- 
pothesis of constancy of the Stanton number is well-justified 
from the point of view of predicting the optimizing mass 
flow 6 as a function of the time ¢. The errors introduced by 
the Stanton number hypothesis in the calculation of skin 
temperatures, however, should be corrected. In this con- 
nection, the following iterative procedure appears to be 
logical: (a) Calculate the optimum burning program in the 
(h,M,m,t) space by assuming a constant value for the pa- 
rameter K,; (b) determine the distribution of temperatures 
T(h) for constant value of Ay; (c) starting from (a) and (b), 
calculate the true distribution of Stanton numbers St(h), 
as supplied by the boundary layer theory; (d) integrate the 
temperature equation by approximate methods. 


Effect of the Propellant Mass on the Temperature at Final 
Point in a Sounding Rocket 


The problem of a sounding rocket is now considered in 
connection with the set of conditions 


= 1, h; =0, T,=1, M, =0.. [82] 
hy; = 12.88, M, =0, 
K,= 020, = 0.19, u=7 


The extremal path includes a pulse-burning sub-are IA, a 
central sub-arc AB flown with regulated thrust and a final 
coasting sub-are BF, and is shown in Fig. 12 for several values 
of the ratio of propellant mass to take-off mass. Inspection 
of Fig. 12 shows that an increase in propellant mass causes 
a shifting of the optimum distribution of speeds toward 
the region of lower velocities. As a consequence, lower final 
temperatures must be expected when m, = 1 — ™, increases. 
This effect is clearly shown in Fig. 13. 
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Fig. 12 Extremal trajectories of type II for different values of 
the propellant mass 


The particular trajectory associated with m,; = 0.484 
(Fig. 12) has the property of minimizing the propellant con- 
sumption for the case where the drag is ideally zero and no 
condition of a thermal nature is imposed. The temperature 
rise associated with such a trajectory, however, is consider- 
able (Fig. 13). In this connection, Fig. 13 points out one 
important concept: A slight increase in propellant mass 
over the minimum value necessary to reach the prescribed 
altitude h, determines a sharp decrease in the skin temperature 
at final point. The use of extra amounts of propellant, 
therefore, may be an effective means for preventing the oc- 
currence of undesirable effects due to aerodynamic heating. 


6 


Conclusions and Recommendations for Further 


Research 
1 


_ The burning program minimizing an arbitrarily specified 
function of the final values of time, mass, altitude, velocity 
and skin temperature is investigated. For problems where no 
time condition is imposed, explicit solutions are obtained 
for the optimizing mass flow as a function of the local coor- 
dinates of the missile. Closed form solutions are derived under 
the assumption of isothermal atmosphere, negligible drag, 
negligible irradiation and constant Stanton number. Nu- 
merical examples show that the Stanton number hypothesis 
originates only minor errors in the prediction of the optimum 
burning program. 

Due to the approximate character of the theory developed 
in Part 2, further research is necessary to reach a full under- 
standing of the problem of the optimum burning program. 
Because of the practical impossibility of obtaining closed form 
solutions for the general case, two lines of approach should 
be explored. In the first place, digital computing equipment 
should be used to integrate the equations of the optimum 
path by accounting for the effects of the aerodynamic drag 
(neglected in Part 2) for irradiation phenomena and by lifting 
the hypotheses of constancy of the Stanton number and the 
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Fig. 13 Extremal trajectories of type II for different values of 
the propellant mass 


assumption of isothermal atmosphere. In the second place, 
an approach should be undertaken with the so-called direct 
methods of the calculus of variations, such as, for instance, 
the Rayleigh-Ritz method. In this connection, the writer 
feels that much can be done by exploring the class of arcs 
composed of sub-ares 6 = 0, sub-arcs 8B = Bmax and sub- 
arcs where the acceleration is assumed to be proportional to 
some power of the Mach number. - 
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An investigation was made of the over-all characteristics 
associated with unbounded mixing of axially symmetric 
compressible coaxial air streams. Experiments were per- 
formed to determine the effects of what are considered the 
fundamental parameters of the problem (the initial Mach 
numbers and initial stagnation temperature ratio of the 
streams). Included in the range of experiments were pri- 
mary Mach numbers of 2.60, 0.97 and 0.68, secondary Mach 
numbers of 0.45, 0.25 and 0.00, and stagnation temperature 
ratios of 1.00 and 1.429. Mass flows in the mixing region, 
and hence mixing rates were computed from the experi- 
mental data for the range of parameters tested. 

The experimental (and to a lesser extent the analytical) 
results exhibited a very marked effect of the secondary 
Mach number. This effect was greatest for the subsonic 
primary Mach numbers. The over-all spreading was less 
for M, = 2.60 and the effect of the secondary Mach num- 
ber was less pronounced. The effect of increased temper- 
ature ratio was to increase the over-all spreading. 

A practical correlation between the spreading rates was 
found for the subsonic Mach numbers. For a constant 
primary Mach number the spreading angle of the mixing 
region seems to depend linearly on the secondary Mach 


number. 


Nomenclature 
AF 
axial distance from nozzle 
radial distance from axis 


Tr = 
u = velocity in axial direction 
v = velocity in radial direction 
= Mach number 
T = temperature 
= density 
b(x) = mixing width = (r; — ri) 
ro = primary nozzle radius 
q = spre: iding angle 
Tr = shear stress 
m = mass flow 
Subscripts 
1 = primary stream 
a = secondary stream 
s = stagnation 
= edge of core region 
] = outside edge of mixing region 7 
Superscripts 
— = time mean 
‘ = difference from mean value 
7; arbitrary integration limit 
Introduction 


HE mixing of two unbounded coaxial compressible 
streams, which are turbulent, is of considerable interest in 
such problems as the spreading of turbojet or rocket engine 
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The Mixing of Unbounded Coaxial Compressible Streams: 


exhausts in flight, and it is the simplest logical extension of 
free jet mixing, a case which has been examined quite exten- 
sively (1-7).4 

The present investigation was not concerned with the de- 
tails of the turbulence, but rather with the over-all effects on 
the mixing of varying the initial Mach numbers and stagnation 
temperature ratio. An understanding of this mixing process 
should assist in the analysis of the more complex problem of 
mixing in a duct. 

Work has been performed on the problem of mixing in a 
duct (8-11). There are several significant differences between 
these investigations and the unbounded mixing studied here. 
In the duct there can be no radial inflow of mass and the 
axial static pressure can vary appreciably. The bounded data 
which have been presented were either for low Mach numbers 
(0.1-0.2) or were difficult to interpret close to the nozzle exit, 
due to shocks or excessive boundary layers. The various 
relative dimensions of the nozzles and ducts used by different 
investigators made correlation of their data difficult. The 
nozzlesused atthis laboratory have unusually large dimensions 
which minimize the boundary layer effects. They were care- 
fully designed and gave to a very good degree shock-free 
parallel flow. Schlieren techniques were employed to check 
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this condition (7). In the idealized problem, used for the 
analytical developments, the secondary stream is infinite and 
the velocity has a “‘step’’ profile at the nozzle exit as shown in 
Fig. 1. The secondary stream is in practice finite, and is 
eventually destroyed by outside mixing. However, all the 
data presented are for points upstream of this interference, 
and it is therefore believed that they will approximate closely 
those with an infinite secondary stream. The primary nozzle 
radius is now the only physical dimension of the appa- 
ratus involved and hence the data presented should have wide 
application. 

Analytical approaches to the mixing problem started with 
Tollmein in 1926 (1). Two basic approaches exist; the first 
attempts to solve the differential equations of motion, the 
second uses a Kérmdn-type momentum integral with an ex- 
plicit assumption for the velocity profile. Both methods as- 
sume similarity of velocity profiles. Different assumptions 
regarding the shear stress have been used. These include 
Prandtl’s mixing length model, Taylor’s vorticity transfer 
model and the assumption of a constant exchange coefficient 
across the mixing region. No single method appears to corre- 
late all the practical data well, although one method may work 
well in one practical case, and another method in another 
case. A detailed discussion has been given in (12), which also 
contains a selective bibliography. 

In this paper a turbulent Prandtl number of unity is as- 
sumed throughout. The energy equation is eliminated by 
means of the Crocco integral which relates density to velocity. 
The momentum equations then reduce to a single integral 
equation in u, p(u), z, rand u’v’. An attempt to find con- 
sistent similar velocity profiles and shear stress expressions led 
to an equation having no physically correct solutions. The 
correct qualitative trends were given by a simple Kaérm4n-type 
approach using a linear velocity profile and the assumption 
that at uw = (u, + wa)/2, u’v’ = k(u, — ua)? where k is con- 
stant. 

The main emphasis of this paper is to present what are be- 
lieved to be new data on unbounded compressible axisym- 
metric mixing, where clean exit conditions allow significant 
measurements to be taken in the core region. 


Experimental Apparatus and Procedures 
Air Supply System 


A schematic of the experimental apparatus is shown in Fig. 
2. The primary air is stored in high pressure bottles at 3000 
psig. It passes through regulating valves, heating ducts, and 
the primary stagnation chamber and exhausts through the 
primary nozzle. The stagnation temperature is altered by 
allowing a controlled amount of the air to bypass the heating 
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Fig. 3 Photograph of pressure-temperature probe showing 
micrometer adjustment mechanism 


ducts. The stagnation chamber is insulated with internal and 
external ceramic liners to reduce heat transfer between the 
streams, and is provided with screens to insure uniform low 
intensity turbulence at the nozzle exit. Two nozzles, one for 
M, = 2.60 and one for subsonic Mach numbers, were de- 
signed, the former by the Foelsch method (13) corrected for 
boundary layers (14), to give parallel shock free flow at the 
exit. Both nozzles have exit diameters of 2.55 in. 

The secondary air is supplied from two centrifugal blowers 
with rated mass flows at 2800 ¢.f.m. each at standard condi- 
tions. The maximum flow from the blowers for the 107.5 sq 
in. annulus corresponds to M, = 0.50. Normal techniques 
such as turning vanes, egg-crates and screens are employed to 
assure uniform secondary flow. The converging section was 
designed to give a linear change of area with axial distance 
and had a radius of 6 in. at its exit. The temperature of the 
secondary air can be varied in the range 80 to 160 F but the 
control response is sluggish, and the temperature cannot be 
varied appreciably during a test run. oe 


Pressure and Temperature Measuring System 


Profiles in the mixing region were taken using a combina- 
tion total-profile, total-temperature probe, positioned directly 
by a micrometer thread (Fig. 3). This arrangement allowed 
swift and accurate positioning of the probe, which was critical 
with mixing regions of the order of 0.2 in. The impact tubes 
were of normal design. The total temperature tubes, consist- 
ing of copper-constantan thermocouples enclosed by stainless 
steel radiation shields, gave high recovery factors, better than 
95 per cent. 

Both stagnation chambers contained static pressure taps 
and thermocouples. (Flow velocities in the stagnation cham- 
bers were sufficiently low for these readings to approximate 
closely the values corresponding to the stagnation conditions.) 

Pressures were measured with mercury manometers or 
Bourdon gages, and stagnation temperatures with recording 
potentiometers. 

Unlike free jet (M. = 0.00) experiments where the tem- 
perature 7’,, is fixed, variations in 7, may occur during a test 


_ run with secondary air flow. It is therefore desirable to have 
a direct record of the ratio (7,/T..) available during an ex- 


perimental run. To this end a servomechanism was designed 
which utilized the recording potentiometers’ property of 


transforming a d-c signal linearly into a mechanical displace- 


ment, and a direct reading of 7,,/7,2 was obtained. 


Range of Experimental Data 


Surveys were taken across the jet at 0, 1, 2, 3 and 4 primary 
nozzle diameters downstream. Total pressure and, where 
applicable, total temperature measurements were taken for 
various combinations of M, = 2.60, 0.97 and 0.68, M, = 0.45, 
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0.25 and 0.00, and 7,,,/T,. = 1.00 and 1.429. A value of ; 
T.,/T.2 = 1.429 corresponds to over 200 F difference in tem- 
perature between the streams. The pressure readings were 1 
recorded photographically and mixing boundaries estimated 
from the profiles. Check runs were carried out to test the re- 
peatability and over-all accuracy of the data. 

In the case of zero secondary flow comparative tests were 
made with the outside of the nozzle lagged with insulating 
material. With M, # 0.00 it was of course not feasible to in- 
sulate the nozzle without disturbing the secondary flow. 

The effect of disturbing the boundary layers on the inside 
and outside of the primary nozzle was also investigated. 
Simple spoilers made of cotton thread were glued to the nozzle 
1.5 in. from the exit. The thickness of the spoilers was 0.02 in. 
compared with an estimated 0.05 in. for the boundary layer. 
Significantly thicker spoilers were found to affect the uniform 
primary stream at exit, presumably due to separation of the 
boundary layer. Profile checks were carried out at 3 diameters 
and compared with the previous data. These results are re- 
ported in the following section. 

When tests were run with M, = 0.00, the secondary con- - = 
verging section was removed so that data could be com- alias 
pared with the previous work using this apparatus (6). There Fig. 4 Total pressure profile at nozzle exit—/, = 0.68, M, = 
appeared to be a small but significant effect of the configura- 0.25, 7'n/T'. = 1.00 
tion just outside the primary nozzle, presumably due to dif- 
ferences in the induced flows. A systematic study of these ‘e 7 
effects was outside the scope of the present investigation but . 
might be of interest in the future. 7 5 
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Accuracy of the Experimental Data 


The pressures and temperatures were recorded to a high 
degree of accuracy by the manometers and thermocouples, 
respectively. The position of the probe was known to 0.003 
in. and normally better. In reading the photographs of the 
manometers errors of +,5 in. Hg could occur. The chief 
source of error arose when estimating the mixing boundaries 
from the pressure and temperature profiles. Close to a 
boundary a small error in temperature or pressure reading can 
alter the estimate of the width of the region appreciably. 

For the isoenergetic case up to four tests were performed 
with the same initial parameters. Two values of the mixing 
width are obtained from each test as measurements were taken 
above and below the axis of symmetry. (Traverses were R 
made in a vertical plane.) The values reported are mean 
quantities from these repeated tests. An idea of the over-all Fig. 5 Over-all total pressure profile—x = 4 radii, \/, = 0.68, | 
accuracy of the data may be obtained by studying the varia- M, = 0.45 
tions in the values of the spreading widths. The expected 
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errors as a percentage of the mean spreading width were at A te 
worst 8 per cent for a single estimate, and 4 percent forthe | 
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mean estimate with M, = 0.00, and 5.0 and 1.6 per cent, re- 
spectively, for M, # 0.00. Most of the errors were consider- 
ably less than these extreme values. 
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Results® 
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Isoenergetic Mixing Data 


The subsonic isoenergetic mixing cases were used to test the 
repeatability of the data. The values presented are the 
means of several tests and form a basis for comparison with 
the other data. 

Figs. 4 to 8 show typical total pressure profiles for various 
axial distances and initial Mach number combinations. They 
represent only a few of the many similar curves obtained in 
over 200 test runs. Fig. 4 demonstrates the high degree of 
uniformity achieved in both primary and secondary streams 
at the nozzle exit. Also clearly indicated in this figure are the 
small relative dimensions of the boundary layers compared 
with the uniform regions. — 2 


Figs. 5 and 7 depict over-all total pressure profiles. In each 
Fig. 6 Detailed total pressure mixing profiles—x = 4 radii, 


inches mercury 


Total pressure 


- More detailed data are reported in (15). M, = 0.68, Mg = 0.45, Ts:/Ta = 1.00 


| 


can be seen clearly the uniform primary and secondary regions, uniform regions required a certain amount of individual dis- 


the mixing zone, and the effect of outside mixing destroying cretion. Radial distances were measured with respect to the 
the uniform secondary region. At z = 8 (primary nozzle) test table. The position of the axis of symmetry was deter- 
radii the uniform secondary region has shrunk considerably. mined by taking the average value of the distances to the four 
For M, = 0.97 and M, = 0.25 it was in fact so small that no mixing boundaries, two above and two below the axis. The 
attempt was made to estimate the mixing boundaries at this values of r; and r; were then obtained except in the case M, = 
point. No data were taken for x greater than 8 radii. 2.60 when r; was not estimated as the weak shock structure 
The mixing regions on both sides of the jet axis are shown in made this impractical. 

more detail in Figs. 6 and 8. The degree of symmetry of the The major experimental results are reported in Figs. 9 and 
mixing is readily demonstrated with this method of presenta- 10. Fig. 9 is a typical example of the several cross plots that 
tion. Fairing in the profiles up to the limiting values of the can be drawn (15). The dependence of the mixing boundaries 
on the initial Mach numbers is clearly demonstrated. For a 
7 0 constant M, increasing M, decreases the spreading angles and 
as the over-all spreading decreases with increasing M,. There is 
a marked effect of the secondary Mach number. For M, = 
150 0.00, the spreading decreases with increasing M,. For M, = 
4 0.25, the spreading is a very weak function of M,, increasing 
a : slightly with increasing M,. For M, = 0.45, the spreading 
definitely increases with M,. On the hypothesis that the 
-? . mixing is primarily determined by momentum transfer this 

; _ phenomena can be rationalized as follows: 
i A larger velocity difference between thestreamscauses higher 


shear stress and a greater momentum difference. The former 
effect tends to increase the spreading, and the latter to have 
the reverse effect. At M. = 0 it appears that the first effect 
predominates but as M, increases the second effect becomes 

more important and is eventually the determining factor. A 
simple approach, described in the Appendix, 
also gives this same inversion effect. The relative positions 
of the outer mixing boundaries (for simplicity the inner 
boundaries are omitted) are shown in Fig. 10 for the com- 
plete spectrum of initial Mach numbers. 

The data are correlated by assuming the spreading linear 
with axial distance. The lines were chosen to give best agree- 
ment with the mass flows in the mixing regions as described 
below. 


ge 


Mass Flow in Mixing Region 


It is of considerable interest to determine the mass flow in 
the mixing region. As there is no experimental measurement 
'@ 2D of the density in the mixing region, it is necessary to make some 

--1"--+ assumption regarding the variation of density with velocity. 


Total pressure Ga 


Fig. 7 Over-all total pressure profile, 7 xz = 6 radii, M, = 


2.60, Ma = 0.25, 7's:/Tsa 1.00 
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Fig. 8 Detailed total pressure profile—z = 6 radii, 14, = 2.60, Fig. 9 Momentum mixing boundaries—M, = 0.68, 7's,/7's, = 
Ma 0.25, T/T = 1.00 ad 1.00 
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Run 1 

M, 0.68 

0.45 

A B 
xz = 2 radii 0.270 0.307 
x = 4radii 0.478 0.496 
z = 6 radii 0.680 0.708 
x = 8 radii 0.950 0.958 


The Crocco integral is again used. The velocity ratio u/w isa 
known function of total pressure (assuming constant at- 
mospheric pressure). Hence density, and the integral m = 
JS’ pur dr may be evaluated. This method, designated 
method A, involves lengthy and tedious calculation. If it 
could be assumed that both primary and secondary streams 
flow axially up to the mixing region, the integral could be 
evaluated simply as 


m= pur dr = pyuy(ro? — ri?) + patta(rz? — 


This method, designated method B, would require only trivial 


computation as 7r;, r; have already been determined. The 
values of m/3piuro? calculated by these two methods were 
compared for three experimental runs, shown in Table 1. 
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Fig. 10 Outer momentum mixing boundaries of primary jet 


Table 1 Comparison of methods for determining mass flow in mixing region 


Table Linear spreading rates 

——————— Momentum boundariee———————. 

rj — To b(x 

68 0.00 1.000 0.1150 0.1010 
68 0.25 1.000 0.0604 0.0694 0.1298 
68 0.45 1.000 0.0295 0.0446 — 0.0741 
97 0.00 1.000 * 0.1040 | 0.0928 1968 
97 0.25 1.000 a 0.0682 0.0739 0.1421 
97 0.45 1.000 0.0434 0.0552 0.0986 
60 0.00 1.000 0.0710 
60 0.25 1.000 0.0672 
60 0.45 1.000 0.0509 
68 0.25 1.429 ra 0.0639 0.0764 0.1403 
68 0.45 1.429 0.0299 0.0537 0.0836 
97 0.25 1.429 0.0807 0.0848 0.1655 
97 0.45 1.429 0.0552 0.0714 0.1266 


Run 2 


0.97 
0.25 


A 
0.701 

0.956 


For M. = 0.25 and 0.45 methods A and B do not disagree 
by more than 10 per cent. As there is no evidence that the — 
Crocco integral is valid to greater accuracy than this, there is — 
little point in using the much longer and more tedious — 
method A. For M, = 0.00, however, considerable disagree-_ 
ment occurs. This is not surprising, as the assumption in| 
method B is obviously bad for this case when there is in fact a 
small but appreciable inflow of mass radially. 

The mean mass flows, where method A is used for M, =— 
0.00, and method B for M, = 0.45 and M, = 0.25, are de- — 
picted in Fig. 11. Here also a strong dependence on M, can > 
be seen. 

The linear boundaries of the mixing regions were now cal-— 
culated by making the sums — 1?) and Y(ro? — 7,7), 
taken over x = 2, 4, 6 and 8 radii, equal for the practical — 
points, and for the points on the straight lines r; = ro + ja | 
and r; = ro — qi, respectively. For comparison purposes | 
values of qi, gj and q; + q; appear in Table 2. : 

The spreading angle of the mixing region b(x)/x, or q; + qj, 
appears to be a linear function of M, for a given value of M,_ 
(Fig. 12); b(z) is assumed zero for M,; = M, to give two more © 
points on the lines. If the variation of spreading angle with — 
M, is known for the free jet case it should then be possible, — 
for subsonic Mach numbers at least, to predict the spreading _ 
angle for M, ¥ 0 according to the formula 


If the extrapolation of the experimental points to M; = M,. 

is truly linear as would appear from the experimental results, | 
then a very interesting conclusion would be drawn, i.e., two — 
high velocity concentric air streams whose Mach numbers are | 
close but different mix very slowly. Another interesting ob-— 
servation which can be made from Fig. 12 is that the crossover 
of the constant M;, lines gives some idea where the changeover — 
of the importance of momentum effects to shear effects takes — 
place for jets in the Mach number range covered. The posi-— 
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Fig. 11 Mass flow in mixing region 


tion of the M, = 0 points for other Mach numbers would 
follow the trend shown, the lower the M, the larger b(x)/z. 
Warren (6) has shown that the data given in this report so 
correlates with that of other investigators. Since the b(x)/x 
given here is different from the spreading angles correlated by 
Warren, the raw data of the other investigators would be 
required to plot their M,, M, = 0 results in Fig. 12. Un- 
fortunately, these data were not available. 


Heated Stream Data T5,/Tsqa = 1.429 


These data are presented in similar fashion to the isoener- 
getic data in Figs. 13, 14 and 15. It is necessary here to esti- 
mate both the momentum and energy boundaries. It ap- 
pears from recent data for M,; = 2.60, M. = 0.00 and 7’;,/T,.. 
= 1.00 that turbulence and conduction due to the large dif- 
ference of static temperature may be responsible for this lack 
of definition. For the subsonic Mach numbers and 7’,,/7,. = 
1.429 the momentum boundaries spread faster than with 7',,/- 
T.< = 1.00, the percentage change in b(x)/x being larger for 
M, = 0.97 than M, = 0.68. The energy boundaries spread 
faster than the momentum boundaries which is in agreement 
with previous observation on lower speed jets (4). If the 
primary stream has a very high Mach number or is very cold 
so that JT. > 7, it appears that this may not be the case 
and the energy core may be longer than the momentum core, 
due to turbulence and the conduction of heat toward the axis 
being appreciable. Indeed stagnation temperatures higher 
than 7’, have been recorded beyond the core for T,, = Ta, 
M, = 2.60, M, = 0.00, and the term isoenergetic appears 
somewhat of a misnomer in this case. 


Effect of Disturbing Boundary Layer and Lagging Nozzle 


Though small, velocity and temperature boundary layers 
do exist at the nozzle exit. If the mixing were largely deter- 
mined by either of these, then the data presented would de- 
pend on the physical dimensions of the apparatus used and 
would have far less general application. Profile checks carried 
out as described in the previous section revealed that within 
the limits of experimental accuracy no change in the over-all 
mixing could be detected with either boundary layer change. 
Small changes may occur at 1 or even 2 radii downstream, but 
were not detected in the over-all mixing further downstream. 
It is therefore thought that the data presented are to a very 
good measure independent of the apparatus used. 


Conclusions 


The dependence of momentum boundaries on initial pa- 
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Fig. 12 Spreading angle vs. secondary Mach number 


rameters can be satisfactorily attributed to two opposing 
effects of increased velocity difference between the streams; 
namely, increase in momentum difference and increase in 
shear stress. For secondary Mach numbers close to zero 
(<0.25) the latter effect is stronger, but the former dominates 
at higher secondary Mach numbers. The isoenergetic sub- 
sonic data shows a linear dependence of spreading angle on 
M, — M. foragiven M,. Heating the primary jet confirmed 
the results of previous lower speed work in that the momen- 
tum boundaries spread faster. 

The relative positions of the energy and momentum 
boundaries appear to depend on both total and static tempera- 
ture ratios (and/or differences) between the streams and upon 
turbulence effects. For the subsonic heated cases where 
T., > T.. and T, > T. the energy boundaries spread faster 
than the corresponding momentum boundaries. However, 
for M, = 2.60 where T,, > T,. but 7; < T. the reverse holds 
at the core boundary, and the energy core extends beyond the 
momentum core. (This could be significant in, for example, 
positioning JATO units.) 

Analysis of the problem using a Ka4rmén-type approach 
was only qualitatively successful as indicated by the predic- 
tion of an inversion effect discussed. For any good quantita- 
tive agreement, a more realistic treatment of the nonsimilar 
velocity profiles, and a better understanding of the turbulent 
shear will probably be needed. 

The experimental conditions were such that they closely 
approximated those of idealized mixing with an infinite 
stream, and thus the results presented are thought to be es- 
sentially a new contribution to the turbulent mixing field. 
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APPENDIX 


Summary of Analytic Approaches 


Despite considerable researches (1-11), the field of turbu- 
lent mixing is still very much open, particularly in the com- 
pressible case. It seems unlikely that any one of the several 
theories expounded will be suitable for all cases of turbulent 
mixing. However, in specific problems, one or more of the 
theories may give a good correlation with the practical data. 

No attempt was made to evolve a new theory, but to adapt 
ali existing theory to the problem of axisymmetric compressi- 
ble dual stream mixing. Any method that gives a good corre- 
lation with existing data, allows easy inter- and extrapolation, 
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Fig. 14 Momentum mixing boundaries, 7;,/7;, = 1.429 
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Fig. 15 Energy mixing boundaries, 7;,/7;, = 1.429 


uses reasonable assumptions and does not involve excessive | 
computation will be considered satisfactory. 

The following assumptions are used throughout: 

1) normal assumptions for jets, small angles, p constant — 

2) the density variation can be represented by a Crocco in- — 
tegral of the type 7, = T + u2/2c, = A’ + B’u giving p/ 
pi = 1/(D + Eu + Fu?) (turbulent Pr = 1) 

3) the velocity profiles in the mixing region are similar for — 


given M,, M,, T,,/7'o settings, i.e. 


7 


Ul — Ua b(x) 

4) u’v’ may be expressed in terms of the mean quantities in _ 
the form h(é). It will be seen that the following expressions 
fall into this category 


: 
‘ 
q 
+M +0.68 
! 
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an With Assumptions 1 and 2 the equation of motion becomes 


Pe) r® — 
2f dr — wif pur dr = —(r*pu'v’)*...[1] 
z Jo Jo 


4 _§ Using Assumptions 3 and 4, an equation in £ and z results. 


The first method tried to find what A(é), f(£), b(x) and r,(z) 
will satisfy this. By assuming power law relations b(x) = 
gr" and r; — ro = ax" it is found that m = n = 1 are the only 
consistent values. By equating coefficients of x° and x! two 
equations in f(£) and h(¢) are obtained; h(£) is easily elimi- 
nated giving one equation in f(t). Assuming a McLaurin 
series for b(x) and r,(x) — ro would also lead to the same two 
equations by equating coefficients of 2° and z'. Very similar 
equations may also be obtained using u’V’ = h(£)/b(x) and 
power law relations for b(x) and ro — r,(z), when m = n = 


4 are found to be the only consistent solutions. The two 
equations may be written 
coeff 7° = P,(f(é), [2] 
coeff x' A(E) = Pi(fl&), &).. 


Elimination of h(£) gives an equation in f(£) and & It can 
be proved that (regrettably) this equation has no solution 
satisfying the physical boundary conditions on f. It there- 
fore appears that one at least of the assumptions used is in- 
correct. 

For x small it may well be argued that Equation [2] is of 
prime importance. The equation is indeed the one that 
would arise in the corresponding two-dimensional problem 
under similar assumptions. To solve this equation it is neces- 
sary to make a further assumption upon the explicit form of 
h(£). None of the more common expressions for shear stress 
nip to an equation which could be immediately solved. With 

= p, = const the problem reduces to Tollmein’s equation 
with tT ~ kb?(x) (du, Or)’. This solution could be extended 
numerically to cover a given case when p/p; = 1/(a + bf + 
cf?). The method did not give better correlations than the 
simplified model used below and, as it was tedious, effort in 
this direction was discontinued. Details of the various meth- 
ods will be found in Appendix B of (15). 

The second method is essentially an engineering approach. 
An explicit form of the shear stress at the point where u — u. 
= (4)(u; — and a velocity profile are assumed. Inte- 
gration of the fundamental equations with r* = © and r* = 
r; + 6/2 gives two equations which can be solved for r; and b. 
As such methods have been found to be insensitive to as- 
sumed velocity profile (7), a linear profile is selected on 
grounds of simplicity. For the shear stress at the half velocity 


point the assumption used was u’V’ = k(u — ua)?. This 
result would be derived from several functional forms in the 
velocity for example 


~ (Ref. 8) 


= kiumax — Umin 


= kone) (2 


| (= 
u'V' = kiu — umin uz) ( *) 


It is known from practice that at the half = point the 
velocity profile is nearly linear. It would therefore appear 
reasonable that at this point the shear stress is a function 
primarily of the velocity gradient, and any terms involving 

24/dr?, etc., would be small at this point. 

On solving the equations (see Appendix B of (15) for de- 
tails), it is found that to a close approximation r; — ro and b are 
both proportional to xz. These results cannot be foreseen and 
were only discovered on substituting numerical values in the 
equations. 

The values of 6(x)/kx found by this method are given in 
Table 3. 


Table 3 Values of b(x)/kx 


M./Ma 0.45 0.25 0.00 
2.60 8.69 = 11.63 14.20 | 
0.97 5.40 9.00 15.88 
0.68 3.10 7.13 15.92 


It will be seen that for a given M, or M, the correct de- 
pendence on M, or M, is predicted. The actual numbers, 
however, are not simply connected with the experimental 
values for spreading rates. The use of an exponential ve- 
locity profile, which approximated the actual profiles more 
closely, also gave the same trends. The effect of M, for M. 
= (0.00 was more pronounced however. 

It is considered that probably the poorest of the several 
assumptions used in these approaches is that of similarity of 
the velocity profiles. The Crocco integral reduces to the be- 
lievable hypothesis of constant stagnation temperature for 
T., = T.«. The assumed expression for shear stress at the 
half velocity point is consistent with forms that have proved 
satisfactory in similar jet mixing problems with similarity 
well established. 
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Solution of the Transient Heat-Conduction Equation for 


Insulated, Infinite Metal Slab 


JOHN H. GROVER? and WILLIAM H. 


Atlantic Research Corp., Alexandria, Va. 

The transient equation for the transfer of heat from hot heat can be assumed and nearly all of the heat transferred 
combustion gases to an insulated, infinite metal slab has to the metal is conducted through the insulation. 
been solved. A graphical representation of the generalized A solution has been developed with appropriate boundary 
solution has been prepared, as well as a graphical represen- conditions for an infinite slab of insulation with finite thermal 
tation of a simplified, approximate solution of sufficient conductivity, in contact at one surface with a slab of metal of 
aceuracy for engineering use. From these curves, the in- infinite thermal conductivity and in contact at the other 
sulation required to protect a metal slab from hot gases surface with hot flowing gases. Additional assumptions 
can be determined. An illustrative example is given. include: (a) Conduction of heat only in the direction normal 


to the plane of the slabs, (b) an initially uniform tempera- 
ture throughout the slabs, (c) no heat transferred to the 
metal slab except by conduction through the insulation, (d) 
no thermal resistance between the two slabs, (e) a finite 


Nomenclature 


a; = positive roots of tan a;l = (bg — a;?)/(b + g)a; 

b = h/k thermal film resistance between the hot flowing gases and the 
(| = heat capacity of insulating slab insulation, (f) temperature of the flowing gases does not de- 
(, = heat capacity of metal slab rease in the direction of flow and (g) thermal diffusivity, 


g pC/PmCmlm h/pC, of the insulation is invariant with temperature. 
h = coefficient of heat transfer between heating or cooling The basic equation of heat conduction 
fluid (hot gases) and surface of insulating slab 


/ = thermal conductivity of insulating slab or ar 
(= thickness of insulating slab 
/,, = thickness of metal slab Ox? 06 
m = k/hl = dimensionless resistance ratio — 
n = pmCnlm/pCl dimensionless heat-capacity ratio has been solved for numerous conditions (1); however, no 
7 = temperature at position x at time 0 analytical solution has been found for the boundary con- 
7, = temperature of metal slab = temperature of insulation at ditions inherent in the above assumptions.® Jaeger (4) has 
position / at time @ solved the equation for heat transfer through a hollow cyl- 
/’’ = temperature of heating or cooling fluid (hot gases) inder with very general boundary conditions. Although his 
/) = initial temperature of insulating slab solution may be modified to yield the conditions of this 
X = a6/l = dimensionless time parameter problem, the numerical solution of the equation is cumber- 
Y ping some and time consuming. It is also possible to solve the 
eels sillins equation graphically by the method of Schmidt (6). How- 
«a = k/pC = thermal diffusivity of slab (insulation) ever, — method is also very laborious and a simpler method 
= a,l =” was sought. 
6 = time from start of heating or cooling 
» = density of insulating slab te = ‘ . ° 
i = density of metal slab Solution of Equation 
A solution of the basic equation of heat conduction was 
Introduction developed, through the use of Laplace transforms with the 


boundary conditions 

ONSIDER a metal slab separated from a transient flow of 

hot gases by a layer of insulation. This paper pre- T= Toat@ =0............... 
sents a new method for predicting the temperature-time 
history of the metal slab as a function of the pertinent var- 
iables. The results are presented as a family of curves which _ a 
interrelate certain dimensionless parameters. These curves —k — 
are similar to those developed by various authors (3, 7)4 Ox 
for the heating of slabs, cylinders, spheres and other shapes. ar ar, 

The curves were developed primarily for use in calculating —k = PnCmlm =1....... 
the amount of insulation necessary to protect rocket motor Ox 
walls from being overheated by hot combustion gases. How- 
ever, they may also be used to determine the insulation 
needed to protect thin-skinned missiles from aerodynamic tT’ —-T 
heating. In addition, the curves could be used to solve T’ —T, 
many transient heating problems where undirectional flow of : 
+ g?)(b sin ax + a; cos air) 


Received Dec. 11, 1956. at 2 2 2 
1 The figures and equations presented in this paper were pre- ((bg a2)? + + + (bg + + 


pared under the sponsorship of the Navy Bureau of Ordnance. a ae 


The solution obtained was i 


2 Chemical Engineer. Mem. ARS. 5 It was later found that the solution given by Mayer (5) for 
3’ Mathematician. the heating of a composite slab can be reduced to the one given 
* Numbers in parentheses indicate References at end of paper. here. 
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+ beg (b + g)ai 
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Re /p [10] 
he The temperature of the metal is obtained by letting x equal 
lin Equation [6] to yield 
\ N in Qhe ~ + g?)(b sin ail + a; cos ail) (11) 
+ ai{((bg + + + (bg + + 9)} 
A XxX es cit Equation [11] may be rearranged in terms of the dimen- 
sionless parameters X, Y, m, n and 6; to give 
\ 
£4 + n){(1 + + n2Bi2) 
n(1 + m2B;2) + m(1 + 
\ \ \ where 
\ - 
m= \ 
+ — Pm Cm en - 
1 3 | 5 
X= 28 /Q" 


Fig. 1 Dimensionless plot of temperature rise vs. hea 
7 _ for an insulated, infinite metal slab 
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8; = ajl, the positive roots of ~~ 


1 — mnB;? 
(m + 
Equation [12] is symmetrical with respect to m and n; = 50.00 
i.e., the values of the two parameters may be interchanged 0.70 — ey rp 
without affecting the equation. Thus m, the surface-re- aso 
sistance ratio, and n, the ratio of heat capacity per unit area SS 
f the metal to that of the insulation, are mathematically Bi SSS .—~s =i 
of the metal to that of the insulation, are mathematically ; SSS SSS 2003 
= + + + 
Graphical Representation of Solution MJ 
Equation [12] was solved numerically to give as a fune- 0.20 \\ NNSS 
tion of X for values of m and n of 0, 0.5, 1, 2 and 6 in all com- \ \ \\ \\ | \S INQ 505 IN 
binations. These results are plotted in Figs. 1 through 5. t NI \ YN iF 
Each figure contains curves for five values of n for a given Ele | \ \ \\ \ NANI b 
value of m. By double interpolation of m and n, the tem- Zit \ \ Je ee Sie 
perature rise of the metal slab for a given set of conditions 0.09 
can be de termined from these curves. 
The family of curves for m = 0 corresponds to Equation AA AAN SS % —~— 
[21] of Carslaw and Jaeger (1), p. 107 for an infinite slab in one TE 
contact with a well-stirred fluid of finite heat capacity on one 0.05 WLLL LA ULYULN VY 
surface, the other surface being at constant temperature. 004 AVA VA K "6, 
This is the solution to the case where there is no gas-film re- \ \ “el 
sistance. Greebler (2) has applied this specific solution to 0.03 st % NX 
the transient heat-flow applications of aircraft insulation. \\ \ 
Letting n = 0, the family of curves is obtained that corre- \\\ \\\\ \ . 2, \ \ 
sponds to the mid-place lines on the charts for large slabs saad | \ i \\ rk x \ NX \ 
given by Gurney and Lurie (3). This is the solution when 08 % Dk \ 
the thickness of the metal slab is zero. we BacaRsses \ INN \ 
Graphical Representation of Simplified, Approxi- o 2 3 4 5 6 7 
mate Solution X= 
f Fig. 6 Dimensionless plot of temperature rise vs. heating time 
As was pointed out above, the use of the curves in Figs. 1 for an insulated, infinite metal slab 
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through 5 may involve double interpolation. The necessity 
of double interpolation can be avoided with only a slight loss 
of accuracy by the adoption of a new empirical parameter 


Although Y is not an exact function of uw only, a maximum 
variation in Y of only about 4 per cent occurs when y is held 
constant and m and n are allowed to vary over their full 
range of values. Values of Y given by Equation [12] are 
closely approximated by 


where @; are the positive roots of 
1 


Curves of Y vs. X for values of u from 0 to 50 are shown in 
Fig. 6. The use of y instead of m and n is adequate for en- 
gineering applications, and it is recommended that Fig. 6 be 
used rather than Figs. 1 through 5, since only a single inter- 
polation is necessary. 

Notice that Equation [19] may be obtained from Equation 
[12] by setting m = 0, and n = uy in the latter equation. 
Thus, Fig. 6 corresponds exactly to Fig. 1 where n is now re- 
placed by u. 


The following example illustrates the use of Fig. 6. 

The temperature of a flat sheet of steel, 0.15 in. thick, 
initially at 80 F, and insulated with a sheet of asbestos from 
hot gases which are at 3500 F, is to be kept below 800 F for 
60sec. The steel has a heat capacity of 0.15 Btu/lb°F, and a 
density of 0.28 lb/in.*. The asbestos has a heat capacity of 
0.25 Btu/Ib°F, a density of 0.0208 lb/in.* and a thermal con- 
ductivity of 0.13 Btu/hr ft °F. The coefficient of heat 
transfer between the hot gases and asbestos is determined 
empirically to be approximately 125 Btu/hr ft? °F. (a) 
What thickness of asbestos is required? (b) What is the 
temperature of the steel after 45 sec of exposure to the gases 
assuming insulation with the required thickness of asbestos? 


Problem 


Data 
= 0.15 Btu/lb °F p = 35.94lb/ft? 
= 0.25 Btu/Ib °F T’ = 3500 F 
= 125 Btu/hr ft? °F T, = 800 F 
0.13 Btu/hr ft °F =80F 
Im == 0.0125 ft 6. = 0.01667 hr 
= 483.8 lb/ft? 6, = 0.0125 br 


a = 0.01447 ft?/hr 


Solution 


Since both X and yu are functions of the insulation thick- 
ness, part (a) must be solved by a trial-and-error process. 

1) Assume 1 = 0.0075 ft 

2) Calculate 


k mbm 
sam = — = 0.1387 n = 13.47 
= hl pCl 
—T 


= 0.7895 


B=m-+n-+mn = 15.48 Y= 


3) Obtain X = 3.8 from Fig. 6 corresponding to the above 
values of wand Y. (The curve for uw = 15.48 is determined 


by interpolation.) es 
4) Calculate = = 0.0080 ft 


5) Repeat steps 1) through 4) assuming / = 0.0080 ft. 

The calculated value of 1 is now found to be 0.0081 ft. 
This is sufficiently close to the assumed value, and 0.0081 ft 
of asbestos is concluded to be adequate insulation. 

The assumed and calculated values of / converge quit, 
rapidly. An initial assumption of 1 = 0.075 ft, for example 
requires only four repetitions of steps 1) through 4) to obtain 
the correct solution. 

The solution to part (b) is obtained directly from the chart. 
For / = 0.0081 ft 


12.46 
aby /l? = 2.774 


m = 0.1283 n 
pw = 14.19 xX 


Y = 0.84 is obtained from the chart for the above values of 
X and yu from which 


T, = T' —(T' — %)Y = 630 °F 
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Technical Notes 


Hypersonic Conical Shocks for 


Dissociated Air in Thermodynamic 


Equilibrium’ 


The Taylor-Maccoll differential equation is shown to be 
integrable in closed form if the assumption is made that 
the flow in the shock layer (defined as the gas enclosed be- 
tween the shock and the body) is incompressible, i.e., the 
density changes are assumed to be negligible.’ This as- 
sumption yields results that, for dissociated air, agree to a — 
high order of accuracy with the numerical integration of 
The calculated 


SAUL FELDMAN? 


AVCO Research Laboratory, Everett, Mass. 


the complete Taylor-Maccoll equation. 
cone surface pressure coefficients for a wide range of cone 
angles and flight velocities between 10,000 and 26,000 fps, 
are shown to be 3 to 6 per cent larger than the Newtonian 
values at a flight altitude of 100,000 ft. 


Nomenclature 


a = speed of sound 

Cy = 2(p — Po)/poV*.: pressure coefficient 

h = enthalpy 

p= pressure 

u = velocity along a ray 7 

v = velocity normal to a ray 

Vo = (u? + v?)'/2: velocity 

Vw = velocity of the flow immediately downstream of the coni- 


cal shock 
> = ratio of specific heats 
6 = angle streamline makes with the free stream flow direc- 
tion when it emerges from the conical shock 
variable defined by Equation [12] 


» = 
6 = angle streamline makes with the free stream flow direc- 
tion 
& = variable defined by Equation [11] 7 
p = density 
= ray angle 
Subscripts .8 
c = cone 
Newtonian 
s = shock; also denotes an isentropic process blo 
= stagnation 
= free stream 


Introduction 

N ORDER to predict with reliability the heating of a body 

flying through the atmosphere at hypersonic speeds, an 

accurate estimate of the surface pressure is of importance. 

Blunt bodies are being advocated for this flight regime, and 

the Newtonian “theory” for predicting pressures seems to be 
very popular in this regime. 


1 This work was sponsored by the Western ee ery Divi- 
sion, Air Research and Development Command, U. S. Air Force, 
under Contract AF 04(645)-18. 

? Principal Research Scientist. Mem. ARS. 

3 Hayes, in Reference 1, has already demonstrated the useful- 
ness of. the constant density assumption in connection with a 
variety of hypersonic flow problems, including the cone problem, 
although his treatment is different from the one presented here. 


Supersonic flow wean a cone 


Fig. 1 


No information is available in the literature, for flow 
Mach numbers larger than 5.8 in air (2),4 as to the agreement 
of this theory with direct pressure measurements. At Mach 
numbers higher than 6 the dissociation of the gas behind the 
shock becomes of importance. Since for pointed bodies the 
pressure distributions can be obtained theoretically, it be- 
comes of interest to compare the Newtonian values with these 
exact theoretical results. This comparison may also throw 
some light on the accuracy of the Newtonian approximation 
for other body shapes. 

Wedges and cones are the simplest bodies to investigate, 
and in dissociated air, they can become quite “blunt” (in- 
cluded angles up to 135 deg for a cone) before the shock 
detaches. 

For wedge angles of 100 to 120 deg included angle, and 
flight velocities between 10,000 and 26,000 fps, the results of 
(3) indicate that the Newtonian pressure coefficient under- 
estimates the actual value by 10 to 20 per cent. The task 
here will be, therefore, to solve the cone problem. 

The Taylor-Maccoll Equation for a Dissociating Gas in 
Thermodynamic Equilibrium 


The coordinate system used here is given in Fig. 1. The 
conservation equations are 


continuity (pv sinw) + 2pusinw = 0............ [1] 

Momentum udu +uvdv+dp/p =0.............. [3] 

The elimination of v between [1] and [2] gives 

p(u” + u’ cotw + 2u) + =0.......... [5] 


where the accent denotes derivatives with respect to w. The 
speed of sound is given by 


evaluated at constant entropy. This may be done by using a 
Mollier diagram for the gas under consideration, or calculated 
directly by statistical mechanical methods. The Mollier 


4 Numbers in parentheses indicate References at end of paper. 


Eprtor’s Note: 


This section of Jer Proputsion is open to short manuscripts describing new developments or offering comments on 


nr previously published. Such manuscripts are published without editorial review, usually within two months of the date of receipt. 
equirements as to style are the same as for regular contributions (see masthead page ‘of this issue). 
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Fig. 2. Conical shock angles in dissociated equilibrium air 


diagram replaces here the usual equation of state. After use 
is made of [2] and [6], Equation [3] may be rewritten as 


p’ can now be eliminated between [5] and [7] to yield the 
Taylor-Maccoll equation 


(a? — u’2)u” = —a%(2u + u’ cotw) + uu”...... [8] 


The boundary conditions are 


atw = = 0 
and atw = ws, u, = V., [9] 
u’ = —V,,sin (w, — 4) 


where 6 is the wedge angle whose shock angle coincides with 
the conical shock angle w,, at the same free stream velocity 
V.. V,,is the velocity over the wedge. 

For dissociated air behind the shock, V,,, w, and 6 are known 
from (3). If the shock angle w, is assumed to be known, the 
only quantity that remains to be determined is the cone angle 
w,. Once w, is found, all other quantities in the flow field 
‘an be determined. 

Note that in deriving [8] the energy equation has not been 
used. In order to integrate [8] numerically, a relationship 
between the speed of sound a, and the velocity u, has to be 
found. This relationship can be obtained by using a Mollier 
diagram that has constant speed of sound lines, thus relating 
speed of sound and enthalpy at constant entropy. The 
energy equation [4], then relates the enthalpy to the velocity. 
The results of calculations for cones flying through the atmos- 
phere at speeds between 10,000 and 26,000 fps at 100,000 ft 
altitude are given in Figs. 2 and 3. 

For a given cone angle in Fig. 2 the difference in ordinates, 
between any constant flight velocity curve and the one labeled 
w,, gives the shock layer angular thickness. Note that in 
dissociated air the thickness is about half the value it would 
be if no dissociation occurred. 

Fig. 3 presents the cone surface pressure coefficients and 
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Fig. 3 Cone surface pressure coefficients. A comparison with 


the Newtonian value 


compares them with the Newtonian and the undissociated 
y = 14, V. = © values. Here again Newtonian theory 
underestimates the pressure coefficient by 3 to 6 per cent for 
the weak solutions, which are the ones of physical interest. 
Note that the pressure coefficient for a chemically frozen 
diatomic gas y = 1.4 does not represent the behavior of air at 
hypersonic speeds. 


Solution of the Conical Shock Equation for an Incompres- 
sible Fluid Downstream of the Shock* 


Since the conical shock layer thickness in dissociated air at 
hypersonic speeds is fairly small and the flow in the shock 
layer is homentropic, one would expect the air density to have 
only small gradients in the w direction. It is reasonable, 
therefore, to solve the conical shock equation assuming that 


the gas downstream of the shock is incompressible. Equation 
[8] then reduces to 
+ u’cotw + 2u =0.............. [10] 
with the boundary conditions given by [9]. 
By letting 
where 


Equation [10] reduces to Legendre’s differential equation of 

the first degree 

(1 — — + 2E O............ [13] 

5 After the manuscript was completed, the author was in- 

formed that Hayes has used, recently, an approach similar to the 
one presented in this section. 
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where the accents denote derivatives with respect to 7. The 
boundary conditions for [13] are 

[14] 

= V, sin (@, — 8)/(1 — 
The general solution of [13] is 
16 


where A and B are constants of integration determined by 


[15]. They are 
| 1 
2 Ne 
and 
= 
where 
V, Si 
E mal — me). [19] 
SIN Ws 


Equation [14] together with [16, 17, 18] give the following 
implicit expression for the cone angle ab ola 


l 1 + Ne 
In 
2 1 — 7? 


J 7 [ _ Vy sin (w, — 4 \- 1 


2 
A comparison is made in Table 1 of the cone angles obtained 
from the above expression, with the angles found by numerical 
integration of [8] the approximate method seems very satis- 
factory, at least for velocities larger than 10,000 fps. 

The present scheme does not give accurate results for a gas 
with a specific heat ratio of 1.4, regardless of the flight Mach 
number. 

The pressure in the flow field can be determined by inte- 
grating [3] with respect to 7. This leads to a lengthy ex- 
pression in 7 which will not be quoted here. However, a 
simple expression for the cone pressure coefficient will be 


Table 1 Comparison of approximate (analytical) and 
exact (numerical) solutions of the conical shock equations. 
(Altitude = 100,000 ft) 

Flight 

veloc- Cone 

ity, pressure 

Shock Coneangle, coefficient, 

10-3 angle, (degrees) Cre 

(ft/ ws (de- Ap- Cpy 
sec) grees) Exact prox. Exact prox. Cp, 
10 66.76 58.10 58.03 1.5245 1.5239 0.9456 
10 78.87 60.83 61.58 1.7437 1.7480 0.8745 
15 37.78 34.96 34.84 0.6784 0.6774 0.9679 
15 48.07 44.62 44.54 1.0173 1.0164 0.9699 
15 59.73 54.88 54.81 1.3816 1.3808 0.9685 
15 63.48 57.81 57.76 1.4826 1.4820 0.9659 
15 67.14 60.50 60.45 1.5739 1.5735 0.9626 
15 71.14 62.92 62.92 1.6596 1.6596 0.9553 
26 25.42 23.99 23.91 0.3390 0.3383 0.9754 
26 40.03 37.86 37.78 0.7708 0.7701 0.9774 
26 58.89 55.29 55.26 1.3833 1.3828 0.9770 
26 69.96 64.28 64.29 1.6703 1.6700 0.9719 


derived in terms of the known pressure coefficient for the 
wedge that under the same flight conditions will have the 
same shock angle (3). The integration of the momentum 
equation [3], between the shock and the cone angle, gives 


where 
— Pa 
Cp, = 22 
Pu (1/2)p..V.2 [22] 
CP» Pw/ Po and V,/V.. are given in (3), while V./V.. is given 
by 
= 
\ (1 +01 —m) ne elf 


k being given by [19]. 

Again, a comparison of the pressure coefficients calculated 
from the exact solution and the approximate incompressible 
one, Table 1, shows very good agreement. 

For a given flight velocity, it is interesting to note, from the 
last column of Table 1, that the Newtonian pressure co- 
efficient is nearest the correct value at a cone angle inter- 
mediate between 6. = 0 and the detachment value. 


1 Hayes, W. D., “Some Aspects of Hypersonic Flows,” 
Ramo-Wooldridge Rept. MR6, Jan. 4, 1955. 

2 Oliver, Robert E., “An Experimental Investigation of Flow 
About Simple Blunt Bodies at a Nominal Mach Number of 5.8,”’ 
Journal of the Aeronautical Sciences, vol. 23, Feb. 1956, pp. 177— 
179. 

3 Feldman, Saul, ‘Hypersonic Gas Dynamic Charts for 
Equilibrium Air,’’ AVCO Research Laboratory, Everett, Mass., 
Jan. 1957. 
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Addendum to ‘Stagnation Point Heat_ 
Transfer for Hypersonic Flow’! 


MARY F. ROMIG? 


Convair, a Division of General Dynamics Corp., San Diego, 


Calif. 


RATHER unusual oversight on the part of the author 
has apparently caused some confusion in the use of an 
empirical formula for stagnation point heat transfer pre- 
sented in the title paper.? 
This oversight was the complete omission of any physical 
units in the empirical formula. To permit correct use of the 


formula 


Ry 
hae 


the units are as follows 


q is in (Btu/ft? sec) 
Pe Is in (lb/ft?) 


Ry is in (ft) 


The correct units were unfortunately not apparent from the 
analysis presented in the paper. 


1 JeT Propuusion, Dec. 1956, pp. 1098-1101. 


2 Staff Scientist, Scientific Research Laboratory. Mem. ARS. 
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Impact Points of Ballistic 


R. E. ROBERSON! 


- Autonetics Div., North American Aviation, Inc., 


Downey, Calif. 


A method of Singer and Wentworth for calculating the 
impact points of ballistic rockets is generalized to the case 
where the terminal point is at a different altitude than the 
initial point. 


INGER and Wentworth recently described a method for 

calculating impact points of ballistic rockets.2 The 
condition they use in determining range angle is that the 
radial distance from the center of the earth at impact be 
equal to the radial distance R, at burnout. An important 
generalization is when the former has a value R; # Rx, par- 
ticularly if the formulas are to be adapted to other applica- 
tions. 

The notation here follows that of Singer and Wentworth 
except for the introduction of R; above, and of v; for the 
velocity at “impact” (i.e., any chosen terminal condition 
which is specified by an assigned R,) and @; for the inclination 
(generally negative) of v,; to the horizontal. Also Ad¢* is 
used rather than Ad, and ¢* rather than ¢, for total range 
angle and total flight time, respectively, to distinguish this 
case. 

From the Singer-Wentworth paper, the elliptical path is 
described by r~! = GMz(Rwx cos 6)~?[1 + € cos with @ 
measured from perigee. The value of ¢ at the initial instant, 
(Rx), is obtained by setting r = R,. The value of ¢, at 
impact is obtained by setting r = R, and properly interpreting 
the resulting inverse trigonometric function. It is convenient 
to write Ap* = Ad, + Ado, where Ad, =x — (Rx), and 
recognize that ¢; = ¢(R1r) + Ad* = Ado +7. Then it is 
easy to see that « cos Ad. = 1 — (Ruz? cos? 6)(G@Mxz)™ 
(R./R;). It follows that 


2 


R R 29 
1 — ——— |... fl 


By iad i r with respect to @ it is easy to show that 


—Rr 
an 6; € sin [2] 


and by using conservation of angular momentum one obtains 
= cos 0)/(R1 [3] 


These last two quantities, #; and v;, may be of considerable 
interest if the trajectory actually continues beyond the 
nominal “impact” point. 

One further quantity to be determined is the time of flight 
t*. Now, rather thant = [1 — 24,{¢(Rz)}]T., we have 


t* = [1 — — — Ade} ]T........ [4] 


A ease of particular interest is that in which R, differs 


somewhat, but not much, from R,. Let Rr = (1 — v)Rz 
with y< 1. Then Equations [1-3] take the form, correct to 
first order in v 


Received June 20, 1957. 
ia Specialist, Guidance Engineering Department. Mem. 
2 Singer, 8. F., and Wentworth, R. C., “A Method of Calculat- 
ing Impact Points of Ballistic Rockets,” JET PROPULSION, vol. 
27, April 1957, pp. 407-409. 
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6; = —0 + v cos? + tan — cos? [6] 


= — VoL X 


{1 + sin cos + sin? — sin 6 cos 0 GM Riz? cos? . [7] 


Also, using the expression for A, from the reference and 


x — Ad. = — Agi + v = + we find 


— Adz] = + 
ee €le + cos $(Rx)] 
+ 2 [1 + € cos |? 
whence 


vT.V1 -é 


ele + cos 1 


{1+ 2 


One final point worth noting is that the change in latitude 
and longitude between burnout and impact is given incorrectly 
in the reference, as can be seen easily from the case y = 0 
with A, ~ 0. Correct expressions are 


= 


= 


7 aresin (sin A@ sin Y cos A, + cos A@ sin Az) — Az.... [9] 


sin Ad cos 
A longitude = arcsin (= — Agr..... [10] 


@s cos A; 


The difference between these expressions and those actually 
used is not apparent i in the example of the reference, in w aeel 
= O and 2; is close to zero. 


Addendum to ‘Heat Transfer to Satellite 
Vehicles Re-entering the Atmosphere”! 


R. W. DETRA?, N. H. KEMP? and F. R. RIDDELL? 


Aveo Research Laboratory, Everett, Mass. 


ER (1)* a correlation formula for the stagnation point heat 


transfer rate to a blunt body of revolution in hypersonic 
flow was given as 


20, 800 |, 3.25 
Jo (“) (1 ) Btu/ft%sec.... [1] 
~ a/R Po he ha 


Here R is the body nose radius, p and pp the ambient and sea- 
level density respectively, wu the flight velocity, u, the satellite 
velocity taken as 26,000 fps, and h,,, and ha the wall and stag- 
nation enthalpies, respectively. While expressions of the form 
of Equation [1] are not exact results of the heat transfer 
theory on which they are based, they are convenient approxi- 
mate formulas. 

Since (1) was written, further refined heat transfer calcula- 
tions have been made (2) and rates have also been measured 
in shock tubes under simulated hypersonic flight conditions 
(3). It is the purpose of this note to give a somewhat refined 
version of Equation [1] and show how it compares with ex- 
periment. 

The new correlation expression which replaces Equation 


[l]is 
_ 17,600 +15 
bee u sl 
2ge 


Received May 10, 1957. 

' The heat transfer investigation of References (2, 3) on which 
this work is based was sponsored by the Western a ae 
Division, Air Research and Development Command, U. 8. Air 
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2 Principal Research Scientist. 

3 Senior Scientist. 

4 Numbers in parentheses indicate References at end of paper. 
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The subscript 300 means evaluated at 300 K. Equation [2] 
agrees with the calculated heat transfer results in the range sea 
level to 250,000 ft and 7000 to 25,000 fps with an accuracy of 
+10 per cent. It differs from [1] in the numerical coefficient, 
the power of u/u,, and the enthalpy factor. The maximum 
change from Equation [1] for a 300 K wall temperature and 
satellite velocity is a reduction of 15 per cent, and the reduc- 
tion decreases with decreasing velocity. 

Shock tube experiments have been performed in air at the 
Aveo Research Laboratory (2) to measure g,. Both the flight 
stagnation enthalpy and flight stagnation pressure are re- 
produced in these tests. The results are plotted in Fig. 1 as 
q, against flight velocity, with an auxiliary scale of the Mach 
number of the moving incident shock in the shock tube (re- 
ferred to a sound speed of 1129 fps). Three groups of points 
are shown for three different initial shock tube pressures py. 
They represent simulation of flight in three altitude ranges: 
~ = 0.1 em Hg is 111,000-127,000 ft; p: = 1.0 is 64,000- 
80,000 ft; and p, = 10 is 11,000-31,000 ft. The expressions of 
Equations [1, 2] can be converted to the shock tube coordi- 
nates of Fig. 1 by means of normal shock relations and the 
properties of equilibrium air (4, 5). The curves of Equations 
{1, 2] are also shown. It can be seen that particularly at high 
Mach numbers, corresponding to high flight velocities, the 
new correlation is significantly nearer the mean of the data 
than the old one. 

The effect on the results and conclusions of (1) of using the 
new correlation rather than the old one is negligible except at 
surface temperature parameters of over 5000 R-ft'/*. The 
reason for this is that the surface temperature parameter varies 
as the fourth root of ¢,, so the new correlation causes a change 
of at most a factor of 1.042 in that parameter, a factor which 
would hardly show on the figures of (1). Of course, some of the 
algebraic formulas of (1) relating to constant heat transfer 
trajectories will change slightly. In problems other than 
those considered in (1), where the heat transfer rate itself is 
important, the new correlation can make a significant dif- 
ference. 
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Studies in a New Type of Flat Flame 


Burner' 


WILLIAM T. BIEDLER III? and H. E. HOELSCHER? 


The John Hopkins University, Baltimore, Md. 


A new type of device for measuring flame velocities has 
been developed. This device is a burner utilizing a stable 
planar flame front which permits the study of stable flat 
flames, free from contact with any surface. Flame velocity 
as a function of concentration of combustible in the un- 
burnt mixture has been determined for an air-propane 
system. The effect of the geometry of the system on flame 
properties is observed and discussed. Temperature and 
concentration measurements were made in the vicinity of 


the flame and are reported. 


Introduction 


HE process of controlled combustion and an understand- 

ing of this process is one of the challenging problems of 
many current research projects. To study this problem, the 
simplest possible physical system is most amenable to sub- 
sequent mathematic analysis of the results. Hence, it is de- 
sirable to avoid complex geometries similar to those of the 
conelike bunsen flame. This led to consideration of a flat- 
flame burner, representing a geometrically simple one-dimen- 
sional system. 

Methods for measuring flame velocities have been surveyed 
by Fiock (2),4 Linnett (4), and others. The general tenor of 
their remarks is that there is no one method capable of use 
over the entire range of flame velocities likely to be encoun- 
tered (from 3 em/see to 1175 cm/sec) and that reported flame 
velocities vary widely according to the method of measuring 
them. One method of determining flame velocities with ac- 
curacy and precision is that used by Powling (5). A disk 
shaped planar flame front perpendicular to the axis of the 
burner was produced above and not in contact with the lip of 
the burner. This burner was modified by Egerton and Thabet 
(1) who measured flame velocities for various mixtures of hy- 
drocarbons. The general features of their burner are pre- 
sented in Fig. 1. This type of burner is excellent for the de- 
termination of flame velocities of very lean combustible mix- 
tures. No other technique exists for examining flame veloci- 
ties below 15 cm/sec. 

A Powling-type burner basically met the requirements of 
the particular investigation contemplated, i.e., for the pro- 
duction of a homogeneous, laminar, premixed, stationary, 
geometrically-simple flame front. However, closer control of 
flame position and greater simplicity of construction seemed 
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Fig. 1 Flat flame burner of reference1 


desirable. A succession of further modifications led to that 
burner shown in Fig. 2. The flat flames produced by this 
burner met all the previously established criteria of simplicity 
with the added advantages that the flame was free of surface 
contacts and could be probed more easily than in previous 
designs. 


Apparatus and Experimental Procedure 


The apparatus used in this study is quite standard. The 
air supply system consisted of two air compressors, a steel 
calming tank to remove pulsations, reducing and flow control 
valves, and a flowmeter with stainless steel float. The ap- 
paratus was interconnected using } in. ID copper tubing. In- 
strument grade propane 99.9 per cent, a two-stage propane 
reducing valve and a series of needle valves for flow control, 
completed the assembly. 

The burner system included a stainless steel mixing cham- 
ber, the burner itself, shown as Fig. 2, and the burner enclosure. 
The burner was constructed of cylindrical copper tubing and 
was 27 cm tall, 10 em internal diam, and 0.63 cm wall thick- 
ness. The movable plate was a circular copper plate 0.63 
em thick and 17.8 em in diam. The burner enclosure was a 
box 12 in. X 12 in. X 24 in. high made of aluminum with a 10 
in. X 12 in. pyrex glass window for observation and photo- 
graphic studies of the flame. A camera with a telephoto lens 
was used to photograph the flame. Peat ee, 


» 
Experimental Procedure @ re 


To initiate a run, the air compressors were turned on and, 
when a steady res uding was observed on the air flowmeter, the 
burner was lighted using a lean mixture to prevent flashback 
onto the top internal screen. The gas flow was increased until 
a stable planar flame resulted Nhe e the burner. 

Air flow rate, gas flow rate, atmospheric pressure, room tem- 
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Fig. 4 Sample flame photograph 


perature, and plate position (i.e., distance above burner lip) 
were measured for each composition yielding a flame. A 
photograph was then taken of the flame. 

In addition to the above measurements, attention was de- 
voted to three special series of measurements. The velocity 
head was measured above the flame through six equidistant 
pressure taps in the plate. Temperature traverses were made 
both vertically and horizontally in the region between the 
flame and the plate through a slot milled in the plate. Finally, 
samples of the gases issuing from the flame region were taken 
and these were analyzed. 


Results 


In the burner shown in Fig. 2, a flat circular disk of flame 
was caused to float between the burner top and the movable 
plate. The flame is stable, thin, and has well-defined edges. 
The flame may be adjusted to a desired vertical position by 
changing the height of the plate above the burner lip. The 
plate may be moved vertically from a position approximately 
3 cm above the burner lip to a position approximately 6 em 
above the burner lip while maintaining a stable flame. The 
flame velocity apparently remains constant regardless of plate 
position. 

Flame velocities of propane-air mixtures were measured to 
test the performance of this apparatus. The result is shown in 
Fig. 3. The data were obtained by photographing the flame 
under different flow rate and composition conditions. The 
flame area was then computed from the measured diameter. 
All volumetric flow rates were reduced to standard conditions 
of 760 mm of mereury and 18 C to compare the calculated 
flame velocities with those presented in (1). The solid line in 
Fig. 3 is a representation of the data from the work reported 
in this reference whereas the points represent measurements 
from this study. .A sample photograph of the flame is repro- 
duced in Fig. 4. It should be noted that the lower limits of 
stability was found at a concentration of 1.98 per cent pro- 
pane. Below this a stable flame could not be maintained in 
the free space. 
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Fig. 5 Horizontal temperature traverse 


The plate was tapped for pressure measurements at six 
equidistant points outward from the center. Absolute and 
differential pressure measurements were made with a micro- 
manometer (sensitivity +0.0001 in. of water). There was no 
detectible difference in velocity head or impact pressure at 
these six points. 

A radial slot 0.63 em wide and 8.3 em long was milled from 
the edge of the plate toward the center. Visual and photo- 
graphic studies indicated that this slot did not effect the flame. 
It was possible to investigate the temperature and concentra- 
tion of the gas in the region between the plate and the flame 
by downstream probing through this slot. In order not to dis- 
turb the flame, upstream probing was done by inserting a 
probe through the enclosure wall. 

Radial temperature traverses were made at various dis- 
tances below the plate in the region between the plate and the 
flame. These data were corrected for radiation error using 
the method described by Friedman (3). Fig. 5 showsa typical 
horizontal traverse, 2.5 cm below the plate and 1 cm above the 
flame. Fig. 6 shows two typical vertical traverses, the first 
at 2 cm from the burner centerline and the second at 5 em from 
the burner centerline. In both of the latter cases the flame 
was again at 3.5 em below the plate. The theoretical adiabatic 
flame temperature for 2.16 per cent propane-air mixture is 
1332 C (this is the gas composition used for all of the data 
shown in the figures and photograph accompanying this 
paper). It is observed that the corrected measured tempera- 
ture at the flame front is very close to this theoretical value. 

The concentration measurements of the burner exit gases, 
made on a gross averaged sample as well as on point samples 
taken between the flame and the plate, indicate that there is 
no unburned propane and no carbon monoxide in the off 
gases. This conclusion is believed to be accurate to within 
+0.03 per cent and is in agreement with reports of previous 
investigators. 

It was suggested that the plate temperature may be exerting 
significant influence on the flame properties. As indicated in 
Fig. 6, the plate temperature was undoubtedly well above 
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500 C. A new plate was constructed with provision for water 
cooling. The plate was hollow and a water stream was in- 
jected into the cavity tangentially and removed axially. 
The flow rate was high enough that the temperature increase 
between inlet and exit stream was less than 2 C. Using this 
water cooled plate, the data taken with the solid plate were 
rechecked. In all cases the data were reproduced exactly. 
The effect of a heated plate remains to be considered. 


Conclusions 


The flat flame burner described is suitable for the measure- 
ment of flame velocities in lean mixtures with a high degree of 
accuracy and sensitivity. These measurements may be made 
in this device under conditions where surface effects are known 
to be absent. This device is limited to the study of slow burn- 
ing flames. 

At the present time, relationships between the following 
parameters have been studied: — 

(a) Plate height above the burner lip. _ 

(b) Flame height above the burner lip. 

(ec) Concentration of combustible in the inlet gas mixture. 

(d) Flame velocity. 

Thus far, only flame velocity and the concentration of com- 
bustible in the unburnt mixture seem to be directly related. 
Egerton and Thabet (1) state that for a given gauze height, 
Fig. 1, the product of the flame height above the burner and 
flame velocity is approximately constant. No satisfactory 
correlation of this type has been found in the present study. 
Over a wide range of condition, the flame height was observed 
to be solely a function of the plate height, all other variables 
being constant, since it permits the flame to be moved (for 
example, relative to a stationary thermocouple) without in- 
troducing extraneous effects on flame velocity. 
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Comments on the Powered Flight 
Trajectory of a Satellite 


THEODORE N. EDELBAUM! 


United Aircraft Corporation, E. Hartford, Conn. a 
Nomenclature 


= acceleration of gravity at the earth’s surface 
= radius of the earth 

horizontal component of velocity 
vertical component of velocity 
time 

burnout time 

altitude 

distance of closest approach of satellite orbit 
ratio of thrust to mass 

angle between thrust vector and horizontal 
= Lagrange multipliers 


= 
ll 


> 
> 
& 
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Discussion 


FEI YWO analyses of the optimum thrust angle variation of a 

satellite launching vehicle have appeared in JET PRoPUL- 
sion (1, 2).2 Both of these analyses have shown that the 
tangent of the angle between the thrust vector and the hori- 
zontal should vary linearly with time. However, neither of 
these analyses includes the effect of centrifugal force, which 
is significant for satellite launching vehicles.* 

The significance of centrifugal force can be demonstrated 
by starting with the equations of motion for a round earth and 
noting what assumptions are necessary to reduce these equa- 
tions to the flat-earth equations assumed in (1, 2). Assuming 
the vehicle to be a point mass moving in a vacuum around a 
spherical, nonrotating earth, the equations of motion may be 
written, in the notation of (1), as follows 


(R, + 


= Nsiny — 


U,U,y 


= |. | 


These equations may be simplified by assuming that the 
altitude during powered flight is small compared to the earth’s 
radius. This assumption allows R, + y to be replaced by R, 
in Equations [1] and also allows the second term in the second 
of Equations [1] to be dropped. This latter simplification may 
be seen by rearranging the second of Equations [1] with the 
aid of the third of Equations [1] 


Uy R, 
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As the altitude is always small compared to R., dy/R, may be 
considered infinitesimal and can be eliminated. With these 
considerations, Equations [1] reduce to Equations [3] 


The equations of motion [3] are equivalent to those used in 
(1, 2) except for the centrifugal acceleration term in the first 
equation. The centrifugal acceleration u,?/R, may be 
eliminated only if it is small compared to the gravitational 
acceleration g or to the vertical component of the thrust to 
mass ratio N sin y. As these assumptions are not usually 
satisfied for satellite launching vehicles, centrifugal accelera- 
tion must be retained in the equations of motion. The follow- 
ing is the derivation of the optimum path including this term. 

The optimum launching trajectory of a satellite, for the 
equations of motion [3], can be found by use of the Calculus 
of Variations. For convenience, the problem is formulated as 
minimizing the time required to reach an orbit with a given 
distance of closest approach. This is equivalent to the prob- 
lem of maximizing the distance of closest approach for a given 
burning time. The variational integral for minimizing time 
for a given functional relationship between altitude, horizontal 
velocity, and vertical velocity at burnout takes the form 


du, 
r= fas) —wN +0) + 


d 
de (' — Ncos v) + Xz dt...... (4] 


The Euler equations for this integral take the following form 


dy, 

(5) 

dt 

Ey = 0 = Neos + AN sin 


These equations may be operated upon to yield _ 
Ai ) 
an 
A = Air + Asr(T — t) 


2r 
= Aor + f —t) 


44 


The values of yr, V7, and Vy, must be such that the re- 
sulting orbit has the required distance of closest approach. 
This can be insured by imposing a transversality condition for 
the distance of closest approach on the ratios of the burnout 
value of the Lagrange multipliers of Equations [6]. For the 
present problem, the transversality conditions are conditions 
(7], where D, := 0D/dy ete.! 
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‘ The expression for distance of closest approach may be found 
in Reference 1. 
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The solution of Equations [3] and [6] with the initial 
boundary conditions and the terminal conditions [7] is the 
solution of the problem for any given variation of thrust and 
mass with time. Unfortunately, the necessary inclusion of 
centrifugal force does not allow any simple closed form solu- 
tion for the optimum path. 
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Some Decomposition Products of — 
Normal Propyl Nitrate 


S. A. GREENE! and L. J. GORDON? 


Aerojet-General Corp., Azusa, Calif. 


ERFORMANCE calculations based on equilibrium con- 
stants were made for normal propy] nitrate as a mono- 
propellant. Table 1 shows a calculated composition and re- 
sultant flame temperature. 
When the material is decomposed in a motor, the flame 
temperature as measured with an unshielded thermocouple is 
found to be 2100 + 100 F, considerably higher than the cal- 
culated temperatures. Gas samples, removed from the motor 
by a water-cooled probe just upstream of the nozzle, were 
analyzed by the technique of gas chromatography. The re- 
sults are shown in Table 2. Not recorded, but qualitatively 
identified, were HCN, NH3, COs, and C,.). No attempt 
was made to measure these quantitatively since it was thought 
that reaction between all components except the last would 
take place in the sample bombs, and representative gas sam- 
ples would not be obtained. 
It can be shown for this system that a mass balance cannot 


Table 1 Equilibrium Calculation for Decomposition of 
Normal Propyl Nitrate with Flame Temperature of 1880 F 


at 300 psia 
? Component Per cent volume 
43.7 
CO 39.3 
Ne 7.6 
CH, 
Cw) 
1.5 


Table 2 Analyses of Some Decomposition Products of 
Normal Propyl Nitrate 


Per cent volume 


Gas 1 Ps 3 4 Average 
He 19.1 21:0 21.4 21:5 20.8 
CO 49.2 48.2 50.7 48.9 49.2 
Ne 9.6 7.8 7.5 Ce 

CH, 13.6 16.7 16.7 15.7 15 
2.5 O7 1.0 LZ 

6.0 5.6 | 
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be made when using the analyses of Table 2 or any combina- 
tion of the qualitatively identified products. Such products 
as NO, N.O, or (C:N:) would permit a mass balance, but the 
NO was looked for and not found; qualitative analyses for the 
latter two products were not attempted. Some flame-tem- 
perature calculations were made on the basis of a postulated 
analysis composed of the experimentally analyzed products 
and (C,N:). The calculated flame temperature was found to 
be 400 F higher than the observed flame temperatures. 

Although this work does not permit quantitative conclu- 
sions, indications are that even at 2200 F, performance cal- 
culations based on equilibrium constants must be considered 
suspect. 


Thermal Conductivity and Prandtl 
Numbers for Nonpolar Polyatomic 
Gases at High Temperatures 


ve 
ELDON L. KNUTH! 
University of California, Los Angeles, Calif. 
Nomenclature 
A* = ratio of Schmidt number for elastic-sphere model to 
Schmidt number for given intermolecular-force model 
Cp = specific heat at constant pressure vt bay 
D = coefficient of self-diffusion 
Pr = = Prandtl number 
Se = u/pD = Schmidt number 
= coefficient of thermal conductivity 
= coefficient of viscosity 
p = density 
Superscript 


’ = refers to the property of a monatomic gas having the same 
molecular weight and the same intermolecular force 


constants as the polyatomic gas 


HE purpose of this note is to present, in a form convenient 

for use by aerodynamic and propulsion engineers, expres- 
sions which permit the calculating of the coefficient of thermal 
conductivity and the Prandtl number for nonpolar polyatomic 
gases at high temperatures, starting only with the specific 
heat of the polyatomic gas and expressions for properties of a 
monatomic gas. 

Molecular transfer of heat in a polyatomic gas is considered 
to take place by transfer (through collisions) of the molecular 
energy associated with translation and by transfer (through 
self-diffusion and inelastic collisions) of the molecular energies 
associated with rotation and vibration (1)? p. 238. If one as- 
sumes that the transfer of translational energy and the self- 
diffusion of molecules are not affected by the presence of in- 
ternal degrees of freedom, and that the transfer of energy from 
molecule to molecule is so rapid that an equilibrium (Max- 
wellian) distribution of molecules corresponding to the local 
temperature is realized even though a temperature gradient 
exists, then the following expression for the coefficient of ther- 
mal conductivity of a polyatomic gas may be derived, see (1) 
p. 238 or (4) p. 501: 


A =X’ + pD(cy — cp’) 


The first term in the right-hand side corresponds to transfer 
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Fig. 1. Prandtl number vs. specific heat for nonpolar gases 


of energy associated with the translational degrees of freedom 
whereas the second term corresponds to transfer of energy 
associated with the internal degrees of freedom.* Since it is 
assumed that the presence of internal degrees of freedom has a 
negligible effect on the coefficients of viscosity and diffusion, 
this equation may be written in the following form also 


Pr’ £5: — <p" 
Se’ 


The value of the Prandtl number for the polyatomic gas may 
be obtained directly now from the relation 

Using results of molecular theory of gases, several suggestions 
for bringing these equations into agreement with experimental 
results have been made. The various suggestions may be con- 
sidered equivalent to assuming a value of 3 for Pr’ and vari- 
ous values for Se’. Eucken (2), without the benefit of much 
of our present knowledge, assumes the equivalent of Se’ = 1; 
Chapman and Cowling (1) p. 238, mention that Sc’ = 0.83 
(result for Bernoulli elastic-sphere model) and Sc’ = 0.64 
(result for Maxwell repulsion model) would give conductivity 
values of the correct order of magnitude; Vines and Bennett 
(7) subscribe to Sc’ = 0.83/A*, where A*, a measure of the 
deviation from results for the elastic-sphere model, is equal to 
1.10 for a wide range of temperatures (result for Lennard- 
Jones attraction-repulsion model), so that Sc’ = 0.76, whereas 
Hirschfelder (5) p. 25, speculates that one should use perhaps 
Pr’/Se’ = 1, implying that Sec’ = 0.67. The resulting equa- 
tions for Prandtl number as a function of specific heat, as well 
as data for 15 nonpolar gases, compiled by Keyes (6), Hilsen- 
rath and Touloukian (3), and Vines and Bennett (7), are 
presented in Fig. 1 and the table. For each gas, the Prandtl 


Pr = Pr 


Values of Prandtl number and specific heat for 15 Non- 
polar gases 


Temp., 
Gas Pr Ref. 
1 He 260 1.00 0.70 3 1.13 
2 A 572 1.00 0.65 6 1.10 
3 He 800 1.42 0.66 3 1.12 
4 CO 260 1.41 0.72 3 1.10 
5 Ne 752 1.47 0.69 6 ee | 
6 NO 260 1.44 0.73 3 1.10 
7 Oz 260 1.45 0.70 3 1.10 
8 N20 212 2.01 0.72 6 1.10 
9 CO. 620 §=2.28 0.68 3 1.09 
10 CH, 212 1.88 0.75 6 1.09 
11 C3Hs 284 4.62 0.75 7 1.10 
12 C2H¢ 284 0.76 7 1.09 
13 320 7.86 0.76 7 1.10 
14 320 5.90 0.76 7 1.10 
15 6068.31 0.76 7 1.10 


Note: Values of A* are determined from Table I-N of 
Reference 4 using values for the maximum energy of at- 
traction tabulated in Table I-A of same book. 
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number used corresponds to the highest temperature for which 
data are given. Examining the table, it is seen that the cal- 
culated values of A* substantiate the use of A* = 1.10 forthe | 
Lennard-Jones attraction-repulsion model. Examining the 
figure, it is seen that the presented data agree best with the 
curve for the Lennard-Jones attraction-repulsion model, es- 
pecially for the gases having many internal degrees of free- 
dom. 

Since interactions of polar molecules are not described ade- 
quately by the simple interaction potentials, such as the Len- 
nard-Jones potential, one would not expect data for polar 
molecules to agree with the presented theoretical curves. attitude angle must be a linear function of the time. 

Water vapor is the polar gas encountered most frequently per- This result may be generalized as follows: Taking horizontal 
haps by aerodynamic and propulsion engineers. Also, since and vertical axes Ox, Oy, respectively, through the launcher 
the assumption of large average distances between molecules, and assuming the trajectory to lie in the plane Ozy, let 
in comparison with molecular diameters, is fundamental to (41,92) be the rocket coordinates at “all-burnt” and let (q:,g2) 
the molecular theory of transport phenomena in gases, one be its velocity components. Then, if # is any function of the 
would expect experimental and theoretical results to agree variables (q1,92,91,92), (€.g., rocket range, maximum altitude or 
better, for a given pressure, at high temperatures than at low eccentricity of subsequent orbit), R will take a stationary 
temperatures. These expectations are substantiated by value (usually a maximum) for a thrust direction program 


Optimal Rocket Trajectories _ 


D. F. LAWDEN! 


Canterbury University College, Christchurch, New 
Zealand 


RIED, B. D., has proved (1)? that, if air resistance and 
variations in the magnitude and direction of gravity dur- 
ing the thrust period of a satellite vehicle can be neglected, 
then, for optimum performance, the tangent of the thrust 


available experimental data, not presented in the accompany- determined by the equation 
ing figure, see (7). 
oR. 
It is concluded that the equations 7 39 (T —t) + — 
tan y = 

an  2ep where y is the angle made by the direction of thrust with Or 

rd oe _ sila 3 cp’ d at time ¢ after launching and T is the total period of thrust. 
Thus tan y is always a bilinear function of the time, viz. 
provide the best available analytical approximations for the eae 
coefficient of thermal conductivity and the Prandtl number | tan yy = — 
for nonpolar polyatomic gases at high temperatures. It is be- ; +d 
lieved that these relations are in a form convenient for use by ane proof will be found in (2). - 7 
aerodynamic and propulsion engineers, Clearly, tan y is a linear function of ¢ if, and only if, R is 

Subsequent to the submission of the manuscript for this References 
technical note, the author discovered that Hirschfelder has 1 Fried, B. D., “On the Powered Flight Trajectory of an 
superseded his speculation with a detailed study (Journal of Earth Satellite,’ Jer Proputston, vol. 27, June 1957, pp. 641- 


Chemical Physics, Feb., 1957). He concludes, as is done in the 643. : : 
present note, that the best analytical expression for the ther- 2 Lawden, D. F., “Dynamic Problems of Interplanetary 
mal conductivity of a nonpolar polyatomic gas is obtained Flight,” Aeronautical Quarterly, vol. 6, 1955, pp. 165-180. 
using results for the Lennard-Jones model (or for the similar 
Buckingham model). The author believes that the present 
publication will hasten acceptance of the modified Eucken 
equation by aerodynamic and propulsion engineers. 


Received August 15, 1957. 
1 Professor of Mathematics. _ 
2 Numbers in parentheses indicate References at end of note. 
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loge (Mo/M,) is a monotonically decreasing function of », the 
thrust factor for circumferential thrust. This penalty in 
mass ratio for low acceleration is caused by leaving the work- 
ing medium at a higher orbital energy than corresponds to the 
altitude of the original orbit. There seems to be a way to re- 
duce this penalty in mass ratio to an arbitrarily small value if 
one is willing to accept a penalty in time instead. 

It appears advantageous to make the escape trajectory a 
succession of short periods of acceleration near the perigee and 
long periods of coasting. If there is no acceleration except 
tangential acceleration near the perigee of an elliptical orbit, 
the altitude of the perigee will not change greatly. In the 
limit as the period of thrust approaches zero, i.e., a pulse at 
the perigee, the trajectory will approach a succession of ellip- 
ses with the same perigee. This result comes from the consid- 
eration that the orbit of the vehicle must be elliptical, both 
before and after the pulse. Since both before and after the 
impulse the velocity is normal to the earth radius, the vehicle 
must be at the perigee of both orbits. In this case the work- 
ing fluid will be given off with the same distribution of orbital 
energies as in the impulsive case and the mass ratio corre- 
sponding to impulsive accelerations to escape velocity can 
be achieved. 


The a a Rocket with 
Tapered Exhaust Velocity 


HOWARD S. SEIFERT 


7 The Ramo-Wooldridge, Los Angeles, Calif 


The effect of varying the exhaust gas velocity of a rocket 


in an earth observer’s frame of reference is examined. 


their burnout velocities if they consume amounts of energy 
which are equal, but distributed in the first one in such a way 
as to cause the exhaust trail to have zero (or small) velocity 
relative to ground, and in the second in the usual manner 
with a fixed velocity relative to the rocket? Presumably the 
first distribution should be more efficient, since more of the 
propellant energy goes into the payload. 

b. Suppose we are given a fixed amount of energy per 
unit payload which may be distributed uniformly into any 
desired fraction ¢ of the mass of a rocket in order to propel 
the remaining payload fraction 1 — ¢. What is the value of 
¢ which will result in the highest terminal velocity of the resi- 
due? The possibility of optimum ¢ is intuitively evident, 
since with fixed energy, if ¢ is too low not much momentum 
will result per unit energy expended while if ¢ is too high the 
consequent low exhaust velocity will reduce performance. 

In both the above problems it has been assumed that any 
arbitrary amount of energy may be associated with a unit 
mass; in other words, temperature is not a limiting factor. 
This is not true at present, but there is no a priori reason to 
believe that the containment of high temperature matter is 
an unsolvable problem. 


’ The Tapered Exhaust Velocity Program 


7 The equation of motion of a rocket moving vertically, and 
possessing at ¢ = 0 an initial exhaust velocity co which in- 
creases by the same amount that the vehicle speed increases, 


is 
- (Mo — mt)(a +9) =m adt+ {1} 


where 


mass flow rate (constant) 


3. 
Il 


so that the exhaust gases possess little or no kinetic energy _ Mo = initial total mass 


Under certain boundary conditions an improvement in_ 
terminal velocity of 10 to 20 per cent is obtained relative to — 


a rocket with fixed exhaust velocity. Given a fixed budget q 


of energy per unit payload mass, the optimum mass of 
working fluid into which to distribute this energy (for 
maximum burnout speed) is calculated as a function of 
gravitational field strength. 


Introduction 


CONSIDERABLE fraction of the propellant energy of a 
rocket is dissipated in the exhaust which, even if ex- 
panded with 100 per cent thermodynamic efficiency, trails 
along after or away from the rocket at high speed, carrying 
energy which would preferably have been retained in the pay- 
load. The conventional chemical propellant must be ejected 
at constant speed with respect to the rocket, a relatively in- 
efficient process. Recent general interest in the use of nuclear 
energy for rocket propulsion suggests the possibility that the 
exhaust velocity of a rocket may become a controlled variable 
rather than a constant as it has been so far. This paper will 


calculate the relative performance of rockets with tapered- — 


velocity exhaust and fixed-velocity exhaust under certain 
boundary constraints, and show that improvement in payload 
velocity can be obtained. An additional calculation will be 


made to show certain optimum circumstances under which a — 


constant but arbitrarily adjustable exhaust velocity may be 
used. 


Statement of the Problems 
Two problems will be examined: 7 ah 


a. Given two rockets with identical mass ratios, mass 
flow rates, and propellant masses, what will be the ratio of 
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This may be reduced by differentiation to 


upward acceleration relative to ground 
g = constant gravitational field acceleration 
co = initial (constant) exhaust velocity 


where to = Mo/m. If the initial acceleration is (cgm/Mo) 
— g], then Equation [2] may be integrated twice to give 


where v’ is the vehicle velocity upward. At burnout 


t, = €Mo/m = (to 


~ 


where 


R-1 M, 
= = 1 
R + M, oading fraction 


= = mass ratio ar 


The velocity at burnout (¢ = ¢,) follows from Equation [3] 


Let us define a gravity loss parameter y, = gt,/c, which is 
essentially the decrease in missile speed at burnout caused by 
the gravitational deceleration measured in multiples of an ex- 
haust velocity c. This exhaust velocity c is that of an equal- 
energy rocket of constant-exhaust velocity. It will presently 
be related to cy. Using y, in Equation [4], we arrive at 


ve! = R — 1) (1 “) 5] 
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The velocity reached by a conventional rocket with the same R 
and burning time ¢, is 


v, = clhn R — gt, = R — yo).......-- [6] 


We may then express the ratio of the burnout velocities of 
the tapered-c vs. constant-c vehicles as 


c(In R — ¥5) 
Re (R 1) ] E R+ ‘| 7] 


The two rockets must expend equal energies, as well as 
have the same values of Mo, m, R, and t,; co and c are therefore 
related. The energy equivalence may be expressed as 


—m = - 
2°" ~ 9 Jo 


where 
c = exhaust velocity of constant-c rocket 
m, = total propellant mass of either rocket 
m = constant mass flow rate of either rocket i 
v’ = instantaneous velocity of tapered-c rocket 


Substituting Equation [3] in [8] and simplifying, after ex- — 

tensive algebra, results in iy 
1 — ;,2/12 
R 

We may now eliminate ¢o/c from [7] by the use of [9] and 
have the comparison v,’/v, in terms of the two fundamental 
parameters R and y, 


ve’ 1 R-1 
1 2/12 — . {10 


A typical value of y, = gt;,//c is 0.25. In efficient rocke ts, it 
is seldom that y, > 1. isietote * we may neglect y,2/12 
relative to unity and simplify [10] to 
R —1)R-**—y 
Vp nR-y 


Values of v,’/v, are shown in Table 1, and plotted in Fig. 1, 
indicating the effect of the tapered-exhaust velocity on v,’/vs. 
It can be seen that in field-free space (yo = 0) at R = 5, an 
11 per cent improvement in relative terminal velocity can be 
obtained. In the presence of gravitational fields, the ratio is 
increased, although both velocities decrease, of course, in abso- 
lute magnitude. 


Maximum Temperature of Tapered Exhaust 
If the exhaust stream is produced by thermal expansion, the 
temperature of the gas reservoir will be proportional to the 
square of the exhaust velocity relative to the rocket. The 
maximum temperature will be reached at burnout in the 
tapered-c case. The ratio of maximum temperatures T,,,’/T,,, 
tapered to constant operation, will be 


+ co)? 

Let us use relations [5] and [9] to evaluate the field-free 
case (y, = 0) 


leo(R — 1) + co}? _ 


[13] 


m 


For example, we see from Table 1 that an 11 per cent im- 
provement in v,(R = 5) results in a fivefold rise in maximum 
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Fig. 1. Effect of mass ratio R and gravity loss 7, on velocity 
ratio vi’ /Vb 
temperature. If Equation [12] is evaluated for R = 5 and 


= 0.25, the ratio T,,’/T,,, becomes 4.4. 
Maximum of Tapered Exhaust 


’ of the tapered-c eal 


Velocity 
~The maximum exhaust velocity c,, 
relative to the fixed exhaust c is given by 
, 
+ Co 


CO 


C max 


Table 1 Effect of R and y») on burnout velocity ratio 
) 
a 0 0.1 0.25 0.5 1.0 
R (vp’/vy) 
1.5 1.005 1.008 1.015 
2.0 1.020 1.023 1.032 1.072 . 
a0 Ti 1.12 1.13 1.16 1.29 
10.0 1.24 1.25 1.27 1.31 1.42 
20.0 1.42 1.43 1.46 1.50 1.62 
_ * For these large values of yo, both v,’ and w are negative; 
i.e., the rocket is sinking rather than rising at burnout, and 
the ratio, v,’/v is not significant. 


Using Equations [4] and [9] and the fact that gt, = cys 


leads to 


iz) 
3 


= — 72/12 — (14] 
Thus for R = 5 and yz = 0.25, c,,’/c 

WR = 2.24. The total span of exhaust velocities 
Cm'/Co = 3.7; whereas if ys = Oc,’/co = R = 5. This indi- 
cates that co is substantially lower in field-free space tha in on 
the earth’s surface. 


Comments on Tapered Exhaust Rocket 


The initial exhaust velocity co must be determined by: equi at- 
ing total energy in the propellant in the uniform and tapered 
programs. If it is desired that the tapered exhaust be abso- 
lutely motionless with respect to an earth observer, it is neces- 
sary to start off the rocket with a boost velocity equal to ¢p. 
This will modify Equation [4] and all subsequent equations in- 
volving [4]. 

The preceding calculations are of interest primarily 
situations where energy is limited in quantity and there is no 
penalty associated with high rate of release of energy (i.e., 
power or temperature not limiting). This is, in fact, the re- 
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verse of conditions governing present chemical rockets. How- 

ever, there are propulsion techniques, such as accelerated par- 

ticle beams, in which the controlling factors may differ from 

those of present conventional engines. 


Optimum Mass Ratio With Fixed Energy Per Unit Payload 


Let us assume that a fixed amount of energy E/m per unit 
payload mp» can be distributed over a propellant mass m, and 
that a constant mass flow rate m and exhaust velocity ¢ will 
be used. What mass ratio R = (m, + mo)/mo will provide 
maximum terminal velocity for mo? 

If the energy £ is distributed uniformly over the working 
fluid mass m,, then the ideal exhaust velocity ¢ is given by 


Inserting this c into the conventional velocity Equation 


[6] gives 
1\ 0.5 
My mo 


We may express this v, in terms of the fundamental quan- 
tities mass ratio R and initial acceleration ratio a = ao/g by 
means of the following relations: 


Mp + Mo 


mo 


R= 


Table 2 Take-off acceleration vs. optimum mass ratio 
ao (™ + 
a= — Rn = {| 
mo opt 
0.0* 12.5 
0.5 9.7 
4.0 6.0 
9.0 5.4 
5.0** 
*qa = 0 means no initial acceleration (rocket momen- 
tarily balanced against g). 
** Field-free case. 


0.5 1 \0-5 0.5 1 \0-5 1 0.5 
mo Mp mo R-1 


(7 — 1\ (my, + mo) 
= (mp + me) 
Glo g R c 


(7 ge 
gt, = 
R 


where F = constant thrust, and ap is the initial acceleration in 
a direction opposite to the gravitational field, with a = ao/g. 
Substituting [19] and [18] into [16] 


We shall now examine Equation [20] to see if v, has a 
maximum with respect to R. The available energy per unit 
payload E/m is presumed to be specified. 


Calculation of Maximum v, With Respect to R 
Taking the derivative of [20] with respect to R, and setting 
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it equal to zero results in 


R, —1 1 2 
——— {1 — — } ]..... 21 


where 2, is the mass ratio for maximum v;. We note that in 
field-free space, @ = d)/g approaches infinity, and that for 
takeoff from earth with 1 g upward acceleration a = 1.0. 
Solving the transcendental Equations [21] numerically yields 
the results shown in Table 2 for the relation between a and 


Conclusion 


We see from Table 2 that in field-free space a rocket with 
mass ratio 5 will give greatest terminal velocity to m5 per unit 
of energy available. To cite an extreme example, (Reference 
1), if mp = 1 ton, and E is the energy available from the fission 
of 3.75 lb of uranium, then the v, available by distributing this 
E into the fission fragments only, could be multiplied more 
than fortyfold by distributing this same £ instead into 4 tons 
of inert working fluid. 


Reference 


1 Seifert, H. S., and Mills, M. M., Physical Review, vol. 71, 
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A Flow Reactor for High Temperature 
Reaction Kineties' 


L. CROCCO,? I. GLASSMAN? and I. E. SMITH 
Department of Aeronautical Engineering, Princeton Uni- 
versity, Princeton, N. J. 


HE purpose of this note is to describe a new type of flow 

reactor designed for the study of the reaction kinetics of 
propellants at high temperatures. In principle the apparatus 
consists of a tubular duct in which a high temperature stream 
of inert gas flows. A secondary cold stream of monopropellant 
or fuel-oxidizer mixture is introduced into and mixed rapidly 
with the inert gas, thus conditions within a rocket motor com- 
bustion chamber are simulated without the complication of 
recirculation. Downstream of the injection point reaction 
commences and continues until all the reactants have been 
consumed. By varying the velocity of the inert carrier gas 
the reaction zone may be extended over any desired length. 

By introducing only a small quantity of the reactants com- 
pared with the mass flow of carrier the decrease in tempera- 
ture resulting from the mixing of the two streams is small, and 
likewise the increase in temperature due to the exothermic 
reaction will also be relatively small. Endothermic decom- 
position reactions which give a subsequent decrease in tem- 
perature may be studied by this scheme, as well. 

Because the zone in which the reaction takes place is large 
and the temperature gradient small, such factors as the heat 
transfer to the unreacted gas and the back diffusion of active 
species, which are dominant effects in self-propagating flames, 
are rendered negligible. In this way an unambiguous deter- 
mination of the reaction kinetics in the presence of inerts is 
possible. 

This scheme has been used to determine the rate of decom- 
position of ethylene oxide with either nitrogen or helium as 


1 This research was supported by the U. 8. Air Force through 
the Air Force Office of Scientific Research of the Air Research and 
Development Command under Contract AF 18(600) 1527. Re- 

roduction in whole or part is permitted for any purpose of the 

. 8. Government. 
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Fig. 1 Cross section of mixing nozzle 
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the carrier gas. Compar: able results have been obtained 


demonstrating the independence of the reaction on the physi- 
cal nature of the inert gas. 

The duct employed is a stainless steel tube 2 in. in diam and 
approximately 6 ft long. In order to maintain a constant 
temperature along its length the duct is heated electrically, 
and the wall temperature may be determined from the power 
input. The inert gas is he ated to the same te mperature by 
passage through a packed ceramic bed which has previously 
been heated by a fuel oil burner. A fine control over this tem- 
perature is obtained by metering a certain amount of cold 
carrier gas into the heated stream. 

The test fuel injection section of the duct is shown in Fig. 1, 
and consists of a sharply convergent section followed by 
parallel throat 2 in. in diam and about 14 in. long. Four small 
bore injection sale branching from a common manifold carry 
the fuel and are arranged so that the streams issuing from 
them impinge on one another, thus a region of high turbulence 
is created. The gas mixture is then brought back to the duct 
section by means “of the 15 deg diffuser. 

The course of the reaction is followed by means of an axi- 
ally-traversing probe carrying at its tip a platinum-platinum/ 
rhodium thermocouple which is 0.001 in. in diam. The output 
from this thermocouple is balanced by a potentiometer for 
zero injection of fuel and during an experiment the out of 
balance EMF due to the change in temperature as the fuel 
reacts is recorded on a sensitive recording potentiometer. A 
single longitudinal transverse of the probe provides a record of 
the temperature change through the reaction zone, from which 
the reaction rate at any specific temperature covered may be 
determined. The determination of the reaction rate over a 
wider temperature range is made by carrying out a series of 
such determinations at different initial temperatures. 

Evidence for the complete mixing of the injected and car- 
rier streams within the nozzle has been obtained by injecting 
cold nitrogen into the hot carrier stream. By employing a 
special probe to enter the throat temperature records of the 
mixing of the streams have been made and these demonstrate 
that a practically homogeneous mixture is obtained in be- 
tween 0.5 and 3 millisecs according to the flow rate of carrier 
gas. 

A temperature record which covers a range from 1020 K to 
1100 K obtained for a typical experiment is given in Fig. 2 
The mean velocity of the gas stream was 1520 cm/sec cor- 
responding to a time scale of 66 millisecs per meter of duct. 

The rate constant for the ethylene oxide reaction has been 
calculated assuming first order kinetics (1)5 from the following 


5 Numbers in parentheses indicate References at end of note. 
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Fig. 2. Typical temperature-distance record for ethylene oxide 
decomposition reaction-range 1020 K to 1100 K 
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Fig. 3 Plot of ethylene oxide decomposition reaction 


rate data 


relationships 
de 
dt 
cQ = CT; — T) 
E/RT 


(7; 


where c is the concentration at any point, 7’ the corresponding 
temperature, 7’; the final — rature, Q the heat of reaction 
for the fuel-inert gas mixture, C, a mean specific heat, ¢ the 
time, & the reaction rate constant, R the gas constant, A the 
collision frequency, and E the activation energy. 

A plot of log k against 1/T is shown in Fig. 3 for a typical 
experiment. The points are best represented by a straight 
line indicating the correct assumption of an over-all-first- 
order reaction mechanism, and the slope corresponds to an 
activation energy of 42-44 keal/mole. This value is sig- 
nificantly lower than that obtained by Mueller and Walters 
(1), who found 52.7 keal/mole over a temperature range of 
623-713 K, though Lossing, Ingold, and Tickner (2) also 
found 42 kcal/mole at 1123 K. 

The complete results of this investigation, including in- 
formation on rates of other gaseous systems and on the in- 
duction or ignition lags of both gaseous and liquid fuels as a — 
function of temperature, will be published at a later data. — 


—T) dt 
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HOW THE 13,200-TON 
PRESS BOOSTS EXTRUSION 
CAPACITY 


FROM 11” TO 28” 


STRUCTURAL SHAPES 


FROM 9” TO 20” 
STIFFENED SECTIONS 


FROM 10” TO 24” O.D. 7 
TUBING 


DOW’S NEW EXTRUSION PRESS PACKS 13,200-TON PUNCH 


24-inch O. D. Tubing 


This mammoth press, newest addition to The Dow Chemical 
Company’s rolling and extrusion mill at Madison, Illinois, 
is the world’s largest magnesium extrusion facility. Its vastly 
increased capacities afford new opportunities for designers 
working with light metals. A wide variety of new magne- 
sium applications for aircraft, missile, military and general 
industrial use are now possible and practical. In addition to 
extruding magnesium, the press is also available for large 
aluminum extrusions. 

Here’s how the big press will increase maximum dimen- 


sions of representative magnesium extrusions: Integrally 
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Here’s a whole new range of large magnesium extrusions: 
28-inch I-Beams 


YOU CAN DEPEND ON 


stiffened sections, from 9 to 20 inches wide; I-Beams from 
11 to 28 inches high; round tubing from 10 to 24 inches 
outside diameter; and maximum lengths of 80 feet. A large 
number of shapes and forms can be produced, limited only 
by the design of the die through which the metal is 
extruded. Many complex shapes that formerly required 
separate operations can now be formed in one operation. 


If your design calls for large magnesium or aluminum 
extrusions, contact the nearest Dow Sales Office or write to 
THE DOW CHEMICAL COMPANY, Midland, Michigan, Depart- 
ment MA 1406F-1. 


80-foot-long sections 
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iS TO KNOW...FOR SURE 
SPECIFY STATHAM 


Accelerometers « Load Cells 


Pressure Transducers 


Catalog, complete with prices, 


available upon request. 
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Throughout the Pattern of Modern Defense... 


AMF has experience you can use — 


e Most defense problems fit into AMF’s “big picture”...a picture drawn from AMF’s 
tremendous backlog of experience in defense work. Wherever you are, there’s probably 
an AMF defense component or integrated system...a product of AMF experience... 
operating near you. e Guided missiles ride on AMF equipment to AMF launching sites 
...to be loaded, fueled and fired by AMF-built equipment. Under the sea, AMF-built 
weapons wait for unfriendly submarines on the prowl. Along our borders, AMF radarscopes 
search the sky for “stranger” aircraft. e There is little room for failure where the job is the 
nation’s defense. And the nation looks to companies like AMF to design, test and produce 
a variety of defense products. With a wide range of experience in the most sensitive fields 
of defense work, AMF may well be the answer to your problems...in defense or industry. 
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letense Products Group 


AMERICAN MACHINE & FOUNDRY COMPANY 
1101 North Royal Street, Alexandria, Va. 
Asbury Park @ Atlanta @ Boston @ Brooklyn @ Dallas @ Dayton @ Los Angeles @ Seattle @ Tucson @ Washington, D.C. 
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NELECTROGRAPH 


MODEL 420 


Visual editing of fape- 
recorded intelligence 
through use of Electro- 
graph records saves 
many hours of valuable 
computer time. 


24-channel Electrophotographic 
Recording Oscillograph 


Now you can edit telemetered tape-recorded information.and select the 
data to be fed into the data reduction computer. As many as 24 channels 
of tape signals connected to the Electrograph galvanometers through ap- 
propriate discriminators may be recorded simultaneously. The analog record 
is instantly readable as it is discharged from the recorder. Proper keying 
of the permanent record to the tape permits visual selection of the data 
to be placed into the computer. 

Records produced by the Electrograph are permanent, requiring no 
further processing, and may be stored indefinitely without loss of trace 
definition. 

For further information regarding the Electrograph, you are invited to 
write, wire or call for bulletin CGC-311. 


Century Electronics & Instruments, Inc. 


1333 No. Utica, Tulsa, Oklahoma 


Book Reviews 


Ali Bulent Cambel, 
Northwestern University, 
_ Associate Editor 


Elements of Gasdynamics, by H. W. Liep- 
mann and A. Roshko, John Wiley and 
Sons, Inc., New York, 1957, xv + 439, 
$11.00 

Reviewed by ELtswortu A. Brown Jr 

Northwestern University 


The authors have presented a timel) 
book on gasdynamics which will be wel- 
come to both the student and the teacher 
It presents an up-to-date coverage of gas- 
dynamics fundamentals in one readily) 
readable volume. 

As indicated by the authors, this book 
was originally planned as a revision of th« 
well-known “Introduction to Aerodynam- 
ics of a Compressible Fluid’? by Liepmann 
and Puckett. However, one has only to 
look into the first chapter of the new vol- 
ume on ‘Concepts from Thermodynam- 
ics’”’ to realize that it is more than a re- 
vision. The vast field of gasdynamics and 
the great strides in aerodynamics have 
required considerable increase in scope. 
Whereas, in the earlier volume the discus- 
sion was generally limited to results based 

upon a perfect gas, the present volume 
includes such additional factors as real 
gas effects, condensation, dissociation and 
multicomponent systems. Many of the 
problems of present day gasdynamics re- 
quire this more general approach. 

Although most of the topics from the 
earlier volume are included in the present 
work, each topic has been rewritten to in- 
clude not only the advances of the past 
10 years, but also expanded to give a more 
complete treatment of the wider field of 
gasdynamics. For example, only optical 
methods of measurement were discussed 
in “Introduction to Aerodynamics of a 
Compressible Fluid’’ This has been ex- 
panded to give a survey of temperature, 
pressure, velocity and Mach number 
measuring techniques, as well as a discus- 
sion of special instruments such as the hot 
wire and X-ray absorption. 

Several new topics have been included. 
Chapter 10 “Similarity Rules of High- 
Speed Flow,” Chapter 11 ‘“Transonic 
Flow,” Chapter 13 “Effects of Viscosity 
and Conductivity” and Chapter 14 ‘“Con- 
cepts from Gaskinetics”’ represent material 
not included or discussed very briefly in 
the earlier volume. 

The job of covering such a large field 
as the fundamentals of gasdynamics in 
one manageable volume obviously requires 
some selectivity in the subject matter. 
One may regret the brief treatment of 
some topics such as rarefied gasdynamics 
which is limited to a few sections in the 
chapter on “Concepts from Gaskinetics.”’ 
But there is no question that this book 
contains a wealth of material in a clear 
and understandable fashion. 

The book is intended as a text at the 
senior or first year graduate level. The 
authors have done an admirable job of 
writing such a book. Without question 
the book will find wide acceptance as a 
text on gasdynamics. 
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Ratio pressure regulators 


Differential pressure 
regulators 
414 types 


CABIN AIR COMPRESSORS 


AiResearch is the largest 
designer and manufacturer of 
pneumatic controls for the air- 
craft and associated industries. 
During the past 10 years more 
than 300,000 units have been 
produced and are in service. 


Designers and manufacturers of aircraft and missile systems and components: rerriceRation systems + 
HEAT TRANSFER EQUIPMENT 


TURBINE MOTORS 
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Compartment regulators 
267 types 


GAS TURBINE ENGINES 


AIRCRAFT 


Pneumatic switches 
10 types 


Flow regulators 
28 types 


Pneumatic engine 
computer-controllers 
6 types 


Temperatures of the fluids 
(including gas and liquids) 
range from —400°F to 
+2000°F at pressures to 
+6000 psig. The units operate 
at any ambient pressure at 
ambient temperatures from 


CABIN PRESSURE CONTROLS 


* Outstanding opportunities for qualified engineers 


AiResearch Manufacturing Divisions 


Los Angeles 45, California... Phoenix, Arizona 


ELECTRO-MECHANICAL EQUIPMENT * 


—300°F to over +1000°F. 
Line diameters range from % 
inch to 15 inches. 

This equipment is developed 
and tested in the finest pneu- 
matic facilities in the world. 
Your inquiries are invited. 


PNEUMATIC VALVES AND CONTROLS + 


TEMPERATURE CONTROLS 
ELECTRONIC COMPUTERS AND CONTROLS 
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IMPREGNATED 
GLASS FABRICS 


tor applications 
60°F-10,000°F 


Send 
for latest data on Tre- 
varno impregnated glass 
fabrics for high tempera- 
ture, high strength, missile 
and aircraft applications. 
Trevarne Silicone and 
Phenolic materials have 
been tested and proved 
to give optimum heat 
resistance and strength 
performance, and exce 
tional electrical and 
mechanical properties as 
well. Also available is 
data on new Trevarno 
chopped F-120 Phenolic 
material for molded 
parts for performance 
at 4,000°F. 
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GLASS FABR 


COAST MANUFACTURING AND SUPPLY COMPANY 


& P. O. Box 71, Livermore, California § 
EB Please send me data on materials for the 
temperature ranges | have checked below: 
60°F to 160°F . oO 
800°F to 1500°F . “@ 
1500°F to 5 to 10,000°F 
New chopped Phenolic . 

Name. 
| 
Comp 
4 Address. a 
City. State. 
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New 


Gas turbine with acoustic surge control 
(2,807,932). Albert G. Bodine Jr., Van 
Nuys, Calif. 

Bladed compressor section in which a 
periodic surge flow phenomena occurs at a 
characteristic sound wave frequency. An 
acoustic wave attenuation means having a 
frequency response for the sound wave 
stall surge frequency, is acoustically 
coupled to the turbine air flow passage. 
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Magazine for releasable loads (2,807,193). 
R. A. Robert and P. P. Matge, Drancy, 
France. 

Rocket magazine for mounting in an 

aircraft fuselage. The bottom wall of the 
magazine is made of a material breakable 
upon impact of the lowermost rocket dur- 
ing its downward advance. The under- 
surface of the suspending member is 
shaped to provide a continuous surface with 
adjacent parts of the fuselage when the 
suspending member reaches its lowermost 
position. 
Propellant supply systems for jet reaction 
motors (2,808,700). Willi Kretschmer, 
Aylesbury, England, assignor to the 
Government of Great Britain. 

A normally closed main valve in the 
propellant supply system arranged to 
open in response to a rise in outlet pres- 
sure to a predetermined value less than its 
inlet. pressure. 

Injector for rocket motor (2,808,701). 
Donald W. Lewis, Los Angeles, Calif. 

Liquid fuel and oxidizer conveying tubes 
spaced from a supporting block through 
which passages supply propellant and 
oxidizer. The tubes, with multiple apera- 
tures, are disposed to define an injector 
face spaced from the block. 

Diffuser cone (2,809,491). Leon R. 
Wosika, San Diego, Calif., assignor to 
Solar Aircraft Co. 

In the front part of the duct for diffusing 

exhaust gases, the cone angle increases to 
more than 13 deg, but less than the angle 
at which the gases stall. In the center 
part, the angle decreases to less than 13 
deg, but not more than zero deg. In the 
rear part, the angle increases to 13 deg at 
its end. 
Apparatus for measuring and/or control- 
ling fuel/air ratio of gas turbines without 
direct gravimetric fuel metering (2,809,- 
492). Seymour M. Arkawy, Peekskill, 
Y., assignor to Simmonds Aerocessories, 
ne. 

Control system for maintaining a pre- 
determined gas temperature by control- 
ling the ratio of mass fuel flow to air flow 
in terms of measured pressure drop and 


George F. McLaughlin, Contributor 


variable pressure coupled to fuel/air ratio 
adjusting means. 

Guide load switch (2,809,558). Sidney 
Hersh, Silver Spring, Md., assignor to the 
U. 8. Navy. 

Guide button extending from a missile 
into a guide rail. A switch, closed when a 
missile is in loaded position, is rocked by the 
button when sliding the missile onto the 
track, latching the missile in place on the 

rail. 

Rocket projector (2,809,559). Charles C. 
Pasadena, Calif., assignor to 
the U. S. Navy. 

Channel-shaped guides for supporting 
electrically ignitible rockets. Adaptor 
guides, mounted in spaced relationship 
within the first guides, support smaller 
rockets. 

Missile firing device (2,809,624). Harold 
Becher and Martin Pepper, New York, 
assignors to Dellenbarger Machine Co. 

A switch actuated by a missile posi- 
tioned in the launching barrel controls a 
valve to permit air to escape. Gas under 
pressure ejects the missile from the barrel. 
Pressure vessel (2,809,762). Carlos J. 
Cardona, North Hollywood, Calif., as- 
signor to Fairchild Engine and Airplane 
Corp. 

A thin-walled sphere containing a liner 
of rubbery material. A reinforced exterior 
layer of glass fibers is wound on the sphere 
and impregnated with a resin. 

Connector ring for two-stage rockets 
(2,809,584). Bernard Smith, Inyorken, 
Calif., assignor to the U. 8. Navy. 

Pair of fi: anges, one fixed to the rear of 
the rocket and the other to the front of the 
booster motor. Grooved rings clamped 

around the flanges are held together by 
three explosive bolts. Actuation of the 
bolts releases the ring, disengaging the 
rocket from the booster. 


Injector for reaction motor (2,810,259). 
Harry W. Burdett Jr., Wycoff, N. J., 
assignor to Reaction Motors, Inc. 

Ports in the combustion chamber wall 
are progressively uncovered and covered 
by a piston as it slides in a slot in a direc- 
tion away from and toward the chamber. 
At least one propellant is supplied under 
pressure greater than that existing in the 
chamber. 


Epitors Note: Patents listed above were selected from the Official Gazette of the U.S. 


Patent Office. 


Printed copies of Patents may be obtained from the Commissioner of 


Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 
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Rocket propulsion method (2,810,257). 
Don R. Carmody, Evan A. Mayerle and 
Alex Zletz, Park Forest, Ill., assignors to 
Standard Oil Co. (Indiana). 

Method consists of separately and 
simultaneously injecting into the combus- 
tion chamber a nitric acid oxidizer and a 
fuel consisting of triethyltrithiophosphite 
and a liquid hydrocarbon mixture which 
boils over the range of 200 to 400 F. 
Aircraft auxiliary jet engine nacelle 
(2,810,534). Paul G. Fandeux, Cour- 
bevoie, France, assignor to SNCASO. 

Streamlined nacelle fixed under the 
fuselage of an airplane. The leading edge 
of the intake mouth is located at a dis- 
tance from the nose tip at least equal to 
half the diameter of the cross section of the 
nacelle. 


(181,317). 
Gasich, 


Edgar 
Los An- 
geles, Calif., assignors to Northrop Air- 
craft Corp. 

Jet propelled airplane with two engine 
pods suspended below the delta wing. 
Method for resizing internally threaded 


Aircraft design 
Schmued and Welko E. 


rocket motor tubing (2,811,730). Conrad 
R. Vegren, Washington, D. C. 

Means for reducing the diameter of the 
end of a tube. Force is applied radially 
inward to the end portion of sufficient 
magnitude to permanently deform it to 
smaller diameter. 

Balanced jet nozzle and cooling means 
therefore (2,811,827). Ralph Kress, La 
Mesa, Calif., assignor to Solar Aircraft Co. 

An actuating mechanism imparts sim- 

ultaneous opposite movements to a pair of 
coaxial members defining an exhaust ori- 
fice. Members can be rapidly moved 
with minimum effort to vary the area of 
the exhaust orifice for various operating 
conditions. 
Bleed means for confinea supersonic flow 
(2,811,828). George H. McLafferty, Man- 
chester, Conn., assignor to United Air- 
craft Corp. 

Ramjet with bleed means composed of 

apertures in the inlet diffuser wall. Ap- 
pertures comprise upstream and down- 
stream lips along the axis of flow. Down- 
stream lips are recessed, forming a re- 
entrant angle with the wall. 
Ramjet employing carbon layer of insula- 
tion for solid carbon propellant (2,811, 
829). Alfred A. Topinka, Cleveland, 
Ohio. 

Combustion chamber lined with solid 
carbon as a fuel and propellant. Carbon 
is protected from the shell by insulation of 
lamp black and powdered carbon. A 
complimentary burning section has an 
additional air inlet and a lining of refrac- 
tory material. 


DECEMBER 1957 


Men who know 


TURBINES, PUMPS, COMPRESSORS 


Get in on the development of the most 
advanced high-speed rotary equipment ever built 


Here’s your opportunity to step up to a new, higher level of 
turbo-machinery technology—the Large Rocket Engine. 

Whatever type of rotating machinery you know best, your 
experience could be extremely valuable in the important develop- 
ments now going on at Rocketdyne. The seasoned and ambitious 
man who has cut his teeth on jet engines, steam or gas turbines, 
or other elements of rotating machinery, is urgently needed to 
apply mechanical principles to meet the increasing demands of 
power plant performance. 

The combination of high speed, light weight, heavy loadings 
and exceptional pressures required in rocket engine work is lead- 
ing to an entirely new breed of high-performance rotating 
machinery ...and a new breed of engineer. You can be one of 
this advance guard of the turbo-machinery field—if you have the 
desire to build your professional status by accepting new 
challenges. 

Rocketdyne is building high-thrust rocket engines for the 
nation’s major missiles. You'll work with the leading producer in 
the nation’s fastest growing industry. You and your work will 
be recognized as a vital part of the overall achievements. Testing 
facilities are among the world’s finest. The power produced is 
beyond anything ever before thought possible. If you would like 
to tackle new assignments working alongside some of the finest 
minds in turbo-engineering, write and tell us about your back- 
ground: A. W. Jamieson, Rocketdyne Engineering Personnel 
Dept. R-12, 6633 Canoga Avenue, Canoga Park, California. 


ROCKETDYNE 


ON OF NORTH AMERICAN AVIAT 


BUILDERS OF POWER FOR OUTER SPACE 
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Now... 
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X=ray of the Falcon’s fin proveshe 


Now you can have structurals die cast to standards as high qf°re- ¢ 
those for primary aircraft structurals*. below s 
Die cast! With all the economy, precision, surface superiori 
ease of assembly that implies. Die cast of aluminum, or of Like 
high-purity magnesium alloy, like the 42” Falcon stabilizer fi e = bs 
a custor 
They said it couldn't be done! They, 


But there you see such a structural! Along with X-ray pro parts ™ 
of soundness, uniformity, primary-aircraft-structural* in w 

Mind you! The Falcon fin is in commercial productio§’°’° 2” 
Doehler-Jarvis has die cast thousands. Each has had to ps It too 
73 inspections before acceptance ... including 100% X-ray. ff both. 


The Falcon is an air-to-air missile built by Hughes Air- ; . : N 
craft Company for the U. S. Air Force. Its light weight In one control, sections of castings are continuously test ow } 
die cast stabilizer fins resist severe strain during course for tensile strength. In this examination of physical propertit ay 

changes. It uses other Doehler-Jarvis die castings also. the castings normally exceed_tensile specifications by 15% oehler- 


: 


heir uniformity and soundness 


more. Quality is so high there has never been a single test *A primary aircraft structural is a major load-bearing component member 


. * any single failure of which would result in loss of the aircraft. 
below specification. 


How did this come about? 


Like many developments in die casting, this one started with Doehler=Jarvis 


a customer’s need... Hughes Aircraft’s. 

They, and the Armed Services, needed such castings. The DIVISION OF NATIONAL LEAD COMPANY 
parts might be made other ways. But at what cost! What sacri- General Offices: Toledo 1, Ohio 
ice in weight — and materials! If it could be die cast, if there 
were a magnesium alloy able to stand supersonic air stresses... In Canada: 

It took a while. But, with close cooperation from engineers Barber Die Casting Co. Limited 
of both Doehler-Jarvis and Hughes, the way was found. Hamilton, Ontario 

Now you can list die castings of primary-aircraft-structural* 
quality among D and J availabilities. Get them from any 
Voehler-Jarvis plant. 
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BINARY 


BINARY 
CODED DECIMAL 


LIBRASCOPE 
SHAFT POSITION-TO-DIGITAL 


CONVERTERS 


Equipped with ANTI-AMBIGUITY 
DOUBLE BRUSH PICKOFFS 


Useful in a wide variety of applications, 
including digital aircraft and missile 
controls, machine tool controls, digital 
readout from strip chart recorders, and 
as the modulator and de-modulator in 
pulse-code modulated radio links. 


GRAY CODE MODEL —Capacity of 8 
binary digits (single brush pickoff), 
BINARY MODEL —Capacity of 7 to 19 
binary digits. 

BINARY CODED DECIMAL MODEL 

- Capacity range from 0-1999 to 0-35,999. 


Units for special codes or 
capacities are built to meet 
specific requirements. 


SHOCK ENDURANCE................... 20g 
TEMPERATURE RANGE. .—50° to 83°C min. 
CODE DISCS....... Rhodium plated phenolic 
Sere Multiple wire brush. 
Two pickoffs /channel 

OTATION...... Continuous, either direction. 


RUGGED—NON-MAGNETIC—LONG LIFE 
MAY BE READ WHILE IN MOTION 


SPECIAL CONVERTERS DESIGNED TO MEET 
YOUR INDIVIDUAL PROBLEMS 


Send for illustrated brochure 
1IBRASCOPE 


A SUBSIDIARY OF GENERAL PRECISION EQUIPMENT CORPORATION 


808 Western Avenue « Glendale, California 
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| Cross Sections, by 


Technical Literature Digest___ 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 


Nuclear Propulsion 


The Lithium Hydride, Deuteride and 
Tritide Systems, by F. K. Heumann and 
O. N. Salmon, Atomic Energy Comm., 
KAPL-1667, Dec. 1956, 54 pp. 

Fast Fuel Rod Design and Analysis, by 
C. R. Stahl, R. J. Futz and R. G. Kenni- 
son, Atomic Energy Comm., KAPL-423, 
Nov. 1950, 47 pp. 

Neutron Energy Distributions Inside 
the Fast Reactor, by Norris Nereson, 
Atomic Energy Comm., LA-1192, Dec. 
1950, 49 pp. 

Pressure Drop Tests on Twisted Ribbon 
Core Assemblies, by A. B. Schultz, Atomic 


| Energy Comm., ANL-5189, Aug. 1953, 42 


pp. 

Reactions Produced by the Electrical 
Explosion of a Metal Immersed in a 
Fluid, by Robert F. Plott, Atomic Energy 
Comm., ANL-5040, Dec. 1950, 35 pp. 

A Correlation of Fast Neutron Fission 
John R. Huizanga, 
Atomic Energy Comm., ANL-5150, Nov. 
1953, 7 pp. 

Studies of Mixing in the Lower Plenum 


Ep:Tor’s Note: Contributions from E. R. 
G Eckert, J. P Hartnett, T. F. Irvine Jr. and 
P. J. Schneider of the Heat Transfer Labora- 
tory, University of Minnesota. are gratefully 
acknowledged 


PLAIN TYPES 


ANALYSIS 


' Stainless Steel 
Ball and Race 


2 Chrome Alloy 
Steel Ball and Race 


Chrome Stee! Ball 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


The James Forrestal Research Center, Princeton University 


Self-Aligning Bearings 


ExT. 
INT. 
PATENTED U.S.A. 
World Rights Reserved 


TERISTICS 


RECOMMENDED USE 


For types operating under high temper- 
ature (800-1200 degrees F.). 


For types operating under radial 
ultimate loads (3000-893,000 


3 Bronze Race and normal loads 
with minimum friction requirements. 

Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. 


SOUTHWEST PRODUCTS CO. 


of a Quarter Scale Flow Model of the 
PWR Reactor, by Herbert R. Hazard and 
Abraham Rotkowitz, Alomic Energy 
Comm., BMI-1172, Feb. 1957, 19 pp. 

Study of a Liquid Metal Neutron Ab- 
sorber for the Control cf a Gas Cooled 
Power Reactor, by Banjamin B. Gordon, 
J. James Stone and Harold M. Epstein, 
Atomic Energy Comm., BMI-1119, July 
1956, 18 pp. 

Enrico Fermi Fast Breeder Reactor 
Plant (Atomic Power Development Asso- 
ciates, Inc.,) Atomic Energy Comm., 
APDA-115, Nov. 1956, 107 pp. 

Nuclear Powered Gas Turbines for 
Light Weight Power Plants, by Frederick 
G. Hammitt and Harold A. Ohlgren, 
ASME Pap. 57-NESC-79, March 1957, 26 

pp. (Nuclear Engineering and Sci. Conf, 
2nd., Philadelphia, March 11-14, 1957). 


Fuel Cycles in Single Region Thermal 
Reactors, by Manson Benedict and 
Thomas H. Pigford, ASME Pap. 57- 
NESC-41, March 1957, 21 pp. (Nuclear 
Engng. and Sci. Conf., 2nd., Philadelphia, 
March 11-14, 1957). 

The Diffusion Theory of Fine Structure 
in Thermal Neutron Reactor Assemblies 
Consisting of Cylindrical Fuel Elements 
Set in a Square Lattice Array, by A. C. 
Clark and D. A. Newmarch, Gt. Brit., 
Atomic Energy Res. Estab., Rep. RP/R 
1657, 1956, 18 pp. 


Theory of the Ideal A.C. Conduction 


“MONOBALL”’ 


Dept. JP-57. 
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Is YOUR JET ENGINE 


RPM MEASUREMENT 


, READS JET ENGINE SPEED to 
GUARANTEED ACCURACY of 
A 10 RPM in 10,000 RPM (=+0.1%) 


The inter-relation of RPM to efficiency and thrust in jet engines is 
fundamental. Proper adjustments for maximum thrust, maximum engine 
life and maximum safety of operation can be made only upon accuracy 
of instrumentation. The TAKCAL tests to guarantee that accuracy. 

The new B & H TAKCAL incorporates a refinement of the frequency 
meter principle. It operates in the low (0 to 1000 cps) range, reading 
the frequency of the tachometer generator on a scale calibrated in 
percent RPM corresponding to the engine’s RPM. It reads engine speed 
while the engine is running with a guaranteed accuracy of +0.1% in 
the range of 0 to 110% RPM. Additionally, the TAKCAL circuit can be 
used to trouble-shoot and isolate errors in the aircraft tachometer 
system, with the circuit and tachometer paralleled to obtain simultaneous 
reading. 


CHECKS TACHOMETER 
“SYSTEM” ACCURACY. 
The TAKCAL’s component parts are identical with those used in the ADAPTS TO ANY FREQUENCY PROBLEM! 
J-Model JeTcat Analyzer. They are here assembled as a separate unit plosion-proof TAKCAL 
tester and for use with all earlier models of the JETCAL Tester. for special applications. 
The TAKCAL operates accurately in all ambient temperatures from Measures 200 to 7500 
—40°F. to 140°F. Low in cost for an instrument of such extreme RPM, direct reading, with 
accuracy, it is adaptable to application in many other fields. ae 
For full information write or wire... 


B & H makes the JETCAL® Analyzer 


B & H INSTRUMENT Co., INC. dye and TEMPCAL® Tester 
<5 3479 West Vickery Blvd. * Fort Worth 7, Texas 


. Sales-Engineering Offices: 


VALLEY STREAM, L. 1: 108 So. Franklin, LO 1-9220 © DAYTON, OHIO: 209 Commercial Bldg., BA 4-4563 © COMPTON, CAL.: 105 N. Bradfield St., NE 6-8970 
ENGLAND: Bryans Aeroquipment Ltd. (Licensee), 15, Willow Lane, Mitcham, Surrey 
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ne Product 


REVERE SPECIALTY: WIRES 


ABRASION 


4) 


For those tough design jobs where ordinary hook-up and 
thermocouple wires die from the heat, get brittle in cold, 
abrade and corrode .. . Revere SPECIALTY wires stand up. 
Built to MIL and customer specifications. Range includes: 


For High Temperatures 


REVCOTHENE — (Extruded Monochlorotrifluoroethylene) 
—40°F to +275°F, AWG 28 to 10, silver-plated copper 
conductors, inert, excellent dielectric strength, no volatile 
plasticizers, non-flammable, thin wall, abrasion and moisture 
resistant. 


HOOK-UP 
WIRES  PERMACODE -— Teflon*— insulated wire with striping down 


- ¥ be. to the conductor for permanent identification, single or mul- 
Te tiple stripes, 15 colors, —130°F to +410°F, AWG 28 to 16, 
silver-plated copper conductors, excellent abrasion and di- 

electric characteristics. 
= MULTI- & A variety of telemetering and other multi-conductor cables 
CONDUCTOR I constructed to customer specifications. Teflon, polyethylene, 
CABLES @ Polyvinyl, nylon, glass, Revcothene, asbestos insulations for 
singles and jackets. Twisting, braiding, shielding, color coding 

to suit conditions. 


THERMO- lron-constantan, copper-constantan, Chromel-Alumel conduc- 
COUPLE tors, AWG 36 to 14, various insulation combinations and pro- 
wires @ ‘tive braids, temperature range from —100°F to +1500°F, 
constructed to rigid tolerances. 


Saturants for flame and abrasion resistance, metallic braids for severe service 
and electrical shielding. Color coding in 15 solid colors and stripes. 


Prompt delivery of standard stock wires. Write for samples 


i and literature on specialty hook-up or thermocouple wire. 
20 


REVERE CORPORATION OF AMERICA 
Wallingford, Connecticut 
A SUBSIDIARY OF NEPTUNE METER COMPANY 


Pump, by W. Murgatroyd, Gt. Brit., 
Atomic Energy Res. Estab., Rep. ED/R 
1566, 1956, 9 pp., 7 fig. 

A oe Generator for High Cur- 
rent Low Voltage D.C. Supply, by D.A. 
Watt, Gt. Brit., Atomic Energy Res. 
Estab., Rep. ED/R 1843. (Declassified 
edition of AERE Rep. CE/R821), 1956, 
24 pp. 

The Time and Temperature Depend- 
ence of Thermal Stresses in Cylindrical 
Reactor Fuel Elements, by Kenneth R. 
Merckx, ASME, Pap. 57-NESC-78, 
March 1957, 8 pp. (Nuclear Engng. and. 
Sci. Conf., 2nd., Philadelphia, March 11 
14, 1957). 

Decay and Storage of Irradiated Fuel, 
by J. W. Ullman and E. D. Arnold, ASME 
Pap. 57-NESC-93, March 1957, 6 pp. 
(Nuclear Engng. and Sci. Conf., ond. 
Philadelphia, March 11-14, 1957). 

UO.-NaK Slurry Studies in Loops at 
600°C., by B. M. Abraham, H. E. Flow- 
tow and R. D. Carlson, ASME Pap. 57- 
NESC-104, March 1957, 12 pp. (Nuclear 
Engng. and Sci. Conf., 2nd., Philadelphia, 
March 11-14, 1957). 

Transient Thermodynamics of Reactors 
and Process Apparatus, by Dale H. Brown, 
ASME, Pap. 57-NESC-81, March 1957, 
11 pp. (Nuclear Engng. and Sci. Conf., 2nd., 
Philadelphia, March 11-14, 1957). 

Open Cycle Takes Nuclear Engine Lead, 
by Robert Cushman, Aviation Week, vol. 
66, May 27, 1957, p. 27. 

A Possible Fusion Reactor, by Goro 
Miyamoto, Gliti Iwata, Sigeru Mori and 
Kenji Inoue, J. Phys. Soc. Japan, vol. 12. 
April 1957, p. 438. 


Combustion, Fuels and 
Propellants 


Interaction of a Free Flame Front with a 
Turbulence Field, by Maurice Tucker, 
NACA Rep. 1277, 1956, 19 pp. (super- 
sedes NACA TN 3407). 

Effect of Standing Transverse Acoustic 
Oscillations on Fuel-Oxidant Mixing in 
Cylindrical Combustion Chambers, by 
William R. Mickelsen, NACA TN 3983, 
May 1957, 49 pp. 

Transverse Wave and Entropy Wave 
Combustion Instability in Liquid Propel- 
lant Rockets, by Sinclaire M. Scala, 
Princeton Univ., Dept. Aeron. Engng., Rep. 
380, 1957, 138 pp. (PhD. thesis). 

Longitudinal Combustion Instability in a 
Duct with Oblique Flame Front, by Marcel 
Vinokur, Princeton Univ., Dept. Aeron. 
Engng., May 1957, 118 pp. (PhD. thesis). 

Generalized Design Equations for an 
Internal Burning Star Configuration Solid 
propellant Charge and Method of Cal- 
culating Pressure Time and Thrust Time 
Relationships, by L. Piasecki and G. 
Robillard, Calif. Inst. Tech., Jet Propulsion 
Lab., Mem. 20-135, Sept. 1956, 33 pp. 

Combustion in Laminar Mixing Regions 
and Boundary Layers, by Donald Allen 
Dooley, Pasadena, Calif. Inst. Tech., 1956, 
244 pp. (Ph.D. thesis). 

Rate of Recombination of Radicals, 
III: Rate of Recombination of Ethyl 
Radicals, by A. Shepp and K. O. Kutschke, 
J. Chem. Physics., vol. 26, May 1957, pp. 
1020-1028. 

A Precision Calorimeter for the Meas- 
urement of Heats of Combustion, by 
Charles E. Holley Jr. and Elmer H. Huber 
Jr., Atomic Energy Commission, LA-2084, 
Jan. 1957, 18 pp. 

Temperature Measurements During 
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Firing of Simulated Submarine Rocket 
Launcher, by Norman E. Jensen Jr. and 
James C. Coe, Naval Air Missile Test 
Center, Tech. Rep. 16, June 1947, 17 pp. 
(Declassified from Confidential by author- 
ity of NAMTC letter A9/Ser. 362, 7/20/ | 
56.) 

Reactions of Gaseous Ions, II: Acety- | 
lene, by F. H. Field, J. L. Franklin and 
Ff. W. Lampe, J. American Chem. Soc., 
vol. 79, June 5, 1957, p. 2665. 

Some Factors Influencing Stability 
Limits of Bunsen Flames, by P. F. Kurz, 
Combustion and Flame, vol. 1, no. 2, June 
1957, pp. 162-178. 

Mixing and Combustion in Free and En- 
closed Turbulent Jet Diffusion Flames, 
by P. D. Sunavala, C. Hulse and M. W. 
Thring, Combustion and Flame, vol. 1, no. 

2, June 1957, pp. 179-193. 

A Law of Similarity for Temperatures in 
Enclosed Turbulent Flames, by C. Code- 
gone, Combustion and Flame, vol. 1, no. 2, 
June 1957, pp. 194-198. 

On the Positive and Negative Ions in 
Diffusion Flames, by Tosiro Kinbara and 
Hideo Ikegami, Combustion and Flame, 
vol. 1, no. 2, June 1957, pp. 199-211. 

Influence of Diborane on Blow-Off 
Limits of Hydrocarbon Flames, by P. F. 
Kurz, Combustion and Flame, vol. 1, no. 2, 
June 1957, pp. 212-216. 

Flame Zone Studies IV—Microstruc- 
ture and Material Transport in a Laminar 
Propane-Air Flame Front, by R. M. Fris- 
tom and A. A. Westenberg, Combustion 
and Flame, vol. 1, no. 2, June 1957, pp. 
217-228. 

Experimental Techniques for the Study 
of Flat Flames by Infra-red Spectroscopy, 
by D. J. Cole and G. J. Minkoff, Combus- 
tion and Flame, vol. 1, no. 2, June 1957, pp. 
241-249. 

Homogeneous Decomposition of Hydro- 
gen Peroxide Vapor, by Charles N. Satter- 
field and Theodore W. Stein, J. Physical 
Chem., vol. 61, May 1957, pp. 537-540. 

Calculation of Experimental Isotope 
Effects for Pseudo-First Order Irreversible 
Reactions, by James Yang-peh Tong and | 
Peter E. Yankwich, J. Physical Chem., vol. | 
61, May 1957, pp. 540-543. } 

Studies of Diffusion Flames, III: The | 
Diffusion Flames of the Butanols, by S. | 
Ruven Smith, Alvin 8. Gordon and May- 
nard H. Hunt, J. Physical Chem., vol. 
61, May 1957, pp. 553-558. 

A Kinetic Approach of the Thermody- 
namics of Irreversible Processes, by O. K. 
Rice, J. Physical Chem., vol. 61, May 1957, 
pp. 622-629. 

Nitric Oxide Decomposition at 2200- 
2400° K., by C. P. Fenimore and G. W. 
Jones, J. Physical Chem., vol. 61, May 
1957, pp. 654-657. 

Thermal Design of Large Storage Ves- 
sels for Liquid Hydrogen and Helium, by 
Russell B. Scott, J. Res., Nat. Bur. Stand- 
ards., vol. 58, June 1957, pp. 317-326. 

Chemical Aspects of High Energy Solid 


Watch CBS 


Propellants, Missiles and Rockets, vol. 2, 
June 1957, pp. 110, 112. : 
Field Methods for Determination of a a 


Unsymmetrical Dimethylhydrazine and 
Water in NIKE Fuels, by R. F. Muraka, 


for the Du Pont 


Why Du Pont Aircraft 
Rivets are used for 


high-temperature 


fastening 
in the J79 


While developing their J79 gas turbine 
engine, General Electric had to find a depend- 
able blind fastener to withstand elevated 
temperatures. 

Du Pont High-Temperature A-286 Super- 
alloy stainless steel rivets solved the problem 
—fastening the diaphragm section inside the 
engine housing. 

From a production engineering standpoint, 
here are the important reasons for General 
Electric’s choice of Du Pont Rivets. 
Temperature Range—Du Pont A-286 Rivets 
will retain their high-strength properties up to 
1400° F. 

Fast setting—Working from the head end of _ 
the rivet, one man can easily set 20-25 rivets — 
per minute. 

No Refinishing—There’s no need for cutting, 
trimming, or buffing after the rivets are set. 
Positive holding power—The rivet shank ex- _ 
pands to fill the hole completely, insuring — 
tight joints. 

Find out how Du Pont Rivets can solve — 
your high-temperature fastening problems. — 
Call your DuPont representative or write: 
E. I. du Pont de Nemours & Co. (Inc.) , 2543 
Nemours Building, Wilmington 98, Delaware. — 


Stephen P. Bango, Lois L. Taylor and 
Julia S.: Whitlick, Calif. Inst. Tech., Jet 
Propulsion Lab., Prog. Rep. 20-315, March 
1957, 11 pp. 

Propellant Stability Testing at Picatinny 
Arsenal; a Résume of Facilities, Methods 


DUPONT AIRCRAFT BLIND EXPANSION RIVETS 


A Product of Du Pont Research 


and Programs, Picatinny Arsenal, Samuel ——— 


Feltman Ammunition Labs., TN 8, June 
1957, 28 pp. 


BETTER THINGS FOR BETTER LIVING... THROUGH CHEMISTRY 
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This Visicorder Oscillograph record* is a sym- 
bol of the leadership that is typical of Honeywell 
engineering. In laboratories all over the world 
the Visicorder’s instantly-readable direct records 
are showing the way to new advances in rocketry, 
control, computing, product design and com- 
ponent test and in nuclear research. 


*reproduced actual size, unretouched 


The Model 906 Visicorder is years ahead of the 
trend. It is the first oscillograph that combines 
the convenience of direct recording with the high 
frequencies and sensitivities of photographic- 
type instruments. The Visicorder alone among 
oscillographs lets you monitor high-speed vari- 
ables as they go on the record. 


record of leadership | ; 
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Some of the general features which give the Visi- 
corder leadership in the direct-recording field are: 


@ Frequencies from DC to 2000 cps without 
peaked amplifiers or other compensation 


@ Six channels plus 2 timing traces on 6” paper 
@ Deflection 6” peak to peak; traces may overlap 


@ Record speeds 0.2, 1, 5, and 25 inches per 
second, minute, or hour 


@ Records require no liquids, powders, vapors, 
or other processing 


Call your nearest Minneapolis- Honeywell Industrial 
Sales Office for a demonstration. 


Reference Data: Write for Visicorder Bulletin 


Minneapolis-Honeywell Regulator Co., Heiland Division, 
5200 East Evans Avenue, Denver 22, Colorado. 
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YIELD STRENGTH/DENSITY/10° 


OFFICES IN NEW YORK, 


o 


7075-16 


302 STAINLESS 


20 400 600 800 1000 
TEMPERATURE ‘F 
‘* Yield strength, weight vs. temperature. 


For High 
Strength-to-Weight 
Ratios at Elevated 
Temperatures — 


TITANIUM and 
MAGNESIUM 


The experience, the specialists, the 
equipment-and-facilities to design 
and produce the very difficult fabri- 
cation and assembly work in these 
and other light metals—you'll find 
at B&P. 


To make this titanium 
pressurized tank, B&P 
welds two cylinders. 
Each is drawn in one 
operation. Tank with- 
stands 5000 psi pressure. 


Magnesium turret 
enclosure for a 
bomber. Design, 
lofting, prototypes, 

production fabrica- 
tion and assembly 
are all done by 
B&P 


ENGINEERING FACILITIES IN DETROIT 
AND LOS ANGELES .. Write for Titanium 
and Magnesium engineering data—also 
folder on B&P’s facilities to handle 
tough jobs. 


BROOKS & PERKINS, Inc. 


1932 West Fort St. 
Detroit 16, Mich. 
Phone: TAshmoo 5-5900 


IN LOS ANGELES: 
11655 Vanowen St. 
North Hollywood, Cal. 
Phone STonley 7-9665 


WASHINGTON, DALLAS 


as Ideal Monatomic 
Gases, by Harwood G. Kolsky, Atomic 
Energy Commission, LA-2110, March 
1957, 138 pp. 

Basic Chemistry of High Temperature 
Inorganic Systems, Semi-Annual Prog- 
ress Report, Jan.-June 1956, by S. J. 
Yosim and T. A. Milne, Atomic Energy 
Commission, NAA-SR-1797, March 1957, 
23 pp. (North American Aviation, Inc.). 


Heat Transfer and Fluid 
Flow 


Nusselt Values for Estimating Turbu- 
lent Liquid Metal Heat Transfer in Non- 
Circular Ducts, by James P. Hartnett and 
Thomas F. Irvine Jr., ASME, Pap. 57- 
NESC-30, March 1957, 8 pp. (Nuclear 
Engng. and Science Conference, 2nd, 
Philadelphia, March 11-14, 1957). 

Gas Coolant for Nuclear Reactors, by 
Michael Silverberg, ASME, Pap. 57- 
NESC-76, March 1957, 4 pp. (Nuclear 
Engng. and Science Conference, 2nd., Phila- 
delphia, March 11-14, 1957). 

A Method of Cascade Design for Two- 
Dimensional Incompressible Flow, by S. 
Rosenblat and L. C. Woods, Australia, 
Aeron. Res. Comm., Rep. ACA 58, March 
1956, 30 pp. 

Bibliography on Shock Waves, Shock 
Tubes and Allied Topics, by K. Dolder 
and R. Hide, Gt. Brit., Atomic Energy 
Res. Estab. Rep. G/R 2055, 1956, 97 pp. 

Bibliography on the Origin of Turbu- 
lence and Related Subjects, by Maurice 
H. Smith, Princeton Univ., James Forrestal 
Res. Center, Lit. Search 13, April 1957, 126 
pp. 

Experiments on the Inhibition of Ther- 
mal Convection by a Magnetic Field, by 
Y. Nakagawa, Proc. Royal Soc., vol. 
A240, April 24, 1957, pp. 108-113. 

The Aerodynamic Variable Nozzle, by 
A. I. Martin, J. Aeron. Sci., vol. 24, May 
1957, pp. 357-362. 

Fluctuations and Hot Wire Anemometry 
in Compressible Flows, by Mark V. 
Morkovin, North Atlantic Treaty Organiza- 
tion, Advisory Group Aeron. Res. Dev.. 
AGARDograph 24, Nov. 1956, 102 pp. 

Shock Waves in a Relaxing Medium, by 
S. P. D’yakov, Russian Transl. 3908, 
(translation of Zhur. Eksp. Teor. Fiz., vol. 
27, no. 6, 1954, pp. 728-734). 

Shock Waves in Binary Mixtures, by 
S. P. D’yakov, Russian Transl. 3922 
(translation of Zhur. Eksp. Teor. Fiz., vol. 
27, no. 3, 1957, pp. 283-287). 

The Stability of Shock Waves, Investi- 
gation of the Problem of the Stability of 
Shock Waves in Arbitrary Media, by S. P. 
D’yakov, Russian Transl. 4295, (transla- 
tion of Zhur. Eksp. Teor. Fiz., vol. 27, no. 
3, 1954, pp. 288-295). 

Two-Dimensional Flow in the Vicinity 
of the Stagnation Point of an Incompressi- 
ble Viscous, Electrically Conducting Field 
in the Presence of a Magnetic Field, by 
Joseph L. Neuringer, Republic Aviation, 
Inc., June 1957, 19 pp. 

Mass Transfer Cooling in High Speed 
Laminar Couette Flow, by E. R. G. 
Eckert and P. J. Schneider, Minnesota 
Univ., Mech. Engng. Dept. TR 12 (AF 
Off. Sci. Res. TN 57-186; ASTIA AD- 
126,481), April 1957, 44 pp. 

Static Pressure Distribution in the Free 
Turbulent Jet, by D. R. Miller and E. W. 
Comings, Indust. and Engng. Chem., vol. 
49, June 1957, p. 985. 

The Ranque-Hilsch Vortex Tube, by 
William A. Scheller and George M. Brown, 


Manufacturing 
at Marquardt 


by 


Roy E. Marquardt 
President 


Another barrier — the PRODUCIBILITY 
BARRIER — is currently being penetrated 
by Marquardt engineers. What do we 
mean by PRODUCIBILITY BARRIER? 


Advanced designs for supersonic ram- 
jet powerplants coming from the draw- 
ing boards call for strength to weight 
ratios and high-precision tolerances 
previously unobtainable. New high- 
temperature alloys are meeting the 
metallurgical demands, but do not 
readily lend themselves to conventiona! 
machining and fabricating techniques. 

Those members of the Marquardt 
team charged with pioneering new 
production methods comprise our Van 
Nuys Manufacturing Division. Here— 
under the direction of John S. Liefeld 
— creativeness and imagination join 
forces with a thorough understanding 
of standard shop practices to produce 
acceptable hardware. But management 
realizes that to do his best work the 
engineer must be supplied with the 
most up-to-date tools of his trade. 

Exemplifying our continuing efforts 
to this end: a specially designed, half- 
million dollar roll-former is being 
added to the company’s ever-expanding 
production facilities. This machine will 
be capable of spinning conical, tubular, 
venturi, and parabolic configurations 
of a size heretofore considered impos- 
sible or impractical to fabricate as a 
single piece. A completely safe, close- 
up view of the actual metal forming 
will be afforded the operator by means 
of two closed-circuit TV cameras 
mounted on the machine frame. 

We are also acquiring other auto- 
matic machines — numerically con- 
trolled units capable of multiplying the 
output of their manually operated 
counterparts several times. Utilizing 
punched and magnetic tape, these 
machines are expected to greatly 
expand the scope of Manufacturing 
Engineering. 

Another of the modern production 
techniques at Marquardt puts to use 
modern optics equipment for the con- 
struction of prototype engines. 

At Marquardt, the engineer will find 
a broad range of challenging assign- 
ments and the opportunity to further 
his career through supplemental educa- 
tional programs. 

Within this Division, engineering 
openings exist now for: 

MErTHOps ENGINEERS Process ENGINEERS 
Toot ENGINEERS AUTOMATION ENGINEERS 
NuMERICAL CONTROL ENGINEERS 

For information about these posi- 
tions and the professional engineering 
environment at Marquardt, we invite 
you to write Jim Dale, Professional 
Personnel, today. 


marquar \AIRCRAFT CO. 
VAN NUYS, CALIFORNIA OGOEN, UTAH 
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Shown here: John S. Liefeld, Director Van Nuys Manufacturing Division 


Marquardt Means Opportunity- Manufacturing engineers no longer need 
feel stymied because of the lack of up-to-date equipment. At Marquardt Aircraft — the com- 
pany where an ENGINEER/BARRIER* has never existed — you will find the most advanced tools 
of your trade. Look to your future by looking to Marquardt, today. Address your inquiries to 
Jim Dale, Professional Personnel, 16556 Saticoy Street, Van Nuys, California. 
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LAVEZZ| MACHINE WORKS 


4635 WEST LAKE ST. 


PRECISE 
CURRENT 
for your 
TESTING 
NEEDS! 


VARIABLE FREQUENCY 
MOTOR GENERATOR SETS 
ADJUST 360 TO 440 CPS. 


Generator mounted controls 
include reset buttons, limit 
switch Motor and generator 
remain stationary. Vari-drive 
ulley adjustment controlled 
small motor. Remote con- 
trol panels available 


Units can be equipped with 
synchronous motor starter and 
magnetic amplifier, auto- 
matic voltage regulator. 


KATO 400 CYCLE 
MOTOR GENERATOR SETS 
NOW UP TO 250 KWI 


KATO MOTOR GENERATOR 
SETS are available in frequencies, es 


speeds and sizes for every special- INPUT 

ized use . . . operating high cycle 60 CYCLES 
tools, testing components end OUTPUT 
electronic equipment 400 CYCLES 


WRITE FOR NEW FOLDER! 


Fin Fecal Machinery Since 1928 


1498 First Ave., Mankato, Minn. 


¥ “=< 


CHICAGO 44, ILL. 


NORTHAM 


DIFFERENTIAL 


PRESSURE 
TRANSDUCER 


A variable inductance instrument for 


measurement of differential pres- 
sures in ranges + 200 to + 5,000 psi 
and line pressures up to 5,000 psi. 
Designed for rocket and jet engine 
testing involving extreme mechani- 
cal disturbances, and for general 
laboratory measurement and indus- 
trial control applications. 

MODEL DP-15 SPECIFICATIONS: 
Pressure Ranges:..+ 200 to + 5,000 psid 
Maximum Line Pressure:......... 5,000 psi 
Excitation Frequency: From 60 to 20,000 cps 
Natural Frequency: From 2,000 cps to 

10,000 cps depending on range 
Acceleration Response: Less than 1%/1,000 
g along most sensitive axis 
(range above 1,000 psi) 
less than 1%/1,000 g axial 
(all ranges) 


WRITE FOR BULLETIN NO JP-8 
Northam Engineering Facilities Are Available 
To Assist You With Any Application Problems. 
NORTHAM PRODUCTS INCLUDE... 


Transducers for pressure, acceleration and 
displacement measurement and auxiliary 
electronic equipment for complete systems. 


NORTHAM ELECTRONICS, INC. 


A Subsidiary of Norris-Thermador Corp. 
2420 North Lake Avenue, Altadena, Calif. 


Indust. Engng. Chem., vol. 49, June 1957, 
pp. 1013-1016. 


Laminar Boundary Layer Flows with | 


Surface Reactions, by Andreas Acrivos 


and Paul L. Chambré, Indust. Engnq. 
Chem., vol. 49, June 1957, pp. 1025-1029, 


Instrumentation, 
Telemetering, Data 
Recording 


A Nonlinear Instrument Diaphragm, by — 


l‘idel Corero, Harry Matheson and Daniel 


P. Jehnson, J. Res. Nat. Bur. Stands., vol, ~ 


58, June 1957, pp. 333-338. 

Redstone In-Flight Telemetering, by 
Charles T. N. Paludan, Missiles and Roci:- 
els, vol. 2, June 1957, pp. 126, 128. 

The Design of an Acceleration Insensi- 
tive Skin Friction Balance for Use in 
Free Flight Vehicles at Supersonic Speeds, 
by Willis Carson Lyons Jr., Teras Unir., 
Defense Res. Lab., CF-3605 (DRL-397), 
June 1957, 42 pp., 18 fig. 

Ultrasonic Flowmeter Monitors Reactor 
Heat Exchange Circulation, by M. Kob- 
lenz, R. Spiegel, S. Kass and H. Sterling, 
ASME Pap. 57-NESC-84; Amer. Inst 
Elect. Engrs., Pap. 57-400, March 1957, 
8 pp. 

A Resistance Temperature Detector for 
Nuclear Service, by F. Ralph Sias ASME 
Pap. 57-NESC-83; March 1957, 6 pp. 

Recent Developments in Nuclear In- 
strumentation at the Knolls Atomic 
Power Laboratory, by Richard S. Stone, 
ASME Pap. 57-NESC-70; Amer. Inst. 
Elect. Engrs., Pap. 57-399, March 1957, 
pp. 

A Wire Activation Technique for Reactor 
Flux Profile Measurements, by Alton | 
KXlickman and Francis R. DeFaleo, ASME 
Pap. 57 NESC-40; Amer. Inst. Elect 
Engrs. Pap. 57-405, March 1957, 4 pp. 

Designing Thermocouples for Response 
Rates, by R. J. Moffat, ASME Pap. 57 
GTP-8, March 1957, 8 pp. 

Methods of Measuring the Thermal 
Conductivity of Solids, by G. L. Cooper, 
Gt. Brit., Atomic Energy Res. Estab., Rep. 
Inf/Bib. 104, 1956, 10 pp. 

Radar Video Integrator Study, by Roy 
W. Hendrick (Internal Research Project 
81-011), Cornell Aeron. Lab., Rep. CAL.- 
78, Jan. 1957, 51 pp. 

Bourdon Elements Require Special 
Tubing, by Charles C. Mathews, /nstru- 
ments & Automation, vol. 30, May 1957, 
pp. 891-893. 

Precision Automatic Manometer 
Reader, by J. Farquharson and H. A. 
Kermicle, Rev. Sct. Instrum., vol. 28, May 
1957, pp. 324-325. 

Thermistor Thermometer Bridge: Line- 
arity and Sensitivity for a Range of Tem- 
peratures, by Kenneth S. Cole, Rev. Sci. 
Instrum., vol. 28, May 1957, pp. 326-328. 

Theory of a New Apparatus for Deter- 
mining the Thermal Conductivities of 
Metals, by S. T. Hsu, Rev. Sci. Instrum., 
vol. 28, May 1957, pp. 333-336. 

X-Ray Powder Diffraction Assembly for 
Studies at Elevated Temperatures and 
High Gas Pressures, by Edward J. Goon, 
John T. Mason and Thomas R. P. Gibb 
Jr., Rev. Sci. Instrum., vol. 28, May 1957, 
pp. 342-344. 

Apparatus for the Measurement of 
Vapor Pressure and Ignition of Decom- 
position Temperature at Pressures Above 
Atmospheric, by A. Greenville Whittaker, 
Rev. Sci. Instrum., vol. 28, May 1957, pp. 
360-364. 
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pressurization equipment proiccts vital 

© gear. A continual program of re- 

nd development creates customized 

ation units that keep the performance 

systems unaffected by altitude and 

ent conditions. Custom units that meet mil- 

specifications help to solve your problems 
recommending electronic components. 

n yeu have a challenging problem to pre- 

pressure, or heat, or moisture, or dust 

_ affecting electronic performance, come to 


East n for complete and creative cngincering 1500 SERIES 
PRESSURIZATION UNIT 


zero to over 70,000 feet at temperatures 

m — 65°F to +160°F. Delivery: 0-3600 cu. 

/min. free delivery, Dischazge Pressure : 0-60 

Standard sub-assemblies and components 

: prmaily are used to create a custom-made de- 

_ m to fit your exact needs. Units may consist 

100 SERIES an air pump and motor assembly, pressure 
PRESSURIZATION UNIT 


INDUSTRIES, INC. 
100 Skiff St., Hamden 14, Conn. 


West Coast Office: 1608 Centinela Avenue 
Inglewood 3, California — Phone ORegon 8-3958 
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Development for the NATIV Servo Sys- W. 
tem, by W. D. Mullins Jr. and K. R. ph 
Jackson, North American Aviation Inc., ] 
N L Rep. AL162, March 1947, 15 pp. (De- Ve 
classified from Confidential by authority Me 
of North American Aviation, Inc., letter ] 
54D-6322, citing WADC letter WCSG, Ra 
Get in now —at the beginning Ra 
/ a 8 Survey of Field Test Range Instrumen- 
of the new era in missiles! tation Requirements for NATIV Firings, 
} by Donald H. Jacobs, Warren Kott and M 
When you join Telecomputing’s ATTRACTIVE SALARIES Seymour Scholnick, North A merican Aviu- Ce 
Engineering Services Division, PROFIT SHARING a tion, Inc., Rep. AL220, June 1947, 48 pp. 
you will be given full scope to (Declassified from Confidential by author- I 
allow you to grow... your talents of ’ Hi 
will be used to the fullest... RELOCATION PAY WCSG, 8/6/54) 
recognition and rewards se be ACCREDITED EDUCATION The NAA Telemetering System, by Br 
yours as a matter of —— GROUP INSURANCE Donald H. Jacobs and Claude W. Chap- 4 
Engineering services is a mem- meetin man, North American Aviation, Inc., Rep. I 
ber of an integrated five-com- . AL-371, Nov. 1947, 59 pp. (Declassified Aci 
pany missile systems corporation A NEW LIFE IN NEW MEXICO'S from Confidential by authority of North lius 
which designs and manufactures FABULOUS “LAND OF ENCHANTMENT" American Aviation, Inc., letter 54D-6233;, pul 
its own data-processing equip- . 2» | citing WADC letter W CSG, 8/6/54.) 8 p 
The NAA Time Standard System, | 
MOUNTAIN SKIING AND DESERT Donald H. Jacobs and Michael May 
RESORTS WITHIN 30 MINUTES! North American Aviation, Inc., Rez Cor 
are engaged in the reduction of 4) AT OL7 
a ave AL372, Nov. 1947, 35 pp. (Declassifiec Cal 
large amounts of flight test data ae from Confidential by authority of North Bee 
being generated by the daily mi<- A WONDERFUL PLACE ' American Aviation, Inc., Letter 54D-6233, anc 
sile firings on the integrated Hol- TO MAKE YOUR HOME — citing WADC letter WCSG, 8/6/54.) ISC 
loman-White Sands range. GRAND COUNTRY TO RAISE KIDS! The NAA Plotting Cameras, by Donald I 
H. Jacobs, North American Aviation Th 
. . Inc., Rep. AL378, Dec. 1947, 24 pp. (De Dir 
Send resume to Director of Technical Personnel boy NE 
of North American Aviation, Inc., lette En 
TELECOMPUTING CORPORATION 54D-6233, citing WADC letter WCSG, Ma 
. e ee 8/6/54.) 
Engineering Services Division (NAA Radar Recording Systems, by 
onald H. Jacobs, North American Avia- 
BOX 447 » HOLLOMAN AIR FORCE BASE » NEW MEXICO tion, Inc., Rep. AL380, Rev. March 1948, D. 
16 pp. (Declassified from Confidential by 103 
| authority of North American Aviation, and 
Inc., letter 54D-6233, citing WADC lette: 14: 
WCSG, Aug. 6, 1954.) 1 
NAA Radar Modifications, by Donald H Ura 
Jacobs, North American Aviation, /nc., Sall 
MISSILE TEST DATA Rep. AL379, Rev. March 1948, 9 pp. AS 
: (Declassified from Confidential by author- 8 p 
TRANSCRIPTION AND ANALYSIS — ity of North American Aviation, Inc., 2nd 
: letter 54D-6233, citing WADC letter E 
The transcription, analysis, and evaluation of ballistic missile test data WCSG, 8/6/54.) Ste: 
is one of the important aspects of R-W’s system engineering and tech= Semiannual Progress Report for the Ene 
nical direction responsibility for the Air Force Ballistic Missile Program: Period Ending June 30, 1955, Nat. Bur. 195 
Utilizing the facilities of R-W’s Data Reduction Center, test data are A meget Off. Basic Instrumentation, Rep. T 
transcribed and reviewed to determine performance characteristics and Ves 
A Simply Contour Gage, by Con 
confidence levels. 3 W. W. Oshe 1 Willis C. st VA 
The following positions in this rapidly expanding field are currently 
VORD, Rep. 5436 (NOTS 1704), April Pip 
open on the R-W staff: 1957, 133 pp. — 7 pis 
Senior analyst to direct mathematicians, engineers, and technical aides 4 Res 
in missile test analysis projects and to participate in overall evaluation a 28 | 
of missile systems performance. The test information originates from coalte = T 
optical, radar, and telemetry instrumentation. Particular emphasis 1S Astrophysics, Aerophy SICS = 
laced upon electronic analysis of vibration data. oo 
Reduction Methods for Photographic VA 
Engineer with wide experience in telemetry operations. The position Meteor Trails, by Fred L. Whipple and NA 
olves planning and directing the conversion of telemetered data into Luigi G. Jacchia, Smithsonian Contribs. T 
A forms appropriate for the analysis of weapon systems performance. Expe= Bt Astrophys., vol. 1, no. 2, 1957, pp. 183 Stre 
ience in the design, maintenance, and operation of telemetry data ie 206. a 
transcription and electronic vibration analysis equipment is desirable: i The Method of Reduction of Short- > 
Trail Meteors, by Gerald S. Hawkins, 
Engineer experienced in the utilization of missile test data transcription = Smithsonian Contribs. Astrophys., vol. 1, 
equipment. This position involves the planning of equipment schedules; no. 2, 1957, pp. 207-214. Ele 
with various research and industrial organizations, for the conversion Of A Rapid Graphical Method of Meteor by 
_ missile test data. A thorough knowledge of the capabilities and limitations ‘ Trial Reduction, by Richard E. McCrosky, Aér 
of telemetry transcription and electronic vibration analysis equipment 5 Smithsonian Contribs. Astrophys., vol. 1, pp- 
s desirable. ie no. 2, 1957, pp. 215-224. 
ah, A Reduction Method for the Motions of At 
Inquiries should be addressed to: Mr. J. H. Armitage Persistent Meteor Trains, by Allan F. 
a Cook and Robert F. Hughes, Smithsonian T 
Contribs. Astrophys., vol. 1, no. 2, 1957, pp. clea 
The Ramo-Wooldridge Corporat = 
00 rl ge orpora for Study of 
ts from hotographic eteor 
5730 ARB ITAE STREET: ANGE! 45. CAL Trails, by Fred L. W 195’ 
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W. Wright, Smithsonian Contribs. Astro- 
phys., vol. 1, no. 2, 1957, pp. 239-243. 

IGY World Warning Agency, Tech. 
News Bull. (Nat. Bur. Standards), vol. 41, 
May 1957, pp. 65-66. 

Meteor Burst Extends VHF Radio 
Range, by Philip J. Klass, Aviation Week, 
vol. 66, June 17, 1957, pp. 96-97, 99-101. 


Materials of 
Construction 


Influence of Crucible Materials on 
High-Temperature Properties of Vacuum- 
Melted Nickel-Chromium-Cobalt Alloy, 
ny R. F. Decker, John P. Rowe and J. W. 
Freeman, NACA TN 4049, June 1957, 
34 pp. 

Inhibition of Titanium in Fuming Nitric 
Acid, by John B. Rittenhouse and Corne- 
ius A. Papp, Calif. Inst. Tech., Jet Pro- 
pulsion Lab., Mem. 20-140, -March 1957, 
8 pp. 

Tabulation, Bibliography and Structure 
of Binary Intermetallic Compounds, II: 
Compounds of Beryllium, Magnesium and 
Calcium, by K. A. Gschneidner, D. J. 
Beernstsen, R. W. Vest, J. A. Kingston 
ind J. F. Smith, Atomic Energy Comm., 
ISC-812, Dec. 1956, 36 pp. 

Internal Friction and Shear Modulus of 
Thorium at High Temperatures, by C. FE. 
Dixon and H. Hori, ASME, Pap. 57- 
NESC-6, March 1957, 4 pp. (Nuclear 
Engng. and Sci. Conf., 2nd., Philadelphia, 
March 11-14, 1957). 

Fabricated Pressure Piping as Related 
to Nuclear Applications, by J. J. Murphy, 
C. R. Soderberg Jr., H. S. Blumberg and 
D. B. Rossheim, ASME, Pap. 57-NESC- 
103, March 1957, 24 pp. (Nuclear Engng. 
and Sci. Conf., 2nd., Philadelphia, March 
11-14, 1957). 

The Effect of Oxygen on Zirconium 
Uranium Epsilon Phase Alloys, by H. A. 
Saller, F. A. Rough and A. A. Bauer, 

ASME, ie 57-NESC-20, March 1957, 
pp. Nuclear Engng. and Sci. Conf., 
hig Philadelphia, March 11-14, 1957). 

Bibliography on the Brazing of Mild 
Steel, by G. L. Cooper, Gt. Brit., Atomic 
Energy Res. Estab., Rep. Inf/Bib 103, 
1956, 4 pp. 

Thermal Stresses in HRT Pressure 
Vessel, by R. E. Aven, Atomic Energy 
Comm., CF-55-2-128, Feb. 1955, 17 pp. 

MTR Hot Cell Tests of High Pressure 
Piping, by M. S. Robinson, F. W. Rider 
and M. H. Bartz, Atomic Energy Comm. 
Res. and Dev. Rep. 1DO-16294, Dec. 1956, 
28 pp. 

Two Factors Influencing Temperature 
Distributions and Thermal Stresses in 
Structures, by William A. Brooks, Jr., 
George E. Griffith and H. Kurt Strass, 
NACA TN 4052, June 1957, 13 pp. 

The Combinations of Thermal and Load 
Stresses for the Onset of Permanent 
Buckling in Plates, by George W. Zender 
and Richard A. Pride, NACA TN 4053, 
June 1957, 10 pp. 

Trends in the Use of Materials in the 
Electronics Industry in the United States, 
by M. Bousquet, Fusées et Recherche 
Aéronautique, vol. 1, no. 2, Oct. 1956, 
pp. 129-130 (in French). ty” 
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The Safety Aspects of the Use of Nu- 
clear Energy, by C. R. Russell, SAEZ, 
Prepr. 4, Jan. 1957, 7 pp., 7 fig. 

Unclassified Engineering Materials List, 
Atomic Energy Comm., TID 4100, Jan. 
1957. 
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7500 SERIES DC 
ELAPSED TIME 
INDICATOR 


Send for Bulletins: AWH ET600 
Basic Elapsed Time Indicators 
7500, 12500, 24200 Series. 


Now the Company offers a com- 
plete line of timing motors and devices 
to record the operating time of any 
electrical or electronic equipment. 
Compact, minimum weight, each unit 
has five digits. They can be used to 
provide daily running time plus a 
total running record, eliminating 
estimating or manual totalizing. They 
increase efficiency and reduce costs, 
since they meme the quickest wa 
to know when to lubricate, pote 
or replace equipment. They give on- 
the-spot of idle waiting 
periods for machinery. Used with 
electronic equipment, they tell when 
to replace tubes, preventing costly 
failures during essential operating 
periods. 

In both AC and DC, continuous or 
manual reset for 50 or 60 up to 400 
cycle line frequency. Will measure 
hours or on down to 10ths of seconds. 
All models can be supplied with Radio 
Interference Filtering to meet MIL- 
1-6181B. 


A. W. HAYDON COMPANY'S 


ELAPSED TIME 


the 


INDUSTRIAL, 


AIRCRAFT 
and 


7500 SERIES RESET TYPE 
ELAPSED TIME INDICATOR 


Send for Bulietins: AWH ET601 
Reset Type Elapsed Time Indica- 
tors 7500, 12500, 24200 Series. 


Preferred Where Performance 
is Paramount. 


; Shown at the left is the new A. W. Haydon Co. catalog ; 
describing the complete line of timing motors and 
. devices; if you haven't received your copy write for it | 
on your company letterhead. 


Company 


248 NORTH ELM STREET, WATERBURY 20, CONNECTICUT 
Design an! Manufacture of Electro-Mechanical Timing Devices 
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Size, weight, performance~a key con- 
sideration in helium pressurization 
projects. Turn to R-F Aeronautical — 
a prime source for helium control 
components. Engineering experience, 
complete helium testing facilities... 
combined to meet your pressurization 
requirements — faster and better. 


TYPICAL — this missile borne Helium 
Pressure Regulator with integral shut-off. 
INLET PRESSURE: 325425 psig, NOMINAL 
OUTLET PRESSURE: 20+0.5 psia (or psig) 
FLOW RATE: 2000 scfm, OPERATING TEMP: 
AMBIENT: —65° to 165°F, GAS: —80° to 
350°F, INITIAL RESPONSE TIME: 20sec. 
WEIGHT: 2.75 Ibs. max. DIMENSIONS: 674” 
x3%”"x5-15/16”, no flutter, chatter, over- 
s shoot, or instability. 


At R-F Aeronautical —a complete line 
helium pressure regulators... 
3 shelf and to specification ... smaller, 
lighter, more reliable. Write for 
technical datey 
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